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I.  INTRODUCTION 


Estimating  the  blast  loading  on  various  targets  is  an  ongoing  effort  at 
the  US  Army  Ballistic  Research  Laboratory  (BRL).  Such  loading  may  be 
estimated  through  tests  in  shock  tubes,  through  high-explosives  (HE)  testing 
in  the  field,  or  through  computation  with  hydrocodes.  One  such  hydrocode, 
which  has  been  used  with  success  at  BRL  and  elsewhere,  is  an  airblast  version 
1  2 

of  the  HULL  *  hydrodynamic  computer  code.  This  code  solves  the  inviscid 
Euler  equations  in  two  dimensions  (2-D)  or  three  dimensions  (3-D)  using  an 
explicit  time-stepping  method.  The  2-D  version  has  a  number  of  options  that 
are  not  available  in  3-D,  but  the  main  disadvantages  of  the  3-D  version  are 
the  much  longer  run  time  and  the  memory  requirement.  The  mesh  spacing  for 
3-D  runs  is  too  frequently  set  much  coarser  than  desirable  in  order  to  fit 
the  storage  for  the  required  computational  field  into  the  computer  or  to 
complete  the  run  in  a  reasonable  time. 

For  some  studies  (simulation  of  an  HE  field  test,  for  example)  it 
would  be  advantageous  to  use  the  2-D  cylindrical  symmetry  version  of  HULL 
from  initiation  to  near  impingement  of  the  blast  wave  on  a  target  position, 
then  continue  the  calculation  with  the  3-D  version  to  study  loading  on 
a  3-D  target  (2-D  cylindrical  HULL  runs  can  be  preceded  by,  and  initiated 
from,  one-dimensional  runs  of  the  SAP  program).  The  2-D  cylindrical  HULL 
could  be  used  as  long  as  the  problem  was  cylindrically  symmetric.  It  would 
be  much  faster  than  the  3-D  version,  even  with  a  finer  mesh,  and  could 
incorporate  some  useful  options  that  are  not  present  in  3-D  such  as  dust  or 
multiple  materials.  However,  the  2-D  HULL  would  not  be  appropriate  for 
examining  the  loading  on  a  3-D  target  such  as  a  building  or  a  vehicle. 

The  hydrodynamic  data  in  a  portion  of  the  2-D  cylindrical  HULL  computational 
field  (actually  three-dimensional)  could  be  transferred  to  a  smaller  3-D  HULL 
computational  field  containing  the  target,  and  the  computation  continued  with 
the  3-D  HULL.  (Some  optional  properties,  multiple  materials  for  example,  do 
net  exist  in  3-D  and  cannot  be  transferred.) 

Such  a  procedure  was  at  least  partially  available  in  our  HULL.  The  grid 
generation  program,  called  KEEL,  will  transfer  results  from  a  2-D  cylindrical 
problem  to  a  new  2-D  or  3-D  problem.  The  user  would  have  had  to  assign  some 
existing,  but  unrealistic,  boundary  conditions.  The  results  from  the  new  run 
would  be  meaningful  until  the  effect  of  whatever  boundary  conditions  were 
imposed  reached  the  target  area.  The  most  likely  boundary  conditions  would  be 
transmissive  or  reflective.  If  realistic  boundary  conditions  from  the  2-D 
run  could  be  imposed,  the  results  would  be  useful  until  the  reflected  waves 
went  from  the  target  to  a  nonreflective  boundary  and  back  to  the  target. 

Using  such  boundary  conditions  would  approximately  double  the  error  free 
time,  from  program  initiation  to  the  arrival  of  false  signals  at  the  target, 
for  a  given  3-D  computational  grid.  To  double  the  uncontaminated  run  time 
using  the  wrong  boundaries,  the  distance  from  the  target  sides  to  the  edges 
of  the  computational  region  would  have  to  be  doubled;  for  a  relatively  small 
3-D  target,  the  computational  region  would  have  to  be  about  13  times  as  large. 

This  report  documents  the  coding  for  this  new  type  of  boundary  for 
HULL,  and  presents  the  results  of  some  limited  tests.  Section  IT  summarizes 
the  older  HULL  boundaries;  the  new  boundary  type  is  discussed  in  Section 
III.  The  program  to  produce  boundary  input  for  a  3-D  HULL  run  from  the 
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restart  output  from  a  2.-X)  cylindrical  HULL  run  is  discussed  in  Section  IV. 

The  results  from  a  series  of  test  runs  is  discussed  in  Section  V.  Appendix  A 
contains  a  tabulation  of  a  sample  runstream  and  the  HULLUP  program  which  will 
produce  boundary  input  for  a  3-D  run.  Appendix  B  includes  listings  of  a 
runstream  to  use  KEEL  for  a  o-D  restart,  a  runstream  for  the  corresponding 
HULL  run,  the  general  UPDATE  corrections  to  HULL,  and  the  HULL  changes  for 
the  new  boundary  conditions. 


II.  HULL  COMPUTATIONAL  GRID  BOUNDARY  CONDITIONS 

In  order  to  understand  the  computational  grid  boundaries  used  in  HULL, 
one  needs  an  overview  of  the  HULL  method  for  solving  the  governing 
equations.  The  basic  equations  used  in  HULL  to  describe  inviseid  fluid  flow 
are: 

dP  ,  -k 

dt  +  9  (V'U)  =  0  (1 


P 


d  u 
dt 


+  Vp  =  -p  ~g 


(2. 


+  V  •  (Pu)  =  -  p\l- g* 


(3: 


P=  F(P,  I ) 


(4 


£=*!+%(  M) 


(5; 


It  +  pyeff  (V'“)  =  0 


wnere , 

v ■ =  the  divergence  operator 

V  s  the  gradient  operator 

p  =  material  density  (g/em^) 
■  2 

p  =  pressure  (dynes/cm  ) 


u  =  the  vector  fluid  velocity  (cm/sec) 
1  =  specific  internal  energy  (ergs/g) 

E  =  specific  total  energy  (ergs/g) 


g 


*=  the  vector  acceleration  of  gravity  (cm/sec  ) 
t  =  time  (sec) 


Y  „  =  1  +  p/  I  is  the  "effective  gamma" 
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Equations  (l),  (2),  and  (3)  are  the  conservation  equations  for  mass, 
momentum,  and  energy,  respectively.  Equation  {h)  is  the  equation  of  state 

(for  this  study,  a  gamma  law  equation  of  state,  p  =  (Y  -  l)pl,  was  assumed 

with  Y  =  l.U).  Equation  (6)  can  be  derived  by  substituting  equations  (5)  and 
(3)  into  equation  (2)  and  assuming  y  =yeff  is  constant. 

The  solution  is  found  by  two  phases  for  each  time  step.  The  first 
phase  is  a  two  step  Lagrangian  solution  in  the  manner  of  Lax-Wendroff . 

The  second  phase  is  a  transport  phase  which  is  also  carried  through  in 
two  steps.  (A  description  of  one  cycle  of  the  HULL  solution  at  a  cell 
for  a  2-D  cylindrically  symmetric  case  is  contained  in  the  SAIL  comments 
at  the  beginning  of  the  HULL  program.  In  fact,  there  are  two  descriptions; 
the  original  HULL  authors  were  in  mild  disagreement.  These  comments  have 
equal  signs  in  column  1  and  are  removed  by  the  preprocessor  POST.  They  are 
in  the  output  of  UPDATE  or  of  program  LIST  which  is  part  of  the  HULL  file. 

A  description  of  the  phase  1  solution  for  3-D  HULL,  with  comments,  is 

contained  in  Section  2  of  Reference  3.) 

At  the  beginning  of  phase  1  (at  time  t)  the  velocity  components,  the 
specific  internal  energy,  and  the  mass  are  known  for  every  interior  cell  of 
the  computational  field.  Routines  and  data  are  available  to  compute  other 
values  such  as  volume,  density,  and  pressure. 

When  the  phase  1  solution  is  to  be  advanced  for  a  particular  cell, 
it  is  assumed  that  the  pressure  and  normal  velocity  components  are  known 
for  each  low  index  side  of  the  cell  at  time  t  +4  t/2.  (It  is  convenient  to 
think  of  these  intermediate  time  values  being  on  the  cell  boundaries.  They 
are  more  closely  Lagrange  solutions  for  points  that  were  between  cell  centers 
at  time  t.)  From  the  values  in  the  cell  and  in  the  adjacent  higher  index 
cells,  the  pressure  and  normal  velocity  components  for  time  t  +4  t/2  are 
computed  at  the  high  index  sides  of  the  cell.  This  completes  step  1  of  phase 
1  for  the  cell.  The  second  step  of  phase  1  uses  these  pressures  and 
velocities  "at  the  cell  boundaries"  to  compute  the  Lagrange  solution  for 
specific  energy  and  velocity  rt  time  t  +4t  for  the  old  cell  center  point. 

For  most  cells  the  values  at  the  low  index  sides  of  the  cell  for  time 
t  +At/2  are  retained  from  the  previous  calculations.  If  the  cell  is  a  low 
index  boundary  cell  (l=l,  J=l,  or  K=l)  these  boundary  values  must  be 
supplied.  The  method  of  supplying  these  values  differs  with  the  condition, 
the  boundary,  the  dimension,  and  the  programmer.  Four  fairly  distinct  methods 
are  used  at  low  index  boundaries: 

1.  A  virtual  external  cell  with  the  desired  properties,  at  time  t. , 
is  assumed  and  the  governing  equations  solved  for  pressure  and 
normal  velocity  at  time  t  +At/2  on  the  boundary. 


2.  The  pressure  and  normal  velocity  on  the  boundary  at  time 
t  +  At/2  are  assigned.  (Some  inconsistency  in  time  was 
discovered  during  this  study.  Changes  to  make  them  consistent 
were  added  to  the  correction  file. ) 

3.  Hydro  values  are  ass.  ined  to  the  boundary  cell  and  computation 
for  that  cell  is  skipped  except  for  computing  the  velocity  and 
pressure  at  t  +  At/2  on  the  high  index  side.  This  is  used 
only  with  2-D  at  the  bottom  boundary.  The  bottom  layer  of 
cells  is  considered  to  be  inside  the  computation  region,  hence 
the  mass  and  energy  changes  in  the  boundary  row  of  cells  is 
added  to  the  system  mass  and  energy. 

k.  Hydro  values  are  assigned  to  the  cell  and  approximate  boundary 

values  are  imposed  on  the  low  index  side.  The  governing  equations 
are  solved  for  the  boundary  cell,  but  the  results  are  in  error 
because  of  the  approximate  boundary  values.  The  error  is  not 
propagated  since  new  values  are  assigned  to  the  cell  for  the  next 
step.  Again,  the  program  must  account  for  system  energy  and  mass. 

Boundary  values  at  high  index  boundaries  are  all  imposed  in  the  same 
manner:  The  cell  with  the  highest  index  (I=IMAX,  J=JMAX,  or  K=KMAX)  is  a 
virtual  external  cell;  appropriate  values  are  assigned  to  these  virtual  cells 
for  the  desired  boundary  condition  and  these  values  are  used  to  compute  the 
pressure  and  normal  velocity  component  at  time  t  +At/2  on  the  high  index 
side  of  the  adjacent  internal  cell. 

The  controls  for  the  phase  2  calculations,  the  transport  phase,  are 
similar.  The  first  step  in  phase  2  computes  the  fluxing  of  mass,  momenta, 
and  energy  through  the  cell  boundaries.  The  second  step  computes  the  result 
of  this  fluxing.  Initially,  for  any  internal  cell,  the  Lagrangian  solution 
for  velocity  and  specific  energy  (at  the  point  that  was  the  cell  center  at 
time  t)  are  known  at  time  t  +  At;  the  mass  (hence  density)  of  the  cell  is 
known  for  time  t.  The  average  of  the  normal  velocity  components  of  two 
contiguous  cells  is  used  as  the  rate  of  mass  transport  between  them.  From 
this,  the  mass  flux,  the  momentum  fluxes,  and  the  energy  flux  between  the  two 
cells  are  computed. 

Just  as  for  phasa  1,  the  values  on  the  low  index  sides  of  the  cell 
are  assumed  to  be  known,  either  from  computation  for  the  cell  on  that  side  or 
from  boundary  condition  assignment.  Computing  the  fluxing  for  the  high  index 
sides  of  the  cell  completes  step  1  of  phase  2.  Step  2  of  phase  2  consists  of 
summing  all  the  fluxes  and  computing  the  resulting  mass,  velocity  components, 
and  specific  internal  energy  for  the  cell  at  time  t  +At. 

Similar  boundary  control  methods  to  those  described  for  the  first 
phase  are  used  for  phase  o.  in  addition,  any  transfer  of  mass  or  energy 
through  the  boundaries  nus4  be  accounted  for.  (In  2-D,  some  of  the  options 
such  as  mult i^at  >ria !  or  radiation  require  additional  computation.) 

The  2-D  and  S-D  codes  ira  really  separate  programs  although  they  share 
some  coding.  They  .*v»n  iit'fer  slightly  in  the  solution  of  the  governing 
equations.  ^'or  example,  2-D  simply  averages  for  values  between  cell 
centers,  4  he  interpol at  es  for  some  of  them.) 


Reading  the  comment  lines  with  the  HULL  coding  could  lead  one  to  believe 
there  were  eight  types  of  boundary  conditions  available  at  every  boundary. 
This  is  not  true.  Of  the  eight  options  listed,  only  one  type,  reflective,  is 
available  at  all  the  boundaries.  Table  1  is  a  summary  of  the  author's 
assessment  of  the  availability  and  reliability  of  the  boundary  conditions  in 
our  alrblast  HULL.  The  eight  types  numbered  0  through  7  are  called 
reflective,  transmissive,  rezone,  square  wave,  LAMB,  SAP,  HULL,  and  oblique 
square  wave,  respectively. 

The  next  section  will  describe  the  new  imposed-solution  boundaries 
developed  at  the  BRL. 

Table  1.  Summary  of  HULL  Boundaries. 


BOUNDARY  TYPE 


BOUNDARY 

0 

1 

2 

3 

k 

5 

6 

7 

AFT 

3-D 

YES 

A,C 

A,C 

A 

A 

— 

— 

YES 

BOTTOM 

S-D 

YES 

A,C 

A.C 

- 

- 

- 

- 

A 

BOTTOM 

2-D 

YES 

A,C 

A,C 

B 

B 

A 

- 

YES 

LEFT 

J-D 

YES 

YES 

A,C 

YES 

YES 

- 

- 

YES 

LEFT 

2-D 

YES 

A,C 

A,  C 

YES 

A 

A 

A 

YES 

FORE 

3-D 

YES 

YES 

- 

- 

- 

- 

- 

- 

TOP 

j-D 

YES 

YES 

- 

-  ■ 

- 

- 

- 

- 

TOP 

2-D 

YES 

YES 

YES 

- 

- 

- 

- 

- 

RIGHT 

o-D 

YES 

YES 

- 

- 

- 

- 

- 

- 

RIGHT 

2-D 

YES 

YES 

YES 

— 

- 

- 

Has  been  used  locally.  Looks  good. 

No  implementation. 

Has  not  been  used  locally.  Appears  to  be  fully  implemented.  May  have 
errors. 

Only  partially  implemented  or  has  known  errors. 

Types  1  and  2  alike  at  each  low  index  boundary.  Believe  they  are 
intended  to  remain  unchanged  ambient. 


(YES) 

(-) 

(A) 

IB) 

(O 
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III.  BOUNDARY  TYPE  9,  IMPOSED-SOLUTION  BOUNDARY 


The  primary  motivation  for  this  project  was  to  transform  a  section 
of  a  2-D  cylindrical  symmetry  HULL  run  into  a  3-D  Cartesian  HULL  run. 
Presumably,  the  2-D  HULL  program  would  be  run  with  a  blast  wave  expanding 
from  some  central  source.  The  3-D  run  would  be  for  some  subsection  of  the 
cylindrical  space  with  a  target  inside.  The  3-D  run  could  be  initiated  just 
before  the  blast  wave  reached  the  target  and  would  continue  with  input  from 
the  2-D  solution  through  the  boundaries  of  the  3-D  grid. 

The  KEEL  program  has  coding  for  copying  a  portion  of  a  2-D  cylindrical 
HULL  problem  into  a  3-D  space  but  the  boundary  conditions  from  the  2-D  run 
could  not  be  satisfactorily  fed  in  after  the  3-D  run  was  initiated.  The 
result  from  such  a  restart  would  be  good  until  a  signal  from  whatever 
boundaries  were  used  reached  the  region  of  interest.  With  the  addition  of 
input  boundaries,  with  data  from  the  2-D  donor  run,  the  results  should  be  good 
until  reflections  from  the  target  reach  a  boundary  and  the  resulting  erroneous 
signal  reaches  the  region  of  interest.  Imposing  a  solution  on  the  boundaries 
of  the  3-D  grid  would  essentially  double  the  real  time  that  could  be  simulated 
by  the  same  grid  with  simple  transmissive  boundaries.  For  comparison,  the 
same  time  doubling  of  simulated  time  could  be  achieved  by  doubling  the 
distances  from  the  target  to  the  simple  transmissive  boundaries  in  the  3-D 
problem,  but  this  could  increase  the  space  and  computer  time  required  for  such 
a  run  by  a  factor  of  eight.  Thus,  this  new  capability  offers  the  opportunity 
to  run  problems  which  might  otherwise  be  too  large  or  expensive. 

In  keeping  with  this  scenario,  the  coding  for  these  new  boundaries  allows 
imposed  values  on  any  side.  (Previous  programming  had  not  added  boundary 
conditions  at  the  high  index  boundaries  even  though  the  existence  of  the  layer 
of  external,  virtual  cells  make  this  relatively  easy.)  These  "imposed- 
solution"  boundaries  will  be  called  B0UND9  boundaries  in  this  report.'  They 
will  be  activated  in  HULL  as  type  9.  (Type  8  has  been  assigned  at  the  BRL  to 
an  undocumented,  partially  successful  exhaust  type  boundary  that  may  be 
reactivated. ) 

After  some  preliminary  study,  it  was  decided  that  the  simplest  way  to 
introduce  these  new  HULL  boundary  conditions  into  the  existing  code  was  to 
supply  input  values  on  a  plane  for  each  boundary,  with  a  new  subroutine  in 
HULL  for  each  of  these  boundaries.  These  subroutines  check  that  the  point  for 
which  data  is  requested  is  on  the  plane.  Then  they  interpolate  in  time,  and 
in  space  on  the  plane,  for  five  hydro  values:  three  velocity  components, 
specific  internal  energy,  and  density.  At  the  low  index  sides  these  planes 
are  on  the  boundaries:  the  "left"  boundary  at  XO,  the  "aft"  boundary  at  YO, 
and  the  "bottom"  boundary  at  ZO.  For  these  boundaries,  the  coding  computes 
and  stores  the  pressure  and  the  normal  velocity  component  at  time  t  +  t/2, 

like  method  (2)  in  the  previous  section.  At  the  high  index  boundaries,  the 
input  planes  for  the  imposed-solution  values  go  through  the  center  of  the 
external  virtual  cells  on  that  side:  the  "right"  boundary  for  maximum  X,  the 
"fore"  boundary  for  maximum  Y,  and  the  "top"  boundary  for  maximum  Z.  The 
coding  inserts  the  "imposed"  hydro  data  into  the  external  cell  (mass  is  stored 
instead  of  density). 

These  new  B0UND9  boundaries  are  designed  to  continue  a  2-D  HULL  run  in 
3-D,  but  any  sort  of  values  could  be  imposed  on  a  side  if  a  file  of  data  in 


'  h'  ,4f  .41  , 


the  proper  form  were  prepared  for  that  side.  (This  assumes  that  the  boundary 
values  are  compatible  with  whatever  is  inside  the  computation  region.) 

The  original  plan  was  to  supply  two  planes  of  hydro  data  at  successive 
times  for  each  B0UND9  side.  An  evenly  spaced  grid  with  extreme  values  at  the 
B0UND9  planes  was  also  imposed.  This  worked  well  for  small  3-D  HULL  grids 
with  evenly  spaced  meshes.  The  examples  cited  later  were  run  with  this 
coding,  but  some  memory  requirements  were  excessive  even  with  a  coarser  mesh 
on  the  input  data  planes  than  was  wanted.  Therefore,  a  new  code  was  prepared 
which  presents  the  imposed-boundary  data  to  the  3-D  HULL  one  row  at  a  time. 

To  avoid  excessive  searching  for  data,  the  hydro  values  for  two  successive 
times  are  stored  together  for  each  row. 

A  description  of  this  B0UND9  input  for  HULL  is  included  in  Appendix  B 
along  with  a  listing  of  the  B0UND9  changes  for  HULL  and  listings  of  the 
runst reams  for  KEEL  and  HULL  for  one  problem. 

IV.  HULLUP 

The  program  to  prepare  the  files  for  the  3-D  HULL  B0UND9  input  boundaries 
from  a  2-D  donor  run  is  a  FORTRAN  5  program  called  HULLUP.  This  program  is 
tabulated  in  Appendix  A.  It  includes  COMMENT  lines  that  hopefully  adequately 
describe  the  input  and  output.  The  main  output  is  the  files  of  data  for  the 
B0UND9  input  to  HULL  described  in  Appendix  B. 

The  original  intent  was  to  allow  for  2-D  and  3-D  input  and  output  at  some 
later  date.  There  are  some  remnants  of  this  remaining,  but  the  coding  is  now 
strictly  for  producing  HULL  input  on  3-D  boundary  planes  from  a  2-D 
cylindrical  HULL  run. 

The  user  must  plan  the  3-D  run  thoroughly  before  knowing  the  input  for 
HULLUP.  If  there  is  a  target,  it  is  one  or  more  orthogonal  parallelepipeds 
(boxes).  If  a  target  is  to  have  flat  sides,  the  3-D  grid  must  have  boundaries 
parallel  to  the  sides  of  the  target.  These  boundaries,  except  for  reflective 
boundaries,  must  be  far  enough  from  the  target  area  to  allow  HULL  to  run 
without  sending  a  false  signal  to  the  target  area  during  the  time  of  interest. 
(For  B0UND9  boundaries,  the  signal  from  the  boundary  will  be  approximately 
correct  until  reflection  from  the  target  reaches  them. )  Any  B0UND9  boundary 
plane  must  be  entirely  inside  the  cylindrical  region  of  the  2-D  donor.  The 
'angle  at  which  the  shock  front  strikes  a  target  is  determined  by  the 
positioning  of  the  target.  The  transfer  from  the  2-D  donor  space,  which  HULL 
denotes  as  (X,Y)  is  clearer  if  one  thinks  of  the  2-D  cylindrical  space  as 
(R,Z),  or  the  3-D  cylindrical  space  (R,0,Z)  at  arbitrary  9  . 

Points  in  the  Cartesian  3-D  space  are  defined  as  X  =  R  cos(0), 

Y  =  R  sin(0),  and  Z  =  Z  (no  rotation  of  Z  plane).  The  user  does  not  need  to 
know  9  .  To  set  the  boundaries  of  the  3-D  Cartesian  space  the  user  must 
envision,  or  draw,  the  X  and  Y  boundaries  of  the  3-D  computational  field  on  a 
constant  Z  plane.  A  line  through  the  origin  parallel  to  the  X  sides  of  the 

computational  field  is  6=  0,  and  9  =  90°  is  parallel  to  the  Y  sides.  The 
choice  of  quadrant  and  rotation  angle  of  9  is  arbitrary.  The  only  restrictions 
are  that  the  left  boundary,  XO  in  HULL,  is  the  minimum  X,  and  the  aft 
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boundary,  YO  in  HULL,  is  the  minimum  Y.  For  the  KEEL  run,  the  initiation  run 
for  HULL,  the  user  needs  to  supply  XO,  YO,  ZO,  and  the  grid  spacing  for  the 
computational  field.  For  HULLUP,  a  more  complete  knowledge  of  the 
computational  grid  is  needed. 

The  input  for  HULLUP  is  through  four  NAMELIST  lines  and  some  other  data  on 
the  INPUT  file,  and  from  the  donor  restart  file.  A  user  must  have  access  to 
file  NUHULUP  and  should  get  specific  instructions  for  program  HULLUP  from  it. 
The  donor  restart  file  must  be  ATTACHed  as  TAPE  9  and  the  output  files  must  be 
CATALOGed  in  the  runstream.  HULLUP  and  an  example  runstream  are  tabulated  in 
Appendix  A. 


V.  TEST  RUNS 

A  limited  number  of  test  problem  runs  have  been  made,  most  of  them  with 
the  original  coding  that  supplied  input  data  for  the  entire  plane.  The  first 
three  problems  were  simply  designed  to  get  the  program  running.  First  a  2-D 
cylindrical  donor  run  was  made.  It  had  a  sphere  of  high  pressure  gas  (E  in 
Figure  l)  with  radius  of  1000  cm  expanding  into  an  (IMAX,JMAX)  *  (32,16)  mesh 
of  toroidal  cells,  80x80  cm  in  cross  section.  In  the  2-D  run  (for  the 
cylindrical  space  A  in  Figure  1  represented  by  plane  B),  the  elevation,  Y, 
went  from  0  to  1200  cm  with  one  additional  layer  of  external  cells  and  the 
radius,  X,  from  0  to  2U80  cm  with  an  external  column  of  cells.  The  sphere  of 
high  pressure  gas  was  centered  at  (0,1200).  The  left,  bottom,  and  top 
boundaries  were  declared  reflective  and  the  right  boundary  transmissive. 
Setting  the  top  boundary  reflective  (i.e.,  symmetric)  is  all  right  for  the 
assumed  constant  atmosphere  with  no  gravity  until  the  reflected  shock  from  the 
ground  reaches  the  top  boundary. 

O 

The  sphere  of  high  pressure  gas  was  given  a  density  of  0.0381201*  g/cm  . 

3 

The  ambient  density  was  set  at  0.00120U12  g/cm  .  The  specific  internal  energy 

was  2.10371*  x  10^  ergs/g  everywhere  and  a  gamma  law  gas  was  assumed  with  gamma 
=  l.U.  Such  a  driver  gas  would  produce  a  3^.7  kPa  (50  psi)  shock  in  a 
straight  shock  tube.  With  spherical  expansion,  the  shock  front  should  be 
about  103  kPa  (15  psi)  at  1800  cm  radius.  The  2-D  HULL  run  was  initiated  at 
time  T  =  0.02  seconds  and  was  run  to  0.05  seconds. 

Two  3-D  problems  with  cubic  cells  80  cm  on  a  side  were  based  on  this  2-D 
donor  computation.  The  first  of  these  had  an  internal  region  (C  in  Figure  l) 
with  both  X  and  Y  from  0  to  1760  cm  and  Z  from  0  to  1200  cm.  This  was  run 
with  the  left,  aft,  top,  and  bottom  boundaries  reflective  and  the  right  and 
fore  boundaries  transmissive.  The  initiation  for  this  problem  was  from  the 
2-D  donor  through  KEEL  at  time  T  =  0.02  sec.  This  was  run  to  check  the 
starting  of  a  3-D  run  from  a  2-D  run  and  to  give  some  idea  of  the  difference 
in  3-D  and  2-D  results  for  this  problem.  (Starting  the  3-D  run  from  a  2-D  run 
will  tend  to  further  distort  the  boundary  of  the  high  density  sphere.)  As 
expected,  the  results  were  bilaterally  symmetric  about  X  =  Y.  Radial  symmetry 
was  reasonable  considering  the  coarseness  of  the  grids. 

The  other  3-D  run  based  on  the  80  cm  grid  donor  also  had  80  cm  cubic 
cells,  but  it  was  for  a  smaller  subspace  (D  in  Figure  l)  with  B0UND9  input 
from  the  donor  on  the  left,  right,  fore,  and  top  boundaries.  The  bottom  and 
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aft  boundaries  were  reflective.  A  9  by  7  by  8  mesh  of  cells  was  selected  with 
the  inner  region  from  X  =  880  to  1520  cm,  Y  =  0  to  1*80  cm,  and  Z  from  0  to  560 
cm.  This  was  also  initiated  at  time  T  *  0.02  sec  from  the  donor  restart  file. 
This  region  was  entirely  outside  the  high  pressure  sphere  defined  for  the  2-D 
run,  so  the  initiation  was  actually  for  an  ambient  region. 

B0UND9  input  planes  were  supplied  through  HULLUP  at  X  =  880,  the  left 
boundary,  X  *  1800,  the  center  of  the  external  plane  of  boundary  cells  on  the 
right,  Z  =  600,  the  center  of  the  top  plane  of  external  cells,  and  Y  *  500, 
the  center  of  the  external  plane  of  cells  on  the  fore  boundary.  The  grid  on 
these  planes  was  such  that  values  were  available  to  HULL  without  further 
smearing  due  to  interpolation  in  the  boundary  planes  (i.e.,  there  were  data 
points  at  the  cell  centers). 

Figures  2-7  show  the  HULL  records  of  overpressure  at  6  stations  (See 
Figure  l)  for  the  three  runs.  Stations  1,  2,  and  3  were  assigned  at  X  *  990, 
1250,  and  1500  cm,  respectively,  Z  =  30  cm,  and  Y  *  30  cm  for  the  3-D  runs 
(0=  0  for  2-D).  Stations  U,  5,  and  6  were  assigned  to  X  =  990,  1250,  and  1500 
cm,  respectively,  Z  -  500  cm,  and  Y  =  30  cm.  Although  assignment  is  to  a 
point  in  space,  the  recorded  values  are  for  a  cell  containing  the  point.  The 
results  are  best  thought  of  as  the  average  values  in  the  cell,  or 
alternatively,  as  the  values  at  the  cell  center.  For  example,  station  1  in 
the  3-D  runs  is  in  the  cell  between  X  =  960  and  10l*0,  Y  =  0  and  80,  and  Z  =  0 
and  80,  with  center  at  (X,Y,Z)  =  (1000,1*0,1*0). 

In  Figures  2-7,  and  the  others  in  this  report,  time  0  corresponds  to  the 
HULL  initiation  time  of  0.02  seconds,  and  of  course,  30  ms  corresponds  to  0.05 
seconds  in  HULL  time.  The  3  curves  are  labeled  DONOR,  3DTRAN ,  and  3DBND9  to 
identify  results  from  the  2-D  donor,  the  3-D  simulation  of  the  2-D  run,  and 
the  3-D  run  in  the  smaller  subspace  using  B0UND9  input,  respectively. 

There  is  a  shift  in  time  between  the  3-D  and  2-D  curves,  most  noticeable 
in  the  shock  arrival  time  and  rise.  Most  of  this  is  from  the  different 
treatment  of  overpressure  recording  for  2-D  and  3-D.  In  2-D  runs,  the 
pressure  at  the  beginning  of  a  time  step  is  recorded  with  the  station  data. 

In  3-D,  the  pressure  is  computed  from  energy  and  density  at  the  end  of  the 
time  step.  For  this  coarse  mesh,  the  time  steps  varied  from  about  0.05  to  0.9 
ms  after  5  ms.  The  time  shift  appears  smaller  for  Figure  1*  since  the  time 
step  is  less  than  0.3  ms  during  the  initial  pressure  increase.  (Incidentally, 
the  time  recorded  with  a  time  step  for  both  2-D  and  3-D  is  the  time  at  the 
beginning  of  the  time  step.  The  time  step  is  not  recorded,  and  not  all  time 
steps  have  output  so  the  correct  time  is  not  known.)  Considering  the 
coarseness  of  the  grid,  the  agreement  among  these  results  is  reasonable. 

Three  similar  runs  were  made  with  cells  approximately  1*0  cm  or.  a  side. 

(Halved  cell  sizes  were  desired  with  the  same  boundaries  as  before).  Since 
HULL  demands  an  odd  number  of  cells  inside  the  boundary  in  the  vertical 
direction,  there  was  a  minor  problem.  In  the  2-D  run  and  the  3-D  full 
simulation,  the  top  boundary  was  located  at  12l*0  cm  and  the  layer  of  cells 
between  1200  and  12^*0  cm  was  declared  to  be  ISLAND  cells.  This  created  a 
reflective  (symmetry)  boundary  at  1200  cm.  For  the  B0UND9  top  boundary  in  3-D 
an  odd  number  of  cells  between  Z  =  0  and  560  cm  was  required.  The  7  internal 
cells  with  DZ  =  80  cm  were  replaced  with  15  cells  with  DZ  =  38,  39,  or  1*0  cm.) 
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Figures  8-13  are  overpressure  records  for  these  "1*0  cm  grid"  runs  at  the 
same  stations  as  for  the  the  "80  cm  grid"  runs.  As  was  expected,  the  time 
shifts  are  less,  the  peaks  are  sharper  and  higher,  and  the  curves  have  more 
detail.  The  peak  pressures  for  the  full  3-D  run  is  always  somewhat  higher 

1* 

than  the  corresponding  2-D  peaks.  This  is  a  consistent  feature  of  HULL 
probably  due  to  minor  differences  in  the  finite  difference  algorithms. 

The  improvement  was  so  encouraging  that  a  2-D  donor  run  and  the 
corresponding  3-D  run  with  B0UND9  input  boundaries  were  run  with  cells 
approximately  20  cm  on  a  side  (overpressure  vs  time  plots.  Figures  1 L  —  1 9 )  and 
for  cells  approximately  10  cm  on  a  side  (Figures  20-25).  Except  for  station 
1,  the  2-D  donor  runs  and  the  B0UND9  3-D  runs  are  more  alike  for  the  20  cm 
grids  than  for  the  1»0  or  80  cm  grids.  The  time  shift  due  to  the  time  step  has 
been  reduced  to  about  0.2  ms  so  the  "smearing"  of  the  shock  front  by 
interpolation  for  and  from  B0UND9  input  boundary  data  is  more  evident  at  those 
stations  near  the  B0UND9  boundaries,  stations  1,  1*,  5,  and  6. 

One  would  expect  further  improvement  with  the  10  cm  grid  runs.  This  is 
not  evident.  There  are  several  reasons,  the  principal  one  being  additional 
smearing  due  to  interpolation  in  the  supplied  B0UND9  data  in  HULL.  For  al1 
the  coarser  grids,  the  HULLUP  coding  used  an  evenly  spaced  mesh  on  every  side 
between  the  B0UND9  values  at  their  edges.  By  choosing  the  appropriate  number 
of  mesh  divisions  there  were  points  in  the  center  of  the  cells  on  the  evenly 
spaced  sides  and  near  the  centers  in  the  Z  direction.  Such  a  fine  grid  was 
not  possible  for  the  10  cm  grid.  A  much  coarser  mesh  was  forced  by  computer 
space  limitations.  Hence,  there  was  smearing  in  both  the  HULLUP  interpolation 
and  the  HULL  interpolation.  This  difficulty  led  to  a  realization  that  a 
revision  was  necessary. 

Another  probable  source  of  differences  at  stations  1,  2,  and  3  is  that 
they  are  closer  to  the  reflective  bottom  boundary  in  2-D  than  in  3-D.  For  2-D 
they  are  in  the  third  layer  up  between  20  and  30  cm,  and  for  the  3-D  they  are 
in  the  fourth  layer  between  28.5  and  38  cm. 

The  notches  in  the  initial  pressure  rise  in  Figures  20  and  21  are  probably 
due  to  an  initial  decay  of  the  incident  shock  followed  by  the  reflected  shock  . 
(The  stations  are  3  or  b  cells  from  the  reflective  bottom. )  An  alternative 
explanation  is  to  assume  a  Mach  reflection  which  needs  time  to  develop  in  the 
3-D  restart  run. 

Station  1  is  near  the  left  B0UND9  input  boundary  in  the  3-D  runs. 
Apparently  it  is  too  close  for  a  sharp  reflected  shock  to  develop. 

Figures  26-31  show  a  comparison  of  the  overpressure  from  the  four  donor 
runs  at  each  of  the  6  stations.  The  plots  for  stations  U,  5,  and  6  (Figures 
29,  30,  and  31 )  are  much  as  one  would  expect:  the  peak  pressure  increases 
with  decreasing  cell  size,  the  rise  and  fall  from  peak  pressure  is  much 
sharper,  and  the  curves  are  more  detailed  as  cell  size  decreases.  The  same 
results  are  not  so  true  at  stations  1,  2,  and  3  (Figures  26,  27,  and  28) 
because  of  the  offset  in  position  and  the  rapid  decay  of  the  peak  reflected 
pressure  near  the  reflective  bottom  boundary.  The  results  seem  to  converge 
ar.d  the  20  cm  and  10  cm  grid  results  agree  fairly  well.  Better  agreement 
would  be  welcomed  at  stations  3  and  6  near  the  end  of  these  short  runs. 
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Figures  s2-j7  show  a  comparison  of  the  overpressures  for  the  four  BQUND9 
driven  o-D  runs  at  the  same  b  stations.  Overpressure  for  the  donor  run  with 
the  10  cm  grid  is  included  as  an  assumed  "correct"  curve.  Here  again  the 
convergence  looks  good  witn  the  20  cm  grid  results  and  10  cm  grid  results 
being  reasonably  alike. 

As  was  stated  earlier,  the  presented  results  were  all  from  an  early 
version  of  HULLUP  and  UPDATE  changes  to  HULL  for  BQUND9  boundaries.  The 
HULLUP  program  was  revised  to  supply  line  by  line  data  on  the  input  planes  at 
grid  points  that  may  be  specified.  The  B0UND9  coding  was  revised  to  use  the 
line  by  line  input.  Spatial  interpolation  was  retained,  but  cell  centered 
input  is  desirable.  The  results  for  the  BQUND9  run  with  the  80  cm  grid  were 
duplicated  exactly  with  the  new  coding. 

The  BUUND9  run  witn  the  10  cm  grid  was  also  rerun  with  the  new  coding. 
Figures  38-41  show  the  overpressure  records  at  the  stations  1  through  4, 
respectively,  for  the  2-D  donor  run  and  tne  3-D  BGUND9  driven  runs  with  the 
early  coding,  and  with  the  revised  coding.  The  revised  coding  did  not  make  as 
much  difference  as  the  author  expected.  There  is  some  separation  of  the  two 
3-D  generated  curves  at  station  1 ,  and  the  peak  values  at  stations  2  and  3  are 
noticeably  different.  The  results  at  stations  4,  5,  and  6  are  nearly 
identical . 

There  is  no  reason  to  doubt  that  the  HULLUP  coding  will  properly  prepare 
B0UND9  input  for  a  o-B  hULL  run  trom  a  cylindrical  2-D  run,  or  to  doubt  that 
the  B0UND9  changes  to  HULL  will  give  a  satisfactory  d-D  continuation, 
however,  there  are  a  number  of  things  tnat  were  not  tested:  there  was  only 
one  computation  field  location,  no  test  with  a  target,  no  test  with  the  bottom 
or  aft  ooundary  a  BUUND9  type,  no  d-D  start  with  a  region  partially  filled 
with  non-ambient  gas  to  check  for  conflict  between  tne  K.EEL  set  up  and  tne 
BuUND9  boundary,  no  testing  of  the  effect  of  changing  grid  size,  no  coordinate 
shifting,  and  only  accidental  testing  of  one  or  two  abort  situations. 


VI.  CONCLUSIONS 

The  HULL  B0UND9  coding  for  j-D  HULL  and  the  HULLUP  program  to  prepare  the 
BOUND9  input  from  a  cylindrically  symmetric  2-D  HULL  run  seem  to  be  working 
properly.  Any  use  of  them  should  be  carefully  monitored;  more  testing  is  in 
order  before  making  any  permanent  changes  to  our  HULL  code. 

The  o-D  run  to  be  made  must  be  thoroughly  planned.  The  locations  of  the 
BOUNDS  input  planes  must  be  exact,  on  the  low  index  boundaries  and  in  the 
centers  of  the  external  cells  on  tne  high  index  boundaries.  A  grid  of  mesh 
points  on  the  BOUND9  input  planes  matching  the  centers  of  the  3-D  HULL  cells 
is  desirable  since  it  reduces  smearing  due  to  interpolation.  A  careful 
analysis  of  expected  velocities  and  desired  run  time  is  needed  to  determine 
the  placement  of  boundaries  for  the  d-D  (and  the  2-D)  run. 

The  HULL  HCuND9  coding  will  accept  input  on  the  boundaries  from  any  source 
as  iong  as  it  is  in  the  proper  format  in  tne  input  files.  For  example,  input 
on  a  side  that  i3  constant  in  space  and  varies  with  time  might  De  useful  and 
would  be  fairly  easy  to  construct.  by  coding  a  program  litce  nbLLUP,  input 
trom  any  nydrocooe  could  be  usea. 
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B0UND9  type  input  for  2-D  codes  has  not  been  ooded.  A  2-D  Cartesian,  or 
new  cylindrical,  run  may  be  initiated  in  KEEL  from  a  2-D  cylindrical  donor.  A 
large  enougn  region  would  be  needed  to  prevent  false  signals  from  necessarily 
Incorrect  boundaries  and  such  a  region  may  not  be  possible.  A  3-D  run  could 
be  made  in  tne  same  way  but  the  region  needed  would  have  to  be  much  larger 
than  one  with  B0UND9  boundaries  for  the  same  uncon tamina ted  run  time. 
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Figure  1.  Computation  Regions. 
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Figure  it.  Comparison  of  Overpressure  From  a  2-D  Donor  Hun 
a  Comparative  3-D  Run,  and  a  3-D  Imposed-Boundary  Run 
for  an  80  cm  Mesh  at  Station  3. 
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Figure  6.  Comparison  of  Overpressure  From  a  2-D  Donor  Run, 
a  Comparative  3-D  Run,  and  a  3-D  Imposed-Boundary  Run 
for  an  80  cm  Mesh  at  Station  5. 
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Figure  8.  Comparison  of  Overpressure  From  a  2-D  Donor  Run, 
a  Comparative  3-D  Run,  and  a  3-D  Imposed-Boundary  Run 
for  a  1*0  cm  Mesh  at  Station  1. 
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Figure  10.  Comparison  of  Overpressure  From  a  2-D  Donor  Run 
a  Comparative  3-D  Run,  and  a  3-D  Imposed-Boundary  Run 
for  a  U0  cm  Mesh  at  Station  3. 
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Figure  11.  Comparison  of  Overpressure  From  i 
a  Comparative  3-D  Run,  and  a  3-D  Imposed 
for  a  UO  cm  Mesh  at  Station  h 
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Figure  13.  Comparison  of  Overpressure  From  i 
a  Comparative  3-D  Run,  and  a  3-D  Imposed 
for  a  1*0  cm  Mesh  at  Station  6 
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Figure  l8.  Comparison  of  Overpressure  From  a  2-D  Donor  Run 
3-D  Imposed-Boundary  Run  for  a  20  cm  Mesh  at  Station  5. 
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Figure  31.  Overpressure  From  2-D  Donor  Runs  at  Station  6. 
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Figure  39.  Comparison  of  Overpressure  From  the  Original  and 
the  Revised  BOUND9  Coding  at  Station  2. 
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Figure  UO.  Comparison  of  Overpressure  From  the  Original  and 
the  Revised  B0UND9  Coding  at  Station  3. 
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Figure  1*1.  Comparison  of  Overpressure  From  the  Original  and 
the  Revised  B0UND9  Coding  at  Station  1*. 


I 


61 


REFERENCES 


M.  A.  Fry,  R.  E.  Durrett,  G.  P.  Ganong,  D.  A.  Matuska,  M.  D.  Stucker,  B. 

S.  Chambers,  c.  E.  Needham,  and  C.  D.  Westmoreland,  "The  HULL 
Hydrodynamics  Computer  Code,"  AFWL-TR-76-1&3,  US  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  NM,  September  1976.  (AD  #B014070L) 

J.  A.  hasdal,  B.  S.  Chambers,  and  R.  W.  Clemens,  "Support  to  BRL:  HULL 
Code  implementation  on  a  CDC  7600,"  SAI-eG-701-AQ,  Science  Applications 
inc.,  McLean,  VA,  August  1979. 

B.  S.  Chambers,  and  J.  D.  Wortman,  "Two-Dimensional  Shore  (Partial  Island) 
Ceils  for  BRL  HULL,"  ARBRL-CR-00497 ,  US  Army  Ballistic  Hesearch 
Laboratory,  Aberdeen  Proving  Ground,  MD,  December  1962.  (AD  #A12j357) 

R.  E.  Lottero,  J.  D.  Wortman,  B.  P.  Bertrand,  and  C.  W.  Kitchens,  Jr., 
"Three-Dimensional  Oblique  Shock  Diffraction  Over  a  Rectangular 
Parallelepiped:  Computational/Experimental  Comparison,"  ARBRL-TR-02443, 

US  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD, 
November  1982.  (AD#A 122254) 


APPENDIX  A 


TABULATION  OF  PROGRAM  HULLUP 

Appendix  A  is  a  listing  of  a  SCOPE  2  runstream  to  use  HULLUP  to  prepare 
BOUND9  boundary  input  with  a  listing  of  file  NUHULUP  which  contains  the  HULLUP 
program. 

Problem  8405.09  is  the  "80  cm  grid"  donor  file  described  in  section  V. 
Problem  0406.20  (or  b4lb.2G)  is  the  corresponding  SQUND9  controlled  3-D  test 
run.  Restart  output  from  problem  8405.09  is  found  in  the  MFZ  file 
HULL8405P09. 

Since  data  was  wanted  at  all  the  available  times  on  tnis  restart  file,  the 
only  input  for  namelist  INDAI1  was  the  donor  problem  number. 

Namelist  NDAT1  and  the  10  lines  following  it  describe  the  locations  of  the 
boundaries  for  the  j-D  run  and  the  number  of  and  locations  of  the  data  points 
on  these  sides.  Since  KNM  was  negative,  the  mesh  points  were  specified  (the 
10  lines  after  namelist  NDAT1).  Equivalent  output  data,  for  this  case,  could 
have  been  formed  by  setting  1NM  =  18,  YNM  =  14,  and  KNM  =  16  in  namelist 
NDA‘1'1 ,  with  no  additional  input.  This  would  give  more  mesh  points,  but  the 
only  effective  ones  would  be  those  at  the  cell  centers.  This  is  possible  only 
because  the  j-D  run  has  equal  cell  sizes. 

Namelist  NDAT2  informs  HULLUP  tnat  B0UWD9  boundary  data  is  to  be  saved  for 
X  =  XN1  and  XN2  (the  left  and  right  boundaries),  Y  =  YN2  (the  fore  boundary), 
and  L  =  ZN2  (the  top  boundary).  Namelist  TABDAT  tells  HULLUP  the  number  of 
times  that  output  is  wanted  for  each  of  these  boundaries. 


**  RUN  STREAM  FOR  HULLUP  ** 


WORTMAN(5TMFZ,T350,  P5) 

ACCOUNT !*****♦)  WORTNAN  B309  X602B 


COMMENT.  PROGRAM  TO  PRODUCE  BOUNDARY  INPUT  FOR  HULL  FROM  AN  OLD  HULL 
COMMENT.  SET  FOR  PROB  8*16.20  FROM  PROB  6*05.09.  8/10/B* 

COMMENT.  THIS  SHOULD  REPRODUCE  8*06.20  UITH  THE  NEW  HULLUP 

COMMENT.  THAT  PRODUCES  OUTPUT  FILES  WITH  ROW  BY  ROW  RECORDS. 

COMMENT. 

COMMENT.  CHECKLIST  FOR  A  NEW  PROGRAM 

COMMENT.  CHANGE  IDENTIFICATION  LINES*  THE  INPUT  FILE  FOR  TAPE9, 

COMMENT.  AND  SET  CATALOGGING  OF  OUTPUT  FILES. 

COMMENT,  ADJUST  NAMELIST  ANO  OTHER  INPUT  BELOW. 

COMMENT.  CHECK  PARAMETER  STATEMENT  VALUES,  FRONT  OF  FILE  NUHULUP. 
COMMENT. 

ATTACH (TAPE 9, HULL8*05P09,  ID* JDW) 

COMMENT.  REQUEST  PERMANENT  FILES  FOR  POSSIBLE  B0UND9  OUTPUT. 

REQUEST! TAPE11, *PF ) 

REQUEST(TAPE12,*PF> 

REQUEST(TAPE13,*PF) 

RE0UEST(TAPE1*,*PF  > 

REQUEST(TAPE15,*PF) 

REQUEST (TAPE16»*PF) 

COMMENT.  TEMPORARY  SET  UP  FOR  COMPILING  HULLUP. 

BEGIN (GETMFA,FILE,LF*HULLUP,PF*NUHULUP, UN* JDW) 

UPDATE! N, I*HULLUPI 
COMMENT.  FTN5( I, L*0, OP  T*  2 ) 

FTN5(I,L0*S/A/R/M, QPT*2) 

LG0!*PL-10000J 

COMMENT.  CATALOG  FILES  OF  BOUNDARY  DATA  HERE. 

CATALOG (TAPE11,P6*1620L8, ID- JDW) 

CATALOG! TAP E12, P8A1620RB* ID- JDW) 

COMMENT.  CATALOG! TAPE  13, P8*1620BB, ID- JDW ) 

CATAL0G!TA»E1*,P8*1620TB,  ID- JDW) 

COMMENT.  CATALOG! T  APE  15, P8A1620AB, ID- JDW) 

CATALOG! TAP E16, P 6* 162 OF B, IO-JOW) 

*EOR 

SIND ATI  PR0BIN-8*05, 09  S 

S NO ATI  XN1-880., XNM-1560. , INM-9, YN1 -0.0, YNM-520.0, JNM- 7, 

ZN1-0.0, ZNH-600. , KNM--BS 
920.0 

80.0  8 

-1 

*0.0 

80.0  * 

80.0  2 

-1 

*0.0 

80.0  7 

-1 

SNDAT2  NBND  -  1,1,0,1,0,18 
STABDAT  LASTPR  -  5,*,0,*,0,3$ 


**  HULLUP  ♦* 


FILE  NUHULUP.  HULLUP  POP.  30  B0UND9  INPUT  FROM  20  DONOR. 


FIRST.  THE  3  COMMON  DECKS  FOR  SUBPROGRAMS  IN  HULLUP. 
THE  FIRST  2  HAVE  DONOR  AND  RECIPIENT  ARRAY  DIMENSIONS. 
PARAMETERS  SET  P  OR  PRDR  8A07.24  FROM  8A07.22. 

THEY  MAY  BE  CHANGED  IF  NECESSARY. 


*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

•COMDECK  C0IN1 

C  VARIABLE  DIMENSIONS  FOR  DONOR  ARRAYS. 

PARAMETER! I  IB IG- 2 57, I J B 16- 123, I  KB IG- 2, I L B IG- 5000) 

COMMON  /CD IN  1  /  XIUIBIG)»DEIXI(IIBIG),YI!IJBIG),0£LYI(IJBI6) 

♦  f HYOROI  (ILB  IG)  , II  (  IKBIGI.OELZK IKBIG) 

•COMDECK  C0NU1 

C  VARIABLE  DIMENSIONS  FOR  RECIPIENT  ARRAYS 

PARAMETER!  NI  JKM  AX-200,  MB  IG-1000) 

COMMON  /CDNU1  /  XNINIJKMAX  ),YN(NIJKMAX), ZN(NIJKMAX),BNO(NLBIG), 

♦  BND2(NLBIG ) 

•COMDECK  CDCOM 


C  OTHER  COMMON  TERMS  FOR  HULLUP. 

C0HM0N/CDIN3 /  IN9, PRQBIN.C START, IIMA X, IIMAX1, UMAX, IJMAX1, 

♦  IKMAX,IKMAX1,IGE0M, INH,  IROWPB, INHPB# 

♦  ITYPE.ITIME, CYCLE I.TIMEIN 

C0MM0N/CDNU3/  XN1,XNM,YN1,YNM, ZN1,  Z NM, INM.JNM, KNM, TSH, X5H»YSH#ZSH, 

♦  NNH.NBND!  t> ) ,  NF  OUT  (  6  )  ,  X  YZ  B(  6),  ICON  V,  NXP  L»  N  YPL»NZPL 
•DECK  HULLUP 

PROGRAM  HULLUP! INPUT , OUT  PUT, T A PE9» T APE  11, T APE12»T APE13»T APE  14# 

♦  TAPE15,T  APE  16, TaPE5« INPUT, TA PE6-0UTPUT, T APE1, TAPE2 ) 

***** 

HULLUP  EXTRACTS  BOUNDARY  INPUT  VALUES  FOR  A  NEW  HULL  RUN  FROM 
A  HULL  RESTART  FILE.  THE  OLD  HULL  RUN  IS  CALLED  DONOR  OR  OLD. 
PARAMETERS  FOR  OP  FROM  IT  MAY  HAVE  PREFIX  OR  SUFFIX  I  OR  IN  . 
THE  NEW  HULL  RUN  WILL  BE  CALLED  RECIPIENT  OR  NEW.  PARAMETERS 
FOR  IT  MAY  HAVE  PREFIX  OR  SUFFIX  N  OR  NU . 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c** 

c 


AN  ATTEMPT  HAS  BEEN  MADE  TO  ALLOW  FOR  FUTURE  CHANGES.  THIS 
FIRST  CODING  WILL  BE  FOP  A  CYLINDRICAL  SYMMETRIC  DONOR  AND 
PREPARE  BOUNDARY  DATA  FOR  A  30  CARTESIAN  RECIPIENT. 

WE  HAVE  MADE  OTHER  ASSUMPTIONS. 

THE  RESTART  DATA  IS  ON  ONE  FILE,  TAPE9.  (WE  WILL  SET  IN9-9 
AND  USE  IN9  TO  MAKE  POSSIBLE  CHANGE  SIMPLER.) 

THE  COORDINATE  POINTS  OF  DATA  ON  IN9  HAY  CHANGE  BETWEEN  DUMPS. 
WE  ASSUME  NO  ISLAND  OR  SHORE  CELLS  ON  NEW  BOUNDARY  PLANES. 
NNH-5  FOR  30.  !IE  5  HYDRO  VARIABLES  OUTPUT) 

FOR  CYLINDRICAL  DONOR ,  RADIAL  VELOCITY  AT  R  ■  0.0  SIT  TO  0.0. 
***** 

WE  HAVE  VARIABLE  PARAMETERS  AND  ARRAYS  IN  COMDECKS. 
VARIABLE  DONOR  PARAMETERS  IN  CDINI. 

VARIABLE  RECIPIENT  PARAMETERS  IN  CDNU1. 

***** 

INPUT 

***** 

INPUT  IS  THRU  NAME  LIST  INPUT,  OR  DATA,  ON  TAPE 5  AND  FROM  IN9. 


NAMELIST  /IN0AT1/ 


(THIS  IS  READ  IN  HULLUP,  I.E.  MAIN) 


IN9 


FILE  FOR  OLD  HULL  RUN  RESTART  DUMPS.  DEFAULT  •  9. 

THIRD  FILL  ON  PROGRAM  CARD.  MUST  ATTACH  IN  RUN  STREAM. 


PROBIN 


OLD  PROBLEM  NUMB fc .  ON  INPUT  FILE.  REQUIRED. 


C5TART 


STARTING  CYCLE  NUMBER  OLD  DATA.  DEFAULT 


0.0. 


"V  .V  .  Ti  \ 


L. 

s 


t 


c« 


¥ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c** 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CSTOP  - 


UPPER  SOUND  FOR  INPUT  CYCLES.  WILL  PROCESS  INPUT 
CYCLES  UNTIL  CYCLEI  GT  CSTQP  OR  INPUT  FILE  ENOS. 
DEFAULT  >  9999999. 


IN9  HAS 
DUMP  IS 


VARIABLE  NUMBER  OF  HULL  RESTART  DUMPS.  EACH 
HEADER.  A  Z-BLOCK,  A  VERTEX  RECORD.  AND  HYORO  DATA. 

THE  HEADER  RECORD  IS  A  NUMBERS. 

-  SIGNAL.  995.0  FOP  A  RESTART  DUMP.  666.0  FOR  ENO  OF  FILE. 

-  DONOR  PROGRAM  NUMBER.  THIS  MUST  BE  THE  SAME  AS  PROBIN. 

-  CYCLE  NUMBER. 

-  TIME. 

THE  Z-BLOCK  IS  A  200  WORD  RECORD.  92  NAMES  AND  VALUES.  AND  ID. 
WE  WILL  EXTRACT  7  VALUES. 


UMAX 

UMAX 

IKMAX 

INH 

ID1HEN 

IGEOM 

IROWPB 


X  DIRECTION. 
Y  DIRECTION. 
2  DIRECTION. 


THE  NUMBER  OF  CELLS  IN  THE 
THE  NUMBER  OF  CELLS  IN  THF 
THE  NUMBER  OF  CELLS  IN  THE 

NUMBER  OF  HYORO  VAR  I  AS L£ S VC ELL .  WE  USE  ONLY  4. 

INPUT  PROBLEM  DIMENSION.  (WE  ASSUME  ONLY  2D  NOW) 

INPUT  GEOMETRY.  1-CARTESIAN,  2-C YL INORIC AL  ASSUMED. 

ROWS  STORED/BLOCK.  (EACH  BLOCK  IS  A  RECORD). 

VERTEX  RECORD 

2D  -  (XI  ( I ),  1*2,  IIMAX+1),  (YKJ  )  ,  J  *1.  IJMAX+l)  •  Xl(l).0.0  ASSUMED. 

30  -  (XI  (I),  I-l,  IIMAX+1),  (YI(J>,J»1,IJMAX*1),(ZI(K),K-1,IKNAX*1). 

HYDRO  VALUES.  INH* I IMAX*IROWPB  VALUES  PER  RECORD. 

FOR  EACH  CELL  IN  2D,  THE  FIRST  5  HYDRO  VALUES  ARE. 

-  PRESSURE  THAT  IS  NOT  USEFUL.  IT  IS  FOR  ANOTHER  TIME. 

-  X  VELOCITY  COMPONENT  (RADIAL  VELOCITY  FOR  CYLINDRICAL).  (CM/SEC) 

-  Y  VELOCITY  COMPONENT  (AXIAL  VELOCITY  FOR  CYLINDRICAL).  (CM/SEC) 

-  SPECIFIC  ENERGY.  (ERGS/GM) 

-  MASS  IN  THE  CELL.  (GM) 


NAMELIST  /NDAT1/ 


(READ  IN  STRTNU) 


THIS  DESCRIBES  THE  NEW  B0UND9  PLANES.  SOME  EXPLANATION  IS  NEEDED. 
FOR  LOW  INDEX  B0UND9  BOUNDARIES,  HULL  NEEDS  VALUES  CENTERED  ON 
THE  OUTSIDE  EDGE  OF  EACH  BOUNDARY  CELL.  FOR  HIGH  INDEX  B0UND9 
BOUNDARIES,  VALUES  ARE  NEEDED  IN  THE  CENTER  OF  THE  EXTERNAL  CELLS. 
HULL  WILL  INTERPOLATE  IN  THE  SPACE  DEFINED  ON  A  PLANE.  IT  WILL 
NOT  EXTRAPOLATE.  ACCURACY  WILL  BE  BEST  IF  THE  CELL  CENTERS  FOR 
THE  NEW  HULL  GRIO  ARE  SPECIFICALLY  INCLUDED  IN  THE  B0UND9  MESH. 
XN1, XNM.YNl.YNM, ZN1, ZNM  ARE  THE  LOCATIONS  OF  BOUNDARY  PLANES. 
(XN(I),I«1»INN>,(YN<4),J«1,4NM),(ZN(K),K»1,KNK>  ARE  THE  MESHES 
ON  THE  BOUNDARY  PLANES. 

2D  DONOR  CELL  VERTICES  ARE  ( XI ( I >, 1-1, IMAX*1>, ( YI { 4 > ,4»1,4NAX*1 ) 
XO  ■  XI  (1 )  IS  0.0.  XKIMAX+1)  AND  YI(IMAX*1)  EXTERNAL. 

THE  HYDRO  IN  EXTERNAL  CELLS  IS  NOT  USEABLE. 

WE  HAVE  THE  FOLLOWING  RESTRICTIONS  FOR  A  3D  RECIPIENT. 

FOR  XMAX-XK  IMAX),  YMAX-  YI(JMAX) ,  X0-XK1),  AND  YO-YI(l), 
(SQRT(XN(I)**2+YN(4)**2)  LT  XMAX  AND  YO  LE  ZN( K)  LE  YMAX. 

ALL  POINTS  ON  AN  OUTPUT  PLANE  MUST  BE  INCLUDED  INSIDE  THE  SPACE 


DEFINED  BY  THE  2D  DONOR  GRID. 

MAKING  2  OR  3  SEPERATE  RUNS  NAY  BE  MORE  CONVENIENT. 


XN1 

XNM 

YN1 

YNM 

ZN1 

ZNM 


-  -  NAMELIST  TERMS  -  - 
LEFT  BOUNDARY.  DEFAULT  -0.0. 

RIGHT  EXTERNAL  BOUNDARY  CELL  CENTER.  REQUIRED. 

2D  BOTTOM  BOUNDARY,  3D  AFT  BOUNDARY.  DEFAULT  ■ 

2D  TOP.  3D  FORE,  EXTERNAL  BOUNDARY  CELL  CENTER. 

3D  BOTTOM  BOUNDARY.  DEFAULT  •  0.0. 

3D  TOP  EXTERNAL  BOUNDARY  CELL  CENTER.  REQUIRED  FOR  3D. 


0.0. 

REQUIRED. 


INM 

JNM 

KNM 


NUMBER  OF  X  POINTS  IN  OUTPUT  ARRAYS. 
NUMBER  OF  Y  POINTS  IN  OUTPUT  ARRAYS. 
NUMBER  OF  Z  POINTS  IN  OUTPUT  ARRAYS. 


DEFAULT 

DEFAULT 

DEFAULT 


200. 

200. 

200. 


c 


c 

TSH  -  SHIFT  IN 

TINE. 

TINENU  •  TINEIN  -  TSH. 

DEFAULT  ■  0.0. 

c 

XSH  -  X  SHIFT. 

XNU  « 

XIN  -  XSH. 

DEFAULT  • 

0.0 

c 

YSH  -  Y  SHIFT. 

YNU  ■ 

YIN  -  2SH. 

DEFAULT  • 

0.0 

C 

ZSH  -  Z  SHIFT. 

ZNU  - 

XIN  -  XSH. 

DEFAULT  - 

0.0 

c  - 

c 

C  THE  OUTPUT  NESH  CAN  BE  DEFINED  IN  2  NUTUALLY  EXCLUSIVE  WAYS 

C  (1)  IF  I NM> O#  JNH>  0»  ANO  KNH>0»  THEN 

C  XN < I )  -  XN1  ♦  ( I— 1 ) * ( XNN-XN1 )  /  ( INM-1 )  FOR  I-l/INN 

C  YNt  J  )  •  YN1  ♦  ( J-1)*(YNN-YN1) /( JNN-1)  FOR  J-l.JNN 

C  ZN( K )  •  ZN1  ♦  (K-l I *( ZNN-ZN1 ) /(KNN-1)  FOR  K-1,KNN 

C 

C  ( 2 )  IF  ANY  OF  INN,  JNN,  OR  KNM  <  0,  THEY  ARE  SET  POSITIVE  AND 

C  THE  FOLLOWING  LINES  ARE  READ  WITH  FORNAT  <E15.8,I5): 

C  XN ( 1  ) 

C  0X1, N1  —  FRON  WHICH  XN(I)  -  XN(I-1)#DX1  FOR  1-2  TO  Nl  +  1 

C  0X2, N2  —  FRON  WHICH  XN(I)  ■  XN(I-1)*DX2  FOR  I-Nl+1  TO  N1+N2+1 

C  . 

C  DXL, NL  —  FOR  XN(I)  UP  TO  I -1-N1*N2+. . . NL  (•?  INN). 

C  ANY# -1  --  TERNIMATOR  FOR  X  GRID. 

C  NUST  HAVE  1*N1>N2*. ..♦NL  -  INN. 

C 

C  FOLLOWEO  BY  SINILAR  INPUT  FOR  YN(J>  AND  THEN  ZN(K>. 

C 

C  THEN, 

C 

C**  3  NANELIST  /N0AT2/  NORE  INPUT  FOR  NEW  PROBLEH.  (READ  IN  STRTNU) 
C 

C  NGEON  -  NEW  GEONETRY.  2  IS  CYLINDRICAL,  1  IS  CARTESIAN. 

C  DEFAULT  IS  1  (CARTESIAN) 

C 

C  NOINEN  -  OUTPUT  OINENSIONS.  2  OR  3.  DEFAULT  •  3. 

C 

C  NNH  -  NUNBER  OF  OUTPUT  HYDRO  VARIABLES  PER  POINT. 

C  DEFAULT  -  5. 

C 

C  (NBND(I), 1-1,6)  -  SIGNAL  FOR  BOUNDARY  1. 

C  FIND  OUTPUT  FOR  BOUNDARY  I  IF  NBND(I)  -  1.  DEFAULT  •  0. 

C  BOUNDARIES  IN  ORDER  ARE  LEFT»R IGHT# BOTTON, TOP, AFT, FORE . 

C  (AT  LEAST  ONE  NBNO(I)  •  1#  OR  A  NULL  REQUEST.) 

C 

C  (NFOUT(I), 1-1,6)  -  FILE  TO  STORE  NBND(I)  RESULTS  IN. 

C  DEFAULT  -  I  ♦  10. 

C  THESE  ARE  THE  ATH  THRU  9TH  FILE  ON  PROGRAN  LINE. 

C  THOSE  WITH  NBNO(I)  -  1  TO  BE  CATALOGED  IN  THE  RUNSTREAH. 

C 

C**  A  NANELIST/TABDAT /  CONTROLS  TAB  OF  B0UND9  FILES.  (READ  IN  REWRIT) 
C 

C  (INITPR(I), 1-1,6)  -  FIRST  DUMP  OF  NFDUT(I)  TO  TAB.  DEFAULT-1. 

C 

C  (LASTPR( I ) , I- 1, 6 )  -  LAST  TINE  DUMP  OF  NFOUT(I)  TO  TAB.  DEFAULT-0 

C 

C  (INITRO(I), 1-1,6)  -  FIRST  ROW  OF  NFOUT(I)  TO  TAB.  DEFAULT  •  1. 

C 

C  ( L ASTRO! I) # I- 1, 6 )  -  LAST  ROW  OF  NFOUT(I)  TO  TAB.  DEFAULT  -  10. 

C  - 

C 

C  *  *  *  SUBROUTINES  •  ♦  * 

C  MOST  OF  THE  PROGRAN  IS  IN  THE  NAIN  PROGRAN  CALLED  HULLUP. 

C  THERE  ARE  NOW  (8/08/8A)  SIX  SUBROUTINES. 

C 

C  STRTIN  -  INITIATION  FOR  READING  THE  INPUT  FILE  (IN9). 

C  READ  HEADER  RECORD  FRON  IN9.  CHECK  THAT  OLD  PROGRAN  NUNBER 
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IS  PROBIN.  READ  THE  FIRST  Z-BLOCK,  EXTRACT  AND  COMPUTE  TERMS 
CHECK  THAT  INPUT  ARRAY  MAXIMA,  (  1 1  B 1 6, I J BIG, I  KB  I G» I  LB  I G ) »  ARE 
LARGE  ENOUGH. 

STRTNU  -  INITIATION  FOR  THE  NEW  HULL. 

READ  NAMELIST  INPUTS  NDAT1  AND  NDAT2.  SET  UP  OUTPUT  POINTS 
AND  PARAMETERS.  CHECK  THAT  OUTPUT  ARRAY  MAXIMA,  (NlJKMAX 
AND  NLBIG),  ARE  BIG  ENOUGH. 

NXTIH(IEND)  -  LOCATES  NEXT  TIME  DUMP  ON  IN9.  (20  ONLY) 

IF  ITIME  «  0,  SEARCH  FOR  C YC LE I -CST ART .  READ  DUMP  HEADER. 
BYPASS  THE  Z-BLOCK.  READ  2D  VERTICES  AND  SET  UP  VOLUME  TERMS 
AND  MIDCELL  VALUES.  SET  I  END  TO  0,  1,  OR  -I.  0  FOR  HYDRO 
DATA  AVAILABLE,  1  FOR  NORMAL  FILE  END,  -1  FOR  READ  TROUBLE. 

FIND I(I»Z,Y«IM)  FIND  I  FOR  LINEAR  INTERPOLATION  IN  Y. 

I«1  IF  Z  <  Y<  2) »  I-IM-1  IF  Z  »  Y(IM-I),  OR  Y( I)  <  Z  <■  Y(I*1) 

NXTBLMJP,  IHB )  -  GET  HYDRO  INPUT  DATA.  (2D  INPUT  ONLY) 

IH8  IS  THE  LAST  HYDRO  BLOCK  READ,  JP  IS  ROW  OF  DATA  NEEDED. 
THE  PROGRAM  CHECKS  FOR  BLOCK  NEEDED.  IT  IS  PUT  IN  HYDROIt— ) 

WRITPL(IBND,NF,X,XS,Y,NY, YS,2,NZ,ZS,B,NL)  -  TEMPORARY  DUMP  TO  NF. 
ROW  DUMP  FOR  A  BOUNDARY  PLANE.  IF  FIRST  ENTRANCE  FOR  THIS 
TIME,  DUMP  HEADER  AND  GRID.  AFTER  NZ  ENTRANCES  FOR  FILE  NF 
THERE  ARE  2*NZ  RECOROS  DN  FILE  NF. 

1  -  HEADER  OF  7  TERMS. 

555. , PROBIN, TIMIN,X,NNH,NY,N2 

2  -  LOCATION  OF  POINTS  IN  THE  PLANE. 

(Y(J),J«1,NY>, (Z(K),K-1,NZ) 

2*NZ  -  HYDRO  DATA  FOR  THE  FOR  KTH  ROW  IN  PLANE,  (K-1,NZ). 

( BND ( L ) , L  *1, NL  )  WHERE  NL  *  NNH*NY . 

***  AFTER  THE  FINAL  TIME  DUMP  A  HEADER  IS  DUMPED,  ALL  666. 

REWRIT ( IBND,NF )  -  REPLACES  CONTENT  OF  FILE  NF  WITH  FINAL  OUTPUT 

THE  CONTENTS  OF  NF  ARE  PUT  INTO  2  FILES  THEN  COLLATED 
AND  PUT  BACK  ON  NF  WITH  2  SUCCESSIVE  TIMES  AT  EACH  DUMP. 

AGAIN  THERE  ARE  2  ♦  NZ  RECORDS  PER  TIME  DUMP. 

EXCEPT  FOR  2  SUCCESSIVE  TIMES,  TI1,TI2,  TERMS  AS  IN  WRITPL. 

1  -  HEADER  OF  8  TERMS. 

555.0, PROBIN,  T II, TI 2, X, NNH, NY, NZ 

2  -  LOCATION  OF  POINTS  IN  THE  PLANE 

(YN( J),J«1,NY),(ZN(K),K-X,NZ) 

2*NZ  -  HYDRO  DATA  AT  Til  AND  THEN  TI2  FOR  ROW  K,  (K»1,NZ>. 
(BND(L),L*1,NL),(BND2(L),L*1,NL)  WHERE  NL*NNH*NY . 

BNO(L)  IS  HYDRO  AT  TIME  Til,  BND2 ( L )  HYDRO  AT  TI2. 

***  AFTER  THE  FINAL  TIME  DUMP  A  HEADER  IS  DUMPED,  ALL  666. 


*  *  *  *  *  GLOSSARY  OF  TERMS  ***** 

SYMBOLS  -  *  NAME  IN  A  COMMON. 

D  NAME  IN  A  COMDECK  (CDIN1,  CDNU1,  OR  CDCOM). 

N  IN  A  NAMELIST  INPUT. 

Z  VALUE  IS  FROM  ZBLOCK  ON  (IN9),  DIRECT  OR  COMPUTED. 
_XXXXX_  DENOTES  PARAMETER  IN  COMDECK  CDIN1  OR  CDNU1. 

*0  BND(_NLBIG_)  -  ARRAY  FOR  ROW  OF  NEW  HYDRO  DATA. 

*D  BND2 (_NLBIG_)  -  ARRAY  FOR  ROW  OF  NEW  HYDRO  DATA. 

*N  CST ART  -  STARTING  CYCLE  TO  BE  MATCHED  ON  FILE  IN9. 

N  CSTOP  -  STOP  PROCESSING  CYCLES  AFTER  CYCLEI  >•  CSTOP. 

*  CYCLEI  -  PRESENT  CYCLE. 

*0Z  DELXK.IIBIG.)  -  VOLUME  CONTRIBUTION  IN  X  DIRECTION  FOR  CELL. 

*DZ  DEL Y I (_I JBIG_)  -  VOLUME  CONTRIBUTION  IN  Y  DIRECTION  FOR  CELL. 

•DZ  DELZ I <_IKBIG_)  -  VOLUME  CONTRIBUTION  IN  Z  DIRECTION  FOR  CELL. 
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c 

CDCOM  - 

c 

CDIN1  - 

c 

CDNU1  - 

c 

FIND  I  - 

c 

*01 

HYDROI!  I 

c 

IBND 

c 

* 

ICQNV  - 

c 

c 

c 

z 

IDIMEN  - 

c 

IEND 

c 

c 

*z 

IGEOM  - 

c 

D 

_ I IP I G_  - 

c 

*z 

UMAX  - 

c 

*1 

UMAX  1  - 

c 

0 

_I Jft I G_  - 

c 

*1 

UMAX  -  , 

c 

*1 

IJMAX1  - 

c 

0 

_IK8 I G_  -  , 

c 

*1 

IKMAX  -  1 

c 

*1 

IKMAX1  - 

c 

0 

_ILBIG_  - 

c 

♦N 

IN<5 

c 

*z 

INH  -  1 

c 

*z 

INHPB  -  1 

c 

N 

INIT  PR( 6) 

c 

N 

INI TRQ(6 ) 

c 

*N 

INM  -  1 

c 

IP 

c 

*1 

IROWPB  - 

c 

* 

ITIME  - 

c 

* 

ITY»E  - 

c 

*H 

JNM  -  1 

c 

JP 

c 

*N 

KNM  -  1 

c 

N 

LAS  T  PR( 6) 

c 

N 

L A$TRO<  6 ) 

c 

♦N 

NBND (6)  - 

c 

N 

NDIMEN  - 

c 

*N 

NFOUT ( 6 ) 

c 

N 

NGEOM  -  ■ 

c 

c 

D 

NIJKMAX. 

c 

D 

_NLBIG_  - 

c 

*N 

NNH 

c 

NTYPE  - 

c 

* 

NXPL  -  i 

c 

NXTBLK  - 

c 

NXTIM  - 

c 

c 

* 

NYPL 

c 

* 

NZPL 

c 

*N 

PROBIN  - 

c 

STRTIN  - 

c 

c 

STRTNU  - 

c 

c 

* 

TIMEIN  - 

c 

Til, T 12  - 

c 

*N 

TSH 

c 

*DZ 

XK.IIBIG 

c 

*D 

XNi .NIJKM 

c 

*N 

XN1 

COMOECK,  C OMHON/CDNU1 /»  PARAMETERS  &  RECIPIENT  ARRATS. 
SUBROUTINE  TO  LOCATE  INTERPOLATION  INDEX. 

LSI G_ )  -  ARRAY  FOR  1  BLOCK  OF  INPUT  HYDRO  VALUES. 

COUNT  ON  POSSIBLE  OUTPUT  FILES. 

CONVERSION  TYPE.  2  IS  CYLINDRICAL  TO  2D  CARTESIAN. 

3  IS  CYLINDRICAL  TO  30  CARTESIAN.  (3  ONLY  8/08/BA) 
COMBINES  IDIMEN,  IGEOM,  NDIMEN,  AND  NGEOH. 

DIMENSIONS  OF  OLD  HULL  RUN. 

SIGNAL  FROM  NXTIM.  0  MEANS  NEW  TIME  ON  FILE  (INO). 

1  MEANS  NORMAL  END  OF  IN<5,  -1  DENOTES  READING  ERROR. 
GEOMETRY  OF  OLD  HULL  RUN.  1  CARTESIAN,  2  CYLINDRICAL. 
ARRAY  SI2E.  NEED  IIBIG  >  UMAX. 

IMAX  FROM  DONOR  2-BLOCK. 

UMAX  ♦  1. 

ARRAY  SI2E.  NEED  IJ6IG  >  UMAX. 

JMAX  FROM  DONOR  2-BLOCK. 

UMAX  ♦  1. 

ARRAY  SIZE.  NEED  IKBIG  >  IKMAX. 

KMAX  FROM  DONOR  2-BLOCK. 

IKMAX  ♦  1. 

ARRAY  SI2E.  NEED  ILBIG  >■  INHPB 


IIMAX*IROWPB*INH 


NUMBER  OF  HYDRO  VARI ABLE S / POINT  IN  OLD  HULL. 

NUMBER  OF  HYDRO  VARI ABL ES / BLOCK  IN  OLD  HULL. 

-  TIME  DUMP  TO  START  TABULATION  OF  NFOUT (I ), 1*1,6. 

-  INITIAL  ROW  TO  START  TABULATION  IN  NFOUT ( I ), I»l, 6, 
NUMBER  OF  OUTPUT  POINT  IN  X  DIRECTION. 

POINTER  FOR  INTERPOLATION  IN  X. 

ROWS/BLOCK  IN  HYDROI. 

TIME  DUMPS  PROCESSED.  (REPEATED  FOR  EACH  BOUNDARY). 
INPUT  TYPE.  ^CYLINDRICAL,  2«2D  CARTESIAN,  3«3D. 
NUMBER  OF  Y  DIRECTION  POINTS  FOR  OUTPUT. 

POINTER  FOR  INTERPOLATION  IN  Y. 

NUMBER  OF  Z  DIRECTION  OUTPUT  POINTS. 

-  TIME  DUMP  TO  STOP  TABULATION  OF  NFOUT (I) , I-l, 6. 

-  LAST  ROW  TO  TABULATE  IN  NFOUTd ),  1*1,6. 

SIGNAL  FOR  BOUNOARY  OUTPUT.  0  IS  NO,  1  IS  YES. 


GEOMETRY  OF  NEW  HULL  RUN.  ONLY  1  (CARTESIAN)  NOW. 

MAY  ADD  2  FOR  CYLINDRICAL  LATER. 

-ARRAY  SIZE.  NEED  NIJKMAX  >-  MAX ( INM, JNM.KNM) . 

ARRAY  SIZE.  NEED  NLBIG  >■  M AX ( NXPL, NYPL, NZPL ) . 

HYDRO  VARIABLES/POINT  FOR  NEW  HULL  PROBLEM  (5  FOR  30). 
NEW  HULL  TYPE.  1-CYLINDR ICAL,  2»2D  CARTESIAN,  3»3D. 
NUMBER  OF  HYDROS/LINE  FOR  OUTPUT  X-PLANES  ( NNH* JNN ) . 
SUBROUTINE  TO  FIND  NEXT  BLOCK  OF  HYDRO  DATA  ON  (IN9). 
SUBROUTINE  TO  FIND  NEXT  TIME  DUMP  ON  FILE  IN9, 

BYPASS  Z-BLOCK,  READ  AND  PROCESS  VERTICES. 

NUMBER  OF  HYDROS/LINE  FOR  OUTPUT  Y-PLANES  ( NNH*INM ) • 
NUMBER  OF  HYDROS /LINE  FOR  OUTPUT  Z-PLANES  ( NNH* INM ) • 


CHECKS  ARPAYS  FOR  INPUT. 

SUBROUTINE.  READS  NAMELIST  INPUT  FOR  NEW  HULL. 
SETS  UP  POINTS.  CHECKS  OUTPUT  ARRAYS  SIZES. 


TIME  SHIFT  FROM  OLO  TO  NEW  HULL. 

_)  -  VERTICES,  THEN  MID-CELL  VALUES,  FROM  OLO  HULL. 

AX_ )  -  OUTPUT  POINTS  IN  X  DIRECTION. 

LEFT  OUTPUT  BOUNDARY.  MAY  BE  MINIMUM  X  FOR  OUTPUT  GRIDS 
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C  *N  XNM  -  RIGHT  EXTERNAL  CELL  CENTER.  NAY  BE  HAX  X  FOR  GRIOS. 

C  *N  XSH  -  SHIFT  IN  X  FROM  OLD  TO  NEW  HULL. 

C  *  XYZB ( 6 J  -  OUTPUT  BOUNDARY  LOCATIONS.  (  XN1, XNM. YNl.YNM.ZNl.ZNM) 

C  *DZ  YH.IJBIG.)  -  VERTICES,  THEN  NIDCELL  VALUES  FROM  OLD  HULL. 

C  *0  YNI.NIJKMAX.I  -  OUTPUT  POINTS  IN  Y  DIRECTION. 

C  *N  YN1  -  AFT  OUTPUT  BOUNDARY  (3D).  BOTTOM  BOUNDARY  (2D). 

C  MAY  BE  MINIMUM  X  FOR  OUTPUT  GRIDS  FOR  Y  AND  Z  PLANES. 

C  *N  YNM  -  EXTERNAL  CELL  CENTER.  FORE  BOUNDARY  3D.  TOP  BOUNDARY  2D. 
C  MAY  BE  MAXIMUM  Y  FOR  OUTPUT  GRIDS  FOR  X  OR  Z  PLANES. 

C  *N  YSH  -  SHIFT  IN  Y  FROM  OLD  TO  NEW  HULL. 

C  *DZ  ZK.IKBIG.)  -  VERTICES.  THEN  MIDCELL  VALUES.  FROM  OLD  HULL. 

C  *D  ZN( _NIJKMAX_  J  -  POINTS  FOR  OUTPUT  IN  Z  DIRECTION. 

C  *N  ZN1  -  BOTTOM  BOUNDARY.  MAY  BE  MINIMUM  Z  FOR  OUTPUT  GRIDS. 

C  *N  ZNM  -  TOP  EXTERNAL  CELL  CENTER.  NAY  BE  MAX  Z  FOR  OUTPUT  GRIDS 

C  *N  ZSH  -  SHIFT  IN  Z  FROM  OLD  TO  NEW  HULL. 

C  - - - 

C  - 

C 

♦CALL  CDIN1 
♦CALL  CDNU1 
♦CALL  COCOM 

C  SET  UP  PARAMETERS  FOR  INPUT  FILE 

NAMELIST  /INDAT1/  IN9, PRO B  IN, C ST ART, CSTOP 
C  INPUT  DEFAULTS 

IN9  -  9 
C START  •  0 
CSTOP  ■  9999999.0 
READ ( 5. INDAT1 ) 

WRITE  <6. INDAT1) 

CALL  STRTIN 

C  SET  UP  PARAMETERS  FOR  THE  NEW  HULL  DATA. 

CALL  STRTNU 

C  READY.  GO  TO  COOING  FOR  SELECTED  CONVERSION. 

IF (ICONV  .EQ.  2 ) GOTO  2000 
IFdCONV  .EQ.  3  )GOTO  3000 
WRITE(6.25)IC0NV 

STOP  •  ABORT  HULLUP  25.  NO  ICONV  CODING.' 

25  FORMAT!//'  ♦♦♦  ABORT.  HULLUP  25.  NO  CODING  FOR  ICONV  •  *.I2) 

C  - - 

C  CODING  FOR  ICONV  -2  NOT  YET  INSERTED. 

2000  CONTINUE 

WRITE(6.2005)IC0NV 

STOP  •  ABORT  HULLUP  2005.  NO  ICONV  CODING.' 

2005  FORMAT ( / / '  ♦♦  ABORT.  HULLUP  2005.  NO  CODING  FOR  ICONV  >  *.12) 

C  - 

C  ICONV  -  3.  INPUT  2D  CYLINDRICAL.  OUTPUT  3D  CARTESIAN. 

C  START  LOOP  ON  BOUNDARIES  FOP  IBND-1.6. 

3000  IBND  -  0 
3010  IBND  -  IBND  ♦  1 

IFUBND  .EQ,  7) GOTO  10000 
IF(NBND( IBND I  .GT.  0)G0T0  3020 
GOTO  3010 

C  START  LOOP  ON  TIME.  ITIME  IS  A  COUNTER. 

3020  ITIME  •  0 

REWINDIIN9) 

C  FIND  NEXT  TIME  DUMP  (INITIALLY  FIND  CSTART )  AND  VERTICES. 

3030  CALL  NXTIM(IEND) 

C  IF  END  OF  READABLE  FILE.  GOTO  END  OF  TINE  LOOP. 

IF( I  END  .NE.  0) GOTO  3900 
C  THE  INPUT  TIME  DUMP  HAS  BEEN  FOUND. 

C  VERTICES  READ,  CELL  CENTERS  AND  VOLUME  VALUES  SET. 

C  BRANCH  ON  BOUNDARY.  1,2  X-PLANE  3.4  Z-PLANE,  5.6  Y-PLANE. 

GOTO! 3050. 3050, 3500, 3500, 3750, 3750) , IBNO 
C 


non  o  o  rt  n  on o  n on 


C  XPLANE  BOUNDARY.  IBND  ■  1  OR  2. 

30S0  CONTINUE 

SET  UP  SOME  PARAMETERS. 

NN  •  ( J N-l ) *NNH  IS  LOCATION  FOR  OUTPUT  HYDROS. 

KN  IS  ROW*  JN  IS  COLUMN  IN  THE  OUTPUT  PLANE. 

IHB  IS  INPUT  HYDRO  BLOCK  NUMBER. 

JP  AND  JP2  ARE  THE  ROW  POINTERS  IN  THE  INPUT  PLANE. 

YI ( JP )  .LE.  ZN(KN)  .LT.  YIUP2).  JP2  ■  JP  ♦  1. 

IHB  -  0 
JP  ■  1 
JP2  •  JP*1 

START  LOOP  ON  ZN. 

DO  3140  KN»1, KNM 
NN*0 

ZP  >  ZN(KN) 

CALL  FINDI(JP*ZP,YI»IJMAX) 

ZR  IS  THE  INTERPOLATION  RATIO  FOR  Z  IN  YKJP)  TO  YI(JP*1). 
ZR  «  (ZP  -  YI(JP))/(YI(JP+1)-YI(JP)  ) 

IF( JP  .LT.  ( IH8-1)*IR0WPB  *1  .OR.  JP  .GT.  IHB*IR0WPB) 

♦  CALL  NXTBLK ( JP* IHB ) 

GET  READY  TO  INTERPOLATE  IN  XI  (FOR  F(R,YI(JP)). 

IP  •  1 

X-PLANE.  START  LOOP  ON  JN. 

XP  •  XYZB(IBND) 

XPSO  -  XP*XP 
DO  3000  JN«1,  JNM 
YP  -  YN(JN) 

YPSQ  ■  YP* YP 
RIN  ■  SORT ( XPSQ  ♦  YPSQ) 

CALL  FINDI(IP,RIN»XI»IIMAX) 

LI  •  ((JP-(IHB-1)*IR0WPB  -1)*IIMAX  ♦  IP-1)*INH 
RR  -  (RIN  -  XI(IP))/(XI(IP*1)  -  XI(  IP) ) 

FINO  F ( RIN*  XI ( IP  )  )  FOR  F  •  VR*VZ*E*RHO. 

RADIAL  VELOCITY*  AXIAL  VELOCITY*  SPECIFIC  ENERGY*  AND 
DENSITY, RESPECTIVELY.  DENSITY  IS  MASS/VOLUME. 

BND(NN+2>  •  HYDRO I( LI+2 )  ♦  (HYDROI ( L 1*2* INH )  -  HYDROI ( LI*2 ) ) *RR 

BND(NN+3)  •  HYDRO I<  L 1*3 )  ♦  (HYDROI ( L 1*3* INH »  -  HYDROI ( Ll*3 )) *RR 

BNO( NN+4 )  «  H YOROI(  L  1*4 )  •*  (HYDROI <  L  I*4*INH)  -  HYDROI ( LI*4> ) *RR 

RHOIP  -  HYDR0I(LI+5)/(DELXI(IP)*DELYI(JP)) 

R HOI  PI  -  HYDR0I(LI*5*INH)/(0ELXI( IP*1)*DELYI (JP) ) 

BND( NN*5)  -  RHOIP  ♦  (RH0IP1  -  RHOIP)*RR 
NN  •  NN  ♦  NNH 
3090  CONTINUE 

C  FINISHFD  INTERPOLATION  FOR  F(R*YI(JP))  FOR  A  ROW. 

3100  NN-0 

IP  •  1 

C  YKJP2)  IS  NOW  THE  INPUT  Y.  Y  GT  ZN.  (JP2«JP*1) 

C  IS  A  NEW  INPUT  HYDRO  BLOCK  NEEDED? 

IF  ( JP2  .GT.  IHB*IROWPB)C ALL  NXTBLK (JP2, IHB) 

C  READY  TO  INTERPOLATE  FOR  F(R»YI(JP2)»  THEN  FOR  F(R,Z). 

C  X-PLANE.  IBND  ■  1  OR  2.  LOOP  ON  YN. 

DO  3130  JN»1, JNM 
YP  -  YN(JN) 

YPSO  •  YP* YP 

RIN  -  SORT (XPSO  ♦  YPSO) 

CALL  FINDI(IP*RIN*XI*IIMAX) 

LI  •  ( ( JP2  -  (IHB-1 >*IROWPB  -1>*IIMAX  ♦  IP-1)*INH 
RR  •  (RIN  -  XI(IP))/(XI(IP*1J  -  X I ( IP) ) 

C  FIND  F(RIN»YI<0)»  THEN  F(RIN*ZN)  FOR  VX* VY* VZ* E* RHO. 

HVDRDP  •  HYDROI ( L 1*2 )  ♦  (HYDR0I(LI*2*INH)  -  HYDROI! LI+2) )*RR 
HYDROP  «  BND( NN*2 )  ♦  (HYOROP  -  BND(NN*2))*ZR 
C  HYDROP  •  VR.  FIND  AND  SAVE  VX  AND  VY  VELOCITY  COMPONENTS. 

IF(RIN  .GT*  0.0)G0T0  3103 
C  ZERO  RADIUS.  SET  VX  •  VY  ■  0.0. 


BND( NN+1)  -  0.0 
BND( NN*2 )  -  0.0 
GOTO  3104 

3103  BND(NN*1>  •  HYDROP*XP/RIN 
BNO( NN*2 >  -  HYDROP*YP/R IN 

3104  rlYOROP  «  H  YDROI  (  L  1*3 )  ♦  (  HYDRO  1  (  LI*3*INH)  -  HYDR0I(LI*3)  >*RR 
BND(NN*3)  •  BNO ( NN*3 )  ♦  (HYDROP  -  BND(NN*3) J *ZR 

HYDROP  •  HYDROI ( L  1*4 )  ♦  ( HYDROI ( L I*4*INH >  -  HYDROI C L I +4) ) *RR 
BND ( NN*4 )  •  BND(NN*4  )  ♦  ( HYDROP  -  BND(NN*4)>*ZR 
RHOI»  -  HYDR0I(LI*5)/(DELXI(IP>*DELYI(4P2>> 

RH0IP1  •  HYDR0I(LI*5*INH)/(DELXI(IP*1I*DELYI(JP2)> 

HYOROP  -  RHOIP  ♦  (RH0IP1  -  RHOIP )*RR 
BNO<  NN*5 )  -  BNO ( NN*  5 )  ♦  (HYDROP  -  BND(NN+5) )*ZR 
NN  *  NN  ♦  NNH 
3130  CONTINUE 

CAU  WRITPL(IBND» 

♦  NF3UT ( IBND ) , X P »  XS H, YN» 4NM, YSH, ZN» KNM,  ZSH, BND.NXPL ) 

3140  CONTINUE 

GOTO  3850 

YPLANE  BOUNDARY.  IBND  »  5  OR  6. 

3750  CONTINUE 

SET  UP  SOME  PARAMETERS. 

NN  -  (IN-1)*NNH  IS  LOCATION  FOR  OUTPUT  HYDROS. 

KN  IS  ROW.  IN  IS  COLUMN  IN  THE  OUTPUT  PLANE. 

IHB  IS  INPUT  HYDRO  BLOCK  NUMBER. 

4P  AND  JP2  ARE  THE  ROW  POINTERS  IN  THE  INPUT  PLANE. 

YI(JPJ  .LE.  ZN(KN)  .LT.  YI(JP2).  JP2  ■  JP  ♦  1. 

IHB  -  0 
JP2  «  4P*1 
4P  •  1 

START  LOOP  ON  ZN. 

DO  3840  KN-l.KNM 
NN  -  0 
ZP  •  ZN(KN) 

CALL  FINDI(4P»ZP.YI,I4MAX) 

ZR  IS  THE  INTERPOLATION  RATIO  FOR  Z  IN  YI(JP)  TO  YI(4P*1> 
ZR  *  (ZP  -  YI(JP))/(YI ( 4P*1 1—Y I ( JP) ) 

IF(JP  .LT.  ( IHB-l J*IROWPB  *1  .OR.  JP  .GT,  IHB*IROWPB) 

♦  CALL  NXTBLK (4  P»  IHB ) 

C  GET  READY  TO  INTERPOLATE  IN  XI  (FOR  F(R,YIIJP>). 

IP  •  1 

C  Y-PLANE.  START  LOOP  ON  IN. 

3770  YP  •  XYZBdBNOJ 
YPSO  -  YP* YP 
DO  3791  IN-1, INM 
XP  •  XN(  IN  ) 

XPSO  -  XP*X» 

RIN  -  SORT (XPSO  ♦  YPSO) 

CALL  FINDI(IP,RIN,XI,IIHAX) 

LI  -  ( (JP-(IHB-1)*IR0WPB  “1 >*I IHAX  ♦  IP-1)*INH 
RR  •  (RIN  -  XI<IP>>/(XI(I»*1>  -  XI(IPI) 

C  FIND  F ( RIN, XI ( IP ) )  FOR  F  ■  VR»VZ,E,RHO. 

C  RADIAL  VELOCITY,  AXIAL  VELOCITY,  SPECIFIC  ENERGY,  AND 

C  DENSITY, RESPECTIVELY.  DENSITY  IS  MASS/VOLUME. 

BND( NN*2 )  •  HYDR0I(LI*2J  *  (HYDROI ( L 1*2* INH)  -  HYDROI ( LI*2 )) *RR 

BND(NN*3)  -  HYDRO I( L I* 3 )  ♦  (HYDROI ( L 1*3* INH )  -  HYDROI ( L I ♦ 3 ) )*RR 

BND(  NN*4)  «  HYDROU  LI*4)  ♦  (HYDROK  LI*4*INN)  -  HYDROI ( L I ♦ 4 > ) *RR 

RHOIP  •  HYDROI(LI*5)/(DELXI(IP)*DELYI(4P>) 

RH0IP1  •  HYOR 01  (L 1*5*1 NH)/(DELXI(IP*1I*DELYI (JPI ) 

BND( NN+5)  -  RHOIP  ♦  (RH0IP1  -  RHOIP)*RR 
NN  •  NN  *  NNH 
3791  CONTINUE 
C 


FINISHED  INTERPOLATION  FOR  F(R,YI(JP>)  FOR  A  ROW 


■\_v~  VWV  V  Tib".  TV  ^ 


YKJP2J  IS  NOW  THE  INPUT  Y.  Y  GT  ZN.  NP2-JP+1) 

IS  A  NEW  INPUT  HYDRO  BLOCK  NEEDED* 

IF  (JP2  .GT,  IH8*IR0WPB)C ALL  N XTBLK ( JP2, IHB ) 

READY  TO  INTERPOLATE  FOR  F(R,YI(JP2),  THEN  FOR  F(R,Z>. 
LOOP  FOR  X-PLANE  OR  Y-PLANE* 

Y-PLANE.  IBND  •  5  OR  6.  LOOP  ON  XN. 

DO  3831  IN-1/ INN 
XP  •  XN(IN) 

XPSQ  -  XP*XP 

RIN  •  SORT ( XPSQ  ♦  YPSO) 

CALL  FINDI(IP*RIN,XI»IIMAX> 

LI  •  ( ( JP2  -  <IHB-1)*IR0WPB  -1>*IINAX  ♦  IP-i)*INH 
RR  -  IRIS  -  XI(IP))/(XI(IP*1)  -  XII  IP) ) 

FIND  F(RlN,YI(«-n,  THEN  F(RIN,ZN>  FOR  VX,VY, VZ,E,RHO. 
HYDROP  -  HYDRO!  (  LI*2  )  ♦  ( HYDRO  H  L  I*  2*INH  >  -  HYDR0l(LI*2)l*RR 
HYDPOP  -  8ND  <  NN  +  2 )  ♦  <  HYDROP  -  BND ( NN*2 ) ) *ZR 

HYOROP  -  VR.  FIND  AND  SAVE  VX  AND  VY  VELOCITY  COMPONENTS. 
IF  ( R  IN  .GT.  O.O) GOTO  3813 

ZERO  RADIUS.  SET  VX  -  VY  -  0.0. 

BND( NN  +  1)  -  0.0 
BND(NN*2>  •  0.0 


BND( NN  +  1)  -  0.0 
BND(NN*2>  •  0.0 
GOTn  3614 

BN0(NN*1J  •  HYOROP*XP/RIN 
BND(NN*2!  -  HYDROP*YP/R IN 

HYDROP  -  HYDROI <  L  1*3 )  ♦  ( HYDRQIlLIOdNH)  -  HYDROI! LI+3) )*RR 

BND( NN-3)  -  BND(NN+3)  ♦  <  HYDROP  -  BND(NN*3 ) ) *ZR 

HYDROP  •  HYDROI  ( L  1*4)  ♦  (HYDRDI(LI*4dNH)  -  HYDR0I(LI*4)  >*RR 

BNDt  NN»4)  -  BND ( NN+4 )  ♦  (HYDROP  -  BN 0< NN+4 J I *ZR 

RHOIP  -  HYDROI(LI»5)/(DELXI(IP»*DELYI(JP2I) 

R HOI PI  -  HYDR0I(LI+5+INH)/(DELXl(IP-l)*DELYI(JP2)) 

HYDPOP  -  RHOIP  ♦  (RH0IP1  -  RHOIP)*RR 

BND( NN*5)  -  BND(NN*5)  ♦  (HYDROP  -  BND(NN-5) )*ZR 

NN  -  NN  ♦  NNH 

CONTINUE 

CALL  WRITPLdBND. 

►  NFOUT(IBND),YP, YSH»XN/INM/X$H» ZN, KNM, ZSH, BND, NY PL) 

CONTINUE 

TINE  DUNP  FOR  AN  X-PLANE  OR  Y-PLANE  BOUNDARY  NOW  COMPLETE. 
CONTINUE 

ITIHE  -  ITINE  ♦  1 

DO  WE  WANT  ANOTHER  TINE  DUNPT 
IF(CYCL£I  .LT.  CSTOP )GOTO  3030 
GOTO  3900 

Z-PLANE  BOUNDARY,  IBND  >  3  OR  4.  (2D  CYLINDRICAL  TO  301 

ZP  -  XYZB(IBND) 

FIRST,  FIND  HYDRO  BLOCK  WITH  ROW  OF  DATA  FOR  YI(JP)  <  ZP. 

IHB  ■  0 
JP  •  1 

CALL  FINDI(JP,ZP,YI,IJMAX) 

FIND  INTERPOLATED  VALUES  AT  Z  -  ZP  ROW  BY  ROW. 

START  LOOP  ON  ROWS. 

FIRST,  FINO  VALUES  AT  Z  •  YI(JP). 


(2D  CYLINDRICAL  TO  3D) 


DO  3560  wN-l»JNM 
CALL  NXTBLMJP,  IHB) 


YP  -  YN(JN) 

YPSQ  •  YP*YP 

START  LOOP  ON  COLUMNS. 

IP  •  l 

DO  3510  IN-1,  INN 
XP  •  XN(IN) 

XPSQ  •  XP*XP 


1 


RIN  ■  SORT (XPSQ  ♦  YPSO) 

CALL  FINDKIP.RIN.XI.IIHAX) 

LI  -  ((JP-(IHB-l)*IROWPB  -1>*IINAX  *  IP-1)*INH 
RR  •  (RIN  -  XI(IP))/{XI(IP*1>  -  XI(IP)) 

C  FIND  F ( RIN*  XI ( IP ) )  FOR  F  «  VR.VZ.E.RHO, 

8ND(NN*2)  •  HYDROI(LI*2)  ♦  (HYDROK LI*2*INH>  -  HYDROI ( LI*2)> *RR 

BND(NN*3)  ■  HYDROI(LI*3)  ♦  (HYDRO I ( LI*3*INH)  -  HYDROI ( LI*3) I *RR 

BND(  NN*4)  -  HY0R0I(LI*4)  ♦  (HYDROI (  L 1*4*  INH I  -  HYDROI  (LI*4)  l*RR 

RHOIP  -  HYDROI(LI*5)/(DELXI(IP)*DELYI(JP) ) 

RHOIP1  ■  HYDROI(LI*5*INHW(DELXI(IP*l)*DELYI(JP)) 

BND(NN*5>  ■  RHOIP  ♦  (RHOIP1  -  RHOIP ) *RR 
NN  •  NN  ♦  NNH 
3510  CONTINUE 

:  NOW.  FIND  INTERPOLATED  VALUES  AT  Z  •  YIUP+1).  THEN  AT  ZP 

JP2  «  JP  ♦  1 

IF ( J  P  2  .GT.  IHB* IROWPB  )  CALL  NXTBLK ( JP2, IHB ) 

:  COMPUTE  INTERPOLATION  PATIO  IN  Z. 

ZR  •  (ZP  -  YI(JP2-1))/(YI(JP2)  -  YKJP2-1)) 

:  LOOP  ON  ROW. 

NN  -  0 

C  LOOP  ON  COLUMN. 

IP  -  1 

DO  3550  IN-1* INH 
XP-XN(IN) 

XPSQ  ■  XP*X* 

RIN  •  SORT (XPSQ  ♦  YPSO) 

CALL  FINDKIP.RIN.XI.IIHAX) 

LI  •  ( (4*2  -  ( IHB-1 )  *IRQWPB  -1)*IIHAX  ♦  IP-1)*INH 
RR  -  (RIN  -  Xl(IP))/(XI(IP*l)  -  XI(IP)) 

C  FIND  F(RIN»YI(*>)»  THEN  FCRIN.ZN)  FOR  VX.  VY.  VZ»  E»  RHO. 

HYDROP  •  HYDROI  ( L  1*2 )  ♦  ( HYDROI ( LI*2*INH)  -  HYDR0I(LI*2)  )*RR 
HYDROP  -  BND ( NN*2  J  ♦  (HYDROP  -  BND( NN*2 ) ) *ZR 
C  HYDROP  -  VR.  FIND  AND  SAVE  VX  AND  VY  VELOCITY  COMPONENTS. 

IF ( P I N  .GT.  0.0 ) GOTO  3543 
C  ZERO  RADIUS.  SET  VX  -  VY  •  0.0. 

BND(NN*1)  -  0.0 
BND(NN*2I  -  0.0 
GOTO  3544 

3543  BND(NN*1)  -  HYDROP*XP/RIN 
BND(NN*2)  •  HYDRQP*YP/RIN 

3544  HYDROP  -  HYDRQI(LI*3)  ♦  ( HYDROI ( LI*3*INH)  -  HYDROI ( LI*3 ) I *RR 
BND( NN*3)  -  BND(NN*3)  ♦  (HYDROP  -  BND(NN*3) )*ZR 

HYDROP  ■  HYOROI ( L !♦ 4 1  ♦  ( HYDROI ( L I* 4*INH)  -  HYDROI (L 1*4) )*RR 
BND( NN*4)  •  BND(NN*4)  ♦  (HYDROP  -  BND(NN*4) )  *ZR 
RHOIP  ■  HYDR0I(LI*5)/(DELXI(IP)*DELYI(JP2)) 

RH0IP1  •  HYDR0I(LI*5*INH)/(DELXI(IP*1)*DELYI(JP2)> 

HYDROP  •  RHOIP  ♦  (RH0IP1  -  RH0IP)*RR 
BND(  NN*5)  •  BND(NN*5>  ♦  (HYDROP  -  BND(NN*5)  )*ZR 
NN  >  NN*  NNH 
3550  CONTINUE 

CALL  WRITPLdBND. 

*  NFOUT ( IBND ) .ZP. ZSH. XN. INM.XSH.  YN»  JNH»  YSH. BND.NZPL ) 

3560  CONTINUE 

ITINE  •  ITIME  ♦  1 

C  00  WE  WANT  ANOTHER  TIME  DUMP? 

IF(CYCLEI  .LT.  CSTOP)GOTO  3030 
GOTO  3900 

C  END  OF  THE  TIME  LOOP  FOR  ICONV  •  3 

3900  GOTO  3013 

C  - 

C  END  OF  HULLUP.  PUT  TERMINAL  SIGNAL  ON  FILES. 

10000  DO  10010  1-1.6 
NF  «  NFOUT ( I ) 

IF(NBNDd)  .GT.  0)WRITE  (NF)  6  66..  666..  666.  .  666.  .  666.  666.  66b 


r>  r»  o  o 


C  TABULATE  A  PORTION  OF  EACH  B0UND9  FILE. 

IF(NBNDU)  ,GT.  0)CALL  REWRIT  (  I. NF  ) 

10010  CONTINUE 

WRITE (6# 10025)  ICONV 
STOP'  NORMAL  END  OF  HULLUP* 

10025  FORMAT ( •  NORMAL  END  OF  HULLUP.  ICONV  *',12) 

ENO 

SUBROUTINE  STRTIN 

CHECK  PROBLEM  NUMBER  ON  FILE  IN9  IS  PROBIN. 

READ  Z-BLOCK  AND  EXTRACT  INPUT  PARAMETERS. 

CHECK  THAT  INPUT  ARRAYS  ARE  LARGE  ENOUGH. 

♦  CALL  CD1N1 
♦CALL  CDCOM 
C 

DIMENSION  0101(41*  I Z ( 200  I* ZBL (200 ) 

EQUIVALENCE  (IZ.ZBL) 

C 

C  CHECK  THAT  INPUT  PROBLEM  NUMBER  IS  PROBIN. 

REWIND  !  IN9) 

NTRY  ■  0 

10  IF(NTRY  ,LT.  10JG0T0  20 
WRITE (6, 15)  IN9 

STOP'  ABORT  PROBLEM  AT  STRTIN  10.' 

15  FORMAT  C // •  ♦♦  ABORT.  INPUT  PROBLEM  AT  STRTIN  10.  FILE  TAPE', 14) 

20  NTRY  •  NTRY  ♦  1 

READ ( IN9) I DI Dili), 1-1,4) 

IF! EOF ( IN9)  .NE.  0) GOTO  10 

IF(  ABS  (OIDI  ( 2  )  -  PROBIN)  .LT.  O.OOOODGOTO  40 
WRITE (6*35)  IN9,DIDI(2), PROBIN 
35  FORMAT ( '  ♦♦  ABORT.  PROBLEM  NUMBERS  DIFFER.'/ 

♦  •  PROBLEM  NUMBER  ON  FILE'. 13.  •  IS  '.E15.6/ 

♦  •  INPUT  PROBLEM  NUMBER.  PROBIN.  IS  '.E15.b/> 

STOP*  ♦♦  ABORT.  STRTIN  30.  PROBLEM  NOS.  DISAGREE.' 

C  PROBLEM  NUMBERS  CHECK.  READ  Z-BLOCK.  SET  UP  IN-PARAMETERS. 

40  READIIN9)!ZBL1I), 1-1,200) 

DO  50  1-1,92 


IF! IZ ! I* 100 ) 

.EQ. 

6H0IMEN 

) 

IDIMEN 

* 

ZBL(Z) 

♦ 

0.5 

IF ( I Z ! 1+100 ) 

.EQ. 

6HGE0M 

) 

IGEOM 

m 

ZBLII) 

♦ 

0.5 

IF! I Z ( 1*100) 

.EQ. 

6HIMAX 

) 

UMAX 

m 

ZBLil) 

♦ 

0.5 

IFIIZIinOO) 

.EQ. 

6H4MAX 

J 

IJMAX 

m 

ZBLII) 

♦ 

0.5 

I F ( I 2 ( 1+100) 

.EQ. 

6HKHAX 

) 

IKMAX 

u 

ZBLII) 

♦ 

0.5 

IF! I Z ! 1*100 ) 

.EQ. 

6HNH 

) 

INH 

m 

ZBLII) 

♦ 

0.5 

IF ( I Z ( 1+100) 

.EQ. 

6HNR0WPB) 

IROWPB 

m 

ZBLII) 

♦ 

0.5 

50  CONTINUE 

INHPR  -  IIMAX*INH 
INHPB  -  INHPRPIROWPB 
I TYPE  •  IDIMEN 
IFUGEOM  . EQ.  2)ITYPE  •  1 
IIHAX1  -  UMAX  ♦  1 
IJMAX1  -  UMAX  ♦  1 
IKMAX1  -  IK.MAX  ♦  1 

C  CHECK  ON  THE  INPUT  ARRAY  LENGTHS 

IFIIIBIG.6E.IIMAX1  .AND.  I4BIG.GE.I4HAX1  .AND.  IKBIG.GE. IKN AX1 

♦  .AND.  ILBIG. GE . INHPB) GOTO  100 

WRITE ( 6,  55) I IBIG,  14 BIG, IKBIG,  ILBIG. I INAX1, 1 4MAX1, IKNAX1, INHPB 
STOP'  ♦♦  ABORT.  STRTIN  55.  DONOR  ARRAY* 

55  FORMAT!//'  ♦♦  ABORT.  STRTIN  55.  DONOR  ARRAY  TOO  SMALL. •» 

♦  /•  IIBIG,  14 BIG,  IKBIG,  ILBIG  -  ',4110, 

♦  /•  IIMAX1.Z4MAX1, IKMAXl, INHPB  >  ',4110) 

100  WRITE (6, 105) PROBIN 

RETURN 

105  FORMAT!/'  READY  TO  PROCESS  PR0B',F10.4> 

END 


no  * 


SUBROUTINE  STRTNU 

SET  JR  OUTPUT  CONTROL  PARAMETERS  FROM  DEFAULT  ANO  INPUT. 


C 
C 

♦CALL  C0NU1 
CALL  CDCOM 


NAMELIST/NDAT1/XN1.XNM* YN1* YNM*ZN1*ZNM* INM* JNM* KNM* TSH* XSH* YSH* ZSH 
NAMELIST  /NDAT2/NGE0M»NDIMEN*NNH*NBN0*NF0UT 

DEFAULT  VALUES  „  . 

NO  ACTIVE  BOUNDARIES.  OUTPUT  FILE  FOR  BOUNDARY  I  IS  10+1. 

DO  10  1*1*  6 
NBND(I)  *  0 
NFOUT(I)  «  10+1 
10  CONTINUE 

C  GEOMETRY  IS  CARTESIAN.  DIMENSION  IS  3. 

NGEOM  -  1 

NDIMEN  *3  ..  . 

C  DEFAULT  NUMBER  OF  NE V  HYDRO  VARIABLES/POINT  IS  5. 

NNH  *  5 

C  DEFAULT  BOUNDARIES  SET  TO  0.0.  DEFAULT  SHIFTS  SET  TO  0.0. 

XN1  *  0.0 
XNM  *0.0 
YN1  *0.0 
YNM  *0.0 
ZN1  -  0.0 
ZNM  *0.0 
TSH  •  0.0 
XSH  •  0.0 
YSH  *0.0 
ZSH  *0.0 

C  DEFAULT  OUTPUT  GRIO  DIMENSIONS  SET  TO  200. 

INM  -  200 
JNM  *  200 

C  KN”  READ°NDAT1  PARAMETER  CHANGES  AND  PRINT  PARAMETERS 
READ15.NDAT1) 

WRITE ( 6.N0AT1 I 

IFdNM.LT.O  .OR.  JNM.LT.O  .OR.  KNM.  LT.OIGOTO  200 
C  GRIO  INDICES  POSITIVE.  SET  EVENLY  SPACED  GRID. 

C  SET  GRID  VALUES.  FIRST*  FORCE  AT  LEAST  2  POINTS. 

IF < INM  .LT.  2)  INM  •  2 
IF( JNM  .LT.  2)  JNM  «  2 
IF<KNM  .LT.  2)  KNM  -  2 
DELP  •  ( XNM  -  XN1JMINH  -  1) 

XN(1 )  «  XN1 
DO  20  1-2* INN 
XN( I )  ■  XN(I-l)  ♦  OELD 
20  CONTINUE 

OELD  •  (YNM  -  YN1>/(JNN  -  1> 

YN( 1 )  •  YN1 
DO  30  J*2* JNM 
YN( J )  •  YN(J-1I  ♦  DELO 
30  CONTINUE 

DELO  -  (ZNM  -  ZNDMKNM  -  1) 

ZN( 1 )  -  ZN1 
DO  AO  K*2»KNM 
ZN(K )  ■  ZN(K-l)  ♦  DELO 
AO  CONTINUE 

2A5  IF( ABS(XN( INMI-XNMJ  .LT.  0.00001*XNMI  XN(INM>  •  XNM 
IF( ABS ( YN( JNM  I- YNM)  .LT.  0.00001*YNM)  YNCJNM)  •  YNM 
1F(ABS(ZN(KNM)-ZNM)  .LT.  0.00001*ZNM)  ZN(KNM)  •  ZNM 
C  PRINT  THE  GRID. 

WRITE (6* 2A6 )  (XN ( I ) *  I* 1* INM I 


WRITE ( 6>  24?)  ( YN ( J ) *  J* 1*  J NM ) 

WRITE ( 6»  2*8 )  ( ZN ( K) »  K» 1*KNH) 

GOTO  50 

246  FORMAT ( •  MESH  FOR  NEW  BOUNDARY  PLANES.  SHOULD  BE  CHECKED. 1 /» 

♦  •  NEW  CELL  CENTERS  MUST  BE  BETWEEN  EXTREMES.  RESULTS  BEST  IF'/* 

♦  •  CELL  CENTERS  ARE  AT  MESH  POINTS  (REMOVES  INTERPOLATION* • //» 

♦  lOX.'X  GRID'/»(1P6E15.?)) 

247  FQRMATI10X* ' Y  GRID* /* ( 1P6E15. ?> > 

248  FORMAT ( 10X* • Z  GRID' /»!1P6E15.7) ) 

50  RE AD( 5.NDAT2) 

WRITE (6* NDAT2 ) 

C  NTYPE  ■  1  FOR  20  CYINDRICAL*  ■  2  FOR  2D  CARTESIAN*  •  3  FOR  3D. 

NTYPE  -  NDIMEN  ♦  1  -  NGEOM 
C  SET  THE  CONVERSION  TYPE*  ICONV. 

C  ICONV  -  3  FOR  2D  CYLINDRICAL  TO  3D  CARTESIAN. 

ICONV  •  3*( ITYPE-1)  ♦  NTYPE 
C  SET  THE  VALUES  FOR  BOUNDARIES. 

XYZB ( 1 )  -  XN1 

XYZB ( 2  )  •  XNM 

XYZB (3)  -  ZN1 

XYZB ( 4  )  •  ZNM 

XYZB ( 5 )  •  YN1 

XYZB ( 6)  •  YNM 

C  30  OUTPUT.  COMPUTE  SIZE  OF  HYDRO  DATA/PLANE. 

NXPL  •  JNM*NNH 
NYPL  ■  INM*NNH 
NZPL  -  INM*NNH 

C  CHECK  THAT  ARRAYS  FOR  OUTPUT  ARE  LARGE  ENOUGH. 

IF(NI JKMAX.GE.INM  .AND.  NI JKMAX.GE. JNH  .AND.  NIJKNAX.GE.KNN 

♦  .AND.  NLBIG.GE.NAXO(NXPL*NYPL*NZPL))GOTO  100 
NAXPL  •  MAXO (NXPL»NYPL*NZPH 
WRITE(6*85)NIJKMAX*NLBIG* INM* JNM*KNM*MAXPL 

STOP'  *♦ ABORT.  STRTNU  85.  RECIPIENT  ARRAY  TOO  SMALL.' 

65  FORMAT!//'  **  ABORT.  STRTNU  85.  RECIPIENT  ARRAY  TOO  SMALL.'* 

♦  /•  NI JKHAX* INN* JNH* KNH» NAXPL  •  '*5110) 

100  WRITE ( 6*105) 

RETURN 

105  FORMAT! •  OUTPUT  PARAMETERS  OK.  READY  TO  COMPUTE  BOUNDARY  VALUES.') 
C  WILL  READ  DATA  AND  FORM  THE  RECIPIENT  GRID. 

200  INM  -  ABS(INM) 

JNM  *  ABS(JNM) 

KNM  -  ABS(KNM) 

C  FORM  XN ( I) *  I «1» INM 

READ ( 5*201 )  XN( 1 ) 

201  F0RMAT(E15.8,I5) 

NI  •  1 

212  READ ( 5*201 )  DXN*NNX 

IF ( NNX  .LT.  0 )G0T0  220 
NX  «  NUNNX 
NI  -  Nl'l 
DO  215  L«N1* N X 
XN! L )  •  XN(L-l)  ♦DXN 
215  CONTINUE 
NI  -  NX 
GOTO  212 

C  FORM  YN< J),4«1,JNM 

220  READ! 5*201)  YN(1) 

NI  -  1 

222  REA0(5»201)  DYN.NNY 

I F ( NNY  .LT.  0)G0T0  230 

NY  ■  NI  ♦  NNY 

NI  ■  Nl+1 

DO  225  L-N1*NY 

YN<L )  «  YN(L-l)  ♦  DYN 
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225  CONTINUE 
N1  •  NY 
GOTO  222 

C  FORM  ZN ( K) *  K«l* KNN 

230  READ ( 5*201 )  ZN(1) 

N1  -  1 

232  READ< 5*201)  DZN*NNZ 

I F € NNZ  .IT.  0 )60T0  240 
NZ  «  Nl+NNZ 
N1  >  Nl+1 
00  235  L-N1.NZ 
ZN(L)  •  ZN(L-l)  ♦  DZN 
235  CONTINUE 
N1  »  NZ 
GOTO  232 

C  RECIPIENT  GRID  FORMED.  CHECK  NX*  NY*  AND  NZ. 

240  IF( NX • EO. INN  .AND.  NY.E0.4NM  .AND.  NZ.EQ.KNM (GOTO  245 
WRITE  <  6»  241) NX* INM*  NY*  JNM*NZ*  KNM 

STOP'  ABORT.  STRTNU  240.  GRIO  INPUT  COUNTS  DISAGREE' 

241  FORMAT t 1  ABORT.  STRTNU  240.  NX.NE.INM  OR  NY.NE.JNM  OR  NZ.NE.KNM* 

♦  /, 'NX* INN, NY, JNM,NZ*KNH»'*3{2I5*5X>) 

END 

SUBROUTINE  NXTIN(IEND) 

FIND  NEXT  INPUT  TIME  DUMP. 

IF  ITIHE  ■  0*  FIND  CYCLE  WITH  DIDK  3)  •  CSTART. 

SET  IENO*CYCLEI*  TIME  IN. 
lEtID  >  0*  DATA  AVAILABLE. 

IENO  -  1,  NORMAL  END  OF  DATA.  DIDI(I)  •  666.0 
IEND  ■  -1*  COULD  NOT  FIND  NEXT  TIME  DUMP. 

BYPASS  THE  Z-BLOCK 

READ  VERTEX  RECORDS*  FIND  VOLUME  CONTRIBUTIONS  ANO  CELL 
CENTERS. 


♦CALL  CDIM1 
♦CALL  COCOM 
C 

DIMENSION  DIDI(4) 

C 

IEND  •  0 
1  NTRY  -  0 

10  READ ( IN9) (DIDI(I)*I>1*4) 

IF l EOF ( IN9  )  .NE.  0)GOTO  20 

C  CHECK  FOR  NORMAL  FILE  TERMINATION  WITH  DIDI ( I ) -666.0 . 

IF ( ABS ( DID  I ( 1 )  -  666.0)  .LT.  0.001 

♦  .AND.  ABS ( 0 ID  I ( 2)  -  666.0)  .LT.  0.001)  GOTO  30 

C  CHECK  FOR  START  OF  TIME  DUMP. 

IFUBStOIDIU)  -  555.0)  .GT.  0.001 

♦  .OR.  ABS  (0101(2)  -  PPOBIN)  .GT.  O.OOOODGOTO  10 

C  WE  HAVE  READ  A  HEADER  RECORD  FOR  A  HULL  TIME  DUMP. 

IF( ITIME  .GT.  OIGOTO  40 
C  INITIATION.  FIND  STARTING  CYCLE. 

IF(ABS(DIDI(3>  -  CSTART)  .LT.  O.DGOTO  40 
C  NOT  STARTING  CYCLE. 

IFIDIDII3)  -CSTART  .LT.  0.0)G0T0  1 
C  REQUESTED  STARTING  CYCLE  NOT  PRESENT. 

WRITE  <6*15)0 ID I<3)*CST ART 
IEND  ■  -1 

STOP'  ♦♦  ABORT  IN  NXTIM,  CANNOT  FINO  CSTART.' 

15  FORMAT!//'  ♦♦  HAVE  READ  CYCLE'*I5*'  WITHOUT  FINDING  CSTART  • ' , 15 
20  NTRY  ■  NTRY  ♦  1 

IFCNTRY  .LE.  10)G0T0  10 
C  CANNOT  READ  ANY  MORE  INPUT. 

IENO  ■  -1 
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RETURN 


C  NORMAL  END  OF  HULL  RESTART  FILE. 

3C  IENO  •  1 
RETURN 

C  HEADER  RECORD  READ.  START  PROCESSING. 

40  CYCLEI  •  DIDI  (3  I 
TIME  IN  ■  DIDI  (4) 

C  BYPASS  THE  Z- BLOCK 

R£AD( IN9IDUNNY 

C  READ  VERTEX  RECORDS. 

IFIIGEOM  .EQ.  2) GOTO  100 

C  3D  VERTEX  INPUT  PURPOSELY  OMITTED  UNTIL  LATER. 

WRITE 1 6#  45 ) IGEOM 

STOP'  **  ABORT.  STRTIM  45.  WRONG  IGEOM.' 

45  FORMAT ( // '  **  ABORT,  STRTIM  45.  NO  CODING  FOR  IGEOM  ••#12) 

C  READ  20  VERTICES.  CHECK  SUITABILITY. 

100  READ (IN9)tXI(I)#Ia2» 1IMAX1 ) #(YI(JI#J«1#IJMAX1) 

XI ( 1 )  >  0.0 

C  CHECK  THAT  INPUT  SPACE  INCLUDES  OUTPUT  SPACE. 

IF ( I CONV  .EQ.  3IG0T0  110 
X F I ICONV  .EQ.  6 ) GOTO  120 
C  INPUT  ANO  OUTPUT  BOTH  2D. 

IF ( XN1  . GE •  X I < 1 )  .AND.  XNM  , LE.  XI(IIHAX) 

♦  .AND.  YN1  .GE.  YI ( 1 )  .AND.  YNM  .LE.  YI ( I J MAX ) ) GOTO  140 
GOTO  130 

C  ICONV  •  3.  2D  CYLINDRICAL  INPUT#  3D  CARTESIAN  OUTPUT. 

110  DSOMAX  ■  AMAX1( XN1*XN1  ♦  YN1*YN1»  XN1*XN1  *  YNM* YNM# 

♦  XNM*XNM  ♦  YN1*YN1»  XNM*XNM  *YNH*YNH ) 

IF(OrONAX  .LE.  XI (IIMAX)**2 

*  .AND.  ZN1  .GE.  YK1I  .AND.  ZNM  .LE.  YI  ( I JMAX )  IGOTO  140 
GOTO  130 

C  ICONV  •  6.  2D  CARTESIAN  INPUT  AND  3D  CARTESIAN  OUTPUT. 

120  IF(  XN 1  .GE.  X  X  f  1 )  .AND.  XNM  .LE.  XIIIIMAX) 

*  .AND.  ZN1  .GE.  YI(1)  .AND.  ZNM  .LE.  YI ( UMAX) IGOTO  140 
GOTO  130 

C  INPUT  GRID  DOES  NOT  INCLUDE  OUTPUT  BOUNDARIES#  ABORT. 

130  WRITE ( 6# 135) 

STOP'  **  ABORT.  NXT1M  135.  GRID  INADEQUATE.' 

135  FORMAT  (•  **  ABORT.  NXTIH  135.  OUTPUT  OUTSIDE  DONOR  BOUNDARY.') 

6RID  OK.  SET  DELYI(J)  •  YIJ*1)-YU>. 

SET  OELXKI)  ■  VOLUME  CELL  II#J  )  FDELYI  ( J  )  . 

SET  X(I)  AND  Y(J)  TO  CELL  CENTERS. 

140  DO  150  I* 1# 1 1  MAX 

DELXI(I)  •  XIII*1>  -  mn 
XIII)  -  0. 5* ( Xl ( 1*1 )  ♦  XIII)) 

IF!  IGEOM  .EQ.  DGOTO  150 

INPUT  GROMETRY  IS  2D  CYLINDRICAL. 

OELXIII)  >  6. 2831853*XI 1 1 )*DE L XI II) 

150  CONTINUE 

DO  160  J*lt  UMAX 
DELYIIJ)  •  Y 1 1  J*1 )  -  YIIJ) 

YIIJ)  •  0.5*1 YI l J)  *  YI I J*l)) 

160  CONTINUE 
RETURN 
END 

SUBROUTINE  NXTBLK I J P» IHB I 

THIS  MOVES  THE  DESIRED  BLOCK  OF  INPUT  HYDRO  DATA  INTO  HYDROI. 

THIS  SUBPROGRAM  ASSUMES  DONOR  HULL  IS  2D. 

IHB  IS  THE  NUMBER  OF  THE  PRESENT  BLOCK. 

JP  IS  THE  DESIRED  ROW. 

♦CALL  CDIN1 
♦CALL  CDCOM 
C 


G2 


onoonnnnnoonononnnon  o  r>  o  r>  o  n  o  o  n  nonnn 


IROWBB  IS  INPUT  ROWS  PER  BLOCK. 

HYDROI  IS  THE  ARRAY  FOR  THE  INPUT  HYDRO  DATA. 

INHPB  IS  THE  HYDRO  VARIABLES  IN  A  BLOCK. 

JLAST  IS  THE  LAST  ROW  IN  THE  PRESENT  BLOCK. 

JLAST  •  IROWPBMHB 

IS  JP  IN  A  HIGHER  BLOCK? 

3FIJP  .GT.  J L AST ) GOTO  30 

IS  JP  IN  THE  PRESENT  BLOCK? 

IF< JP  .GT.  JL AST-IRQWP8  JRETURN 

HOVE  TO  BEGINNING  OF  PRESENT  BLOCK. 

BACKSPACE  ( IN9) 

HOVE  TO  BEGINNING  OF  BLOCK  IHB-I. 

10  BACK$PAC£< IN9) 

IHB  •  If!B-l 

JLAST  •  JLAST  -  IROWPB 

IS  ROW  JP  IN  A  LOWER  BLOCK? 

IF( JP  .LE.  JLAST-IROWPBJGOTO  10 

JP  IS  IN  THIS  BLOCK.  READ  HYDRO  DATA. 

20  READ  ( IN9H  HYDRO  I  (L>»L«1#  INHPB) 

RETURN 

CHECK  FOR  ROW  JP  IN  THE  NEXT  HIGHER  BLOCK. 

30  IHB* IHB+1 

JLAST  -  JLAST  ♦  IROWPB 

IS  JP  IN  THIS  BLOCK? 

IFIJP  .LE.  JLASTJGOTO  20 

JP  GT  JLAST t  SO  BYPASS  HYDRO  BLOCK. 

READ! IN9)DUHHY 
GOTO  30 
END 

SUBROUTINE  WRITPL (IBND»NF»X/XS»Y»NY»YS»ZtNZ»ZS# BOUND. NL ) 

THIS  ROUTINE  STORE S  THE  OUTPUT  DATA  FOR  A  PLANE  IN  A  FILE. 
STORE  1  ROW  AT  EACH  ENTRANCE. 

TEHPOR ARY  FORHAT.  FINAL  STORAGE  IN  REWRIT.  (8/OS/BA) 
THE  NOTATION  IS  LIKE  AN  X-PLANE#  BUT  ROUTINE  IS  FOR  ALL  3. 
X-PLANE  OROER  (X> Y«Z ) 

Y-PLANE  ORDER  (Y,X,Z) 

Z-PLANE  ORDER  (Z*X# Y) 

IBND  IS  FILE  COUNT  USEO  HERE  TO  TELL  INITIAL  ENTRY. 

NF  IS  THE  FILE  NUHBER 
X  IS  LOCATION  OF  THE  PLANE. 

XS  IS  THE  SHIFT  IN  THE  PLANE  LOCATION. 

Y  IS  THE  ARRAY  OF  COLUHN  POSITIONS  IN  THE  PLANE. 

NY  IS  THE  NUHBER  OF  COLUHNS. 

YS  IS  THE  SHIFT  IN  COLUHN  POSITIONS. 

Z  IS  THE  ARRAY  OF  ROW  POSITIONS. 

NZ  IS  THE  NUHBER  OF  ROWS. 

ZS  IS  THE  SHIFT  IN  ROW  POSITIONS. 

BOUND  IS  THE  LOCATION  OF  THE  HYDRO  DATA  FOR  THE  PLANE. 

NL  IS  THE  AHOUNT  OF  HYDRO  DATA  FOR  THE  PLANE. 

•CALL  CDCOH 

DINE  NS  ION  Y(1)»Z(1)»B0UND(1J 
DIMENSION  TBNDI 6) 

DATA  TBND/ 6*-999999, 999/ 

C 

C  SHIFT  DATA  FOR  OUTPUT 

TIHEIN  -  TIHEIN  -  TSH 
X  •  X  -  XS 

IF ( YS  .EQ.  0.0) GOTO  30 
DO  20  J*l» NY 
Y<J)  ■  Y(J)  -  YS 
20  CONTINUE 

30  IF( ZS  .EQ.  O.OIGOTO  A5 
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Z  (K )  •  Z  (K)  -  ZS 
40  CONTINUE 

4$  IFtTIMEIN  .EQ.  TBNO ( IBND)  ) GOTO  50 
TBND(IBND)  -  TIMEIN 

WRITE (NF)  55 5.0, PROBIN, TIME  IN, X,NNH,NY,NZ 
WRITE (NF)  (YO), J«1,NT),(Z(K),K-1,NZ) 

50  WRITE(NFI  < BOUND (  L) ,  l- 1, NO 
C  SHIFT  OAT A  BACK. 

TIHEIN  •  TINEIN  «■  TSH 
X  •  X  ♦  XS 

IF(YS  .£0.  0. 0) GOTO  70 
DO  60  J>1,NY 
YU)  -  Y ( J  )  ♦  YS 
60  CONTINUE 

70  IF ( ZS  .EQ.  O.OJGOTO  90 
DO  80  K-1,NZ 
Z (K )  •  Z(K )  ♦  ZS 
00  CONTINUE 
90  RETURN 
ENO 

SUBROUTINE  FINDII I, ZP, YI, IN) 

FINO  I  SUCH  THAT  I  •  1  IF  ZP  IT  YI(2I 

I  •  IN-1  IF  ZP  GE  YI ( IN-1 ) 

OR  YI< I)  <  ZP  <-  YI(I*1>. 

I  IS  THE  CURRENT  POINTER. 

ZP  IS  A  VARIABLE  (FOR  3D  OUTPUT,  THE  LOCATION  OF  A  PLANE), 
rl  IS  A  TABLE  (AN  ARRAY  OF  CELL  CENTERS  FRON  OLD  HULL  RUN). 
IN  IS  THE  NUNBER  OF  ENTRIES  IN  THE  TABLE. 

DIMENSION  rill) 

C 

C  FIRST,  CHECK  THAT  0  <  I  <  IN. 

I F ( I  .LT.  1  .OR.  I  . GE .  IN-1 ) I *1 
C  IS  ZP  >  YKDI 

IFIZP  .GT.  YI (I) ) GOTO  20 

C  Z®  <■  Till).  CHECK  FOR  I  «  1,  THEN  REDUCE  I. 

10  I F ( I  .EQ.  1) GOTO  30 
I  -  I  -  1 

C  IS  ZP  <«  YI(I)T 

IFIZP  .LE.  YI (I ) ! GOTO  10 
GOTO  30 

C  ZP  »  YI  (I )  •  CHECK  FOR  I  •  IN-1. 

20  IF( I  .EQ.  IN-1) GOTO  30 
C  CHECK  FOR  ZP  <■  YIII*1>. 

IFIZP  .LE.  YI ( 1*1 )) GOTO  30 
I  -  I  ♦  1 
GOTO  20 

C  DESIRED  I  FOUND.  RETURN 

30  RETURN 
END 

SUBROUTINE  REWRIT(I8ND,NF ) 

C  SUBROUTINE  TO  COMBINE  TINES  AND  PRODUCE  FINAL  B0UND9  FILES. 

♦CALL  CDNU1 
C 

C  WRITTEN  AS  IF  FOR  AN  X  PLANE,  BUT  FOR  ANY  3D  B0UND9  PLANE. 

DIMENSION  INITPR(6)»LASTPR(0),INITRO(6),LASTRO(6) 

NANEL IST/T  ABDAT /INITPR, LASTPR, IN1TR0,LASTR0 

DATA  ISW, I NIT PR, LAST PR, INITRO, LASTR0r0»6*i»6*0, 6*1, 6*10/ 

IF( ISW  . GT .0 ) GOTO  10 

C  INITIAL  ENTRY.  READ  TABULATION  CONTROLS. 

ISW  •  1 

READ! 5,TA8DAT) 


WRITE  1 6, TABO AT) 

READ  FILE  NF.  PUT  IN  FILE  1  AND/  SKIPPING  FIRST  DUMP/  IN  2 
10  REWIND  NF 
REWIND  1 
REWIND  2 
SW  «  0.0 
WRITE (6/110)NF 

110  FORMAT ( • 1  BEGIN  COMBINING  TIMES  FOR  FILE'/I5) 

20  RE AO (NF I  HE  ADI /PROBIN/ Til/ X/NNH/NY/NZ 
IF(E  OF  (NF  )  .HE.  0)G0T0  <90 
IF(ABS(H£ADl-666.0)  .LT.  O.OODGOTO  AO 
IF ( ABS(HEADl-555.0)  .GT.  O.OODGOTO  91 
WRITE (1)  HEA01/ PR  OB  IN, Til /X/NNH/NY/NZ 
IF ( S W  .GT.  0.0)WRITE(2)  HEAD1/ PR0BIN/TI1 
R£AD( NF) (YN( J)/ J«1/NY), ( ZN (K) / K*l/N Z > 

NL  -  NNH*NY 
DO  30  K-l/NZ 

READ ( NF )  ( BNO (L ) / L* 1/ NL  ) 

WRITE(1)(BN0(L)/L«1/NL) 

IF ( S W  .GT.  0.0)  WRITE ( 2 )  ( BND ( L ) / L- 1 /NL ) 

30  CONTINUE 
SW  «  1.0 
GOTO  20 

FILE  NF  HAS  BEEN  READ  AND  STORED  IN  FILES  1  AND  2. 

NOW  WE  COLLATE  DUMPS  WITH  PROPER  HEADINGS. 

40  WRITE (1)666.  0/666. 0/666.0/ 666. 0/666/666/666 
WRITE (2)666.0/666.0/ 666.0 
REWIND  1 
REWIND  2 
REWIND  NF 
NOUT  -  1 

50  READ ( 1 )  HEAD1. PROBIN. Til/ X/NNH/NY/NZ 
IF(ABS(HEADl-666.0>  .LT.  O.OODGOTO  92 
IF  (  ABS  ( HE  ADI  -  555.0)  .GT.  O.OODGOTO  92 
READ ( 2)  HE AD2/DUM/TI2 
IF(ABS(HEAD2-666.0)  .LT.  O.OODGOTO  70 
IF(HE ADI  .NE.  HEAD 2  .OR.  DUM  .NE.  PROBIN)GOTO  93 
WRITE(NF)HEAD1/PR0BIN/TI1/TI2/ X/NNH/NY/NZ 
IF( IN  I TPR ( IBNO)  .LE.  NOUT  .AND.  LAS  TPR ( IBND )  .GE.  NOUT) 

♦  WRITE (6, 900)  HEAD1/ PROBIN, Til, TI2,X,NNH,NY/NZ 
WRITE (NF)  (YNUI/J-l/NY),  (ZN(K),K-1,NZ> 

IF(INITPR( IBNO)  .LE.  NOUT  .AND.  LASTPR(IBND)  .GE.  NOUT) 

♦  WRITE (6/901 1  (YN(J)/4*l/NY)/(ZN(K)/K»l/NZ) 

DO  60  K*lz  NZ 

READ ( 1 )  (BND(L)zL-lzNL) 

RE  AD ( 2 )  ( BND2 (L ) / L* 1/ NL ) 

WRIT E(NF)( BNO (L )/ L- 1/NL )/ ( BND2 U ) , L -1/NL > 

IF ( IN ITPR ( IBND)  .LE.  NOUT  .AND.  LASTPR(IBND)  .GE.  NOUT  .AND. 

♦  INITRO(IBND)  .LE.  K  .AND.  LASTROI IBND )  .GE.  K) 

♦  WR I TE ( 6/ 902  )  K, (BND ( L >, L-l/N L >, ( BND2 ( L > , L- 1/NL ) 

60  CONTINUE 

NOUT  ■  NOUT  ♦  l 
GOTO  50 

FILE  NF  NOW  STORED  IN  FINAL  FORM  EXCEPT  FOR  TERMINAL  RECORD 
70  WRITE (NF )  666.0/666.0/666.0/666.0/666.0/666/666/666 
RETURN 

ABORT  PRINTS 

90  WRITE (6,190)NC 

STOP'  ABORT.  REWRIT  90.  NO  DATA  ON  FILE.' 

190  FORMAT ( / / '  ABORT.  REWRIT  90.  NO  DATA  ON  FILE', IS) 

91  WRITE (6/191 )NF, HE  ADI /PROBIN/ Til/ X/NNH/NY/NZ 
STOP'  ABORT.  REWRIT  91.  MIX  UP  IN  FILE  NF.' 

191  FORMAT (// '  ABORT.  REWRIT  91.  HIXUP  IN  FILE  NF.',I5/ 

♦  '  HEAD1/PR0B IN/ Til/ X/NNH/NY/NZ  • » /, 1P4E 15 .7/ 3 110 ) 


92  IT  ■  1 
GOTO  94 

93  IT  ■  2 

94  WRITE (6*194) IT# HE  ADI » PR08IN#TI1#X#NNH#NY»NZ» HE AD2»DUM#TI2 
STOP*  ABORT.  REWRIT  94.  NIXUP  IN  FILE  1  OR  2.' 

194  FORMAT ( • ABORT .RE WRIT  94.  NIXUP  IN  FILES  1  OR  2.  IT»'#I5/# 

♦  •  HEA01#P ROBIN# Til# X#NNH#NY#NZ  •' /# 1P4E1 5. 7* 311 0/» 

♦  •  HE AD2* DUN  »T I 2  »*#1P3E15.7> 

900  FORMAT (20X# • HEADER* //IP 5E 15. 7/ 3110) 

901  FORMAT! 10X# 'VERTICES'  /#  (1P6E15.7) ) 

902  FORNATUOX# 'PART  OF  HYDRO  DATA  STARTING  AT  ROW  •»  15/#  UP5E15.7  J  > 


APPENDIX  B 


LISTINGS  FOR  KEEL  AND  HULL 

Appendix  B  has  listings  of  the  runstreams  for  KEEL  and  HULL  for  problem 
8U i6 . 20 .  The  listing  of  the  HULL  runstream  is  followed  by  listings  of  the 
correction  file  for  HULL,  file  HULCORR,  and  the  changes  for  B0UND9  input,  file 
NEWBND9.  These  two  files  are  in  MFA. 

Note  that  file  HULL81405P09,  the  same  file  that  was  used  in  HULLUP,  is 
attached  for  initial  data  in  the  KEEL  run.  (This  is  actually  only  an 
illustrative  gesture  for  this  problem.  The  entire  3-D  region  was  ambient 
initially.)  Note  also  in  the  HULL  runstream,  that  files  for  the  four  B0UND9 
input  sides  are  attached  to  appropriate  local  files  and  the  MFA  files  NEWBND9 
and  HULCORR  are  included  for  UPDATE.  These  attaches  and  the  declaring  of 
B0UND9  boundaries  and  the  GENERATE  AIR  FIREIN  HULL  sequence  in  KEEL  are  the 
only  additional  input  at  this  stage.  A  lot  of  planning  is  needed  to  set  up 
the  HULLUP  run  to  produce  the  desired  input  in  the  files  for  the  BOUND9 
boundary  i nput . 

KEEL  Changes  for  B0UND9  Boundaries 

There  were  no  changes  needed  to  input  for  KEEL  for  the  new  type  9 
boundaries.  A  few  minor  changes  were  made  in  subroutine  HULLIN  to  add  some 
problem  aborts  with  explanatory  printing  and  to  use  the  ambient  atmosphere  for 
normalizing.  If  a  region  from  a  2-D  cylindrical  run  is  to  be  mapped  into  a 
3-D  (or  new  2-D)  region,  the  restart  file  from  the  donor  run  must  be  attached 
as  TAPE  UL  and  the  GENERATE  AIR  FIREIN  HULL  input  sequence  must  be  used  for  a 
region  that  includes  all  cells  in  the  new  space  that  overlap  the  non-ambient 
part  of  the  2-D  donor  space  at  the  initiation  time.  The  initiation  time  is 
the  first  time  on  the  donor  restart  file  greater  than  the  input  time  T.  T  is 
replaced  by  this  time.  The  user  may  set  the  boundary  conditions  in  KEEL,  but 
they  have  no  affect  until  the  HULL  run. 

If  one  must,  the  reference  location  for  the  new  run  can  be  shifted  through 
input  of  XOB,  YOB,  and  HOB  in  KEEL.  Time  cannot  be  shifted  in  KEEL  without 
some  recoding. 

HULL  Changes  for  B0UND9  Boundaries 

The  coding  changes  to  use  the  B0UND9  boundaries  for  3-D  HULL  runs  is  in 
the  MFA  file  NEWBND9,  ID=JDW.  This  coding  may  be  made  a  permanent  part  of  the 
HULL  file  when  the  check-out  for  these  boundaries  is  completed.  If  one  wishes 
to  use  this  coding  with  a  3-D  HULL  run,  file  NEWBND9  must  be  picked  up  in  the 
SC0PE2  runstream  and  READ  as  part  of  the  UPDATE  input.  The  user  must  also 
have  available,  and  attach,  an  input  file  for  each  B0UND9  boundary.  Any  of 
the  6  boundary  planes  may  be  declared  a  type  9  boundary.  The  NEWBND9  coding 
assumes : 

TAPE31  -  Left  boundary  at  X0 

TAPE32  -  Right  boundary  at  X^^  -  1/2  DX^  ^ 

TAPE33  -  Bottom  boundary  at  ZO 


TAPE j4  -  Top  boundary  at  2^^  -  1/2  D2RMAX 

lAPLio  -  Aft  boundary  at  YO 

TAPcob  -  Fore  boundary  at  ¥lljA„  -  1/2 

vJMAA  JMAX 

Notice  that  the  planes  of  imposed  data  for  the  low  index  boundaries  are  on 
the  boundaries  and  the  planes  of  imposed  data  for  the  high  index  planes  are 
through  the  centers  of  the  external  plane  of  cells. 

The  data  on  these  files  is  a  series  of  "dumps'*  of  unformatted  data  for 
successive  time  intervals.  Eaoh  dump  is  a  header  record,  a  grid  reoord,  and 
the  hydro  data  at  both  ends  of  the  time  interval  one  row  per  reoord. 

The  header  record  for  an  X-plane  is  8  numbers,  5  reals  and  i  integers: 

HEAD  1 ,  OLPROB,  Til,  T12,  XIN,  NNH,  JN,  KN 

where, 

HEAD1  =  555.0 

QLPKOB  =  The  donor  problem  number 

'i'll  s  The  first  time 

Il2  s  The  last  time 

XXN  s  X  value  for  the  plane 

NNH  s  5,  the  number  of  hydro  values  per  point 

JN  =  The  number  of  columns  of  data 

KN  =  The  number  of  rows  of  data. 

The  grid  record  for  an  X-plane  is: 

(YIN(J),  J=1 , JN) ,  (ZIN(K),  K=1 ,KN) . 

This  defines  a  JN  x  KN  grid  of  points  on  the  X-plane  at  which  hydro  data 
will  be  defined. 

This  is  followed  by  KN  records  of  hydro  data  for  K=1 ,2, . . . ,KN: 

(((HB(L,J,N),N=1 ,NNH) , J=1 , JN) ,L=1 ,2) . 

The  5  hydro  values  are  in  order 

U  -  The  x  velocity  component  (cm/sec) 

V  -  The  Y  velocity  component  (cm/sec) 

w  -  The  L  velocity  component  (cm/sec) 

E  -  internal  specific  energy  (ergs/g) 

P  -  Density  (g/cnr) 

The  5  hydro  values  are  given  for  J=1 ,2, . . . , JN  for  time  Til  (L=1)  and  then 
for  time  TI2. 

For  a  Y-piane,  the  boundary  plane  location  is  YIN  and  there  are  IN  columns 
of  data  at  XIN(i),  1=1,  IN  on  KN  rows  at  ZIN(K),  K=1,  KN. 


for  a  2-plane,  the  boundary  plane  is  at  21N  and  there  are  IN  columns  of 
data  at  XxN(I),  1=1,  IN  on  JN  rows  at  YIN(J),  J=1,  JN. 

The  program  Interpolates  In  this  input  data  in  time  and  space  on  the 
boundary  plane.  Any  attempt  to  extrapolate  in  either  time  or  space  causes  a 
program  abort.  The  program  also  aborts  if  the  plane's  location  does  not  agree 
witn  the  requested  position. 

There  is  no  further  new  input  other  than  declaring  the  type  9  boundaries. 
There  are  a  few  coded  in  values  that  might  be  changed:  Arrays  are  set  for  a 
maximum  of  200  words  for  IN,  JN,  and  K.N,  and  4000  words  for  the  hydro  data. 
This  uses  400  words  of  SCM  and  4000  words  of  LCM  for  each  BGUND9  boundary.  It 
might  be  necessary  to  increase,  or  decrease,  these  arrays.  Also,  the 
programming  allows  shifting  in  time  and  position  by  TSHIFT,  XSH1FT,  YSHiFT, 
and  ZSH1FT,  all  of  which  are  set  to  zero.  If  needed,  these  could  be  changed. 

There  is  a  separate  subroutine  to  obtain  the  data  for  each  of  the  6 
boundaries.  Each  is  compiled,  and  called  in  HULL,  if  the  corresponding 
boundary  is  declared  type  9  and  the  dimension  is  3  (e.g.,  LB0UND=9  and 
DIMENSIONS)  • 


♦  ♦  RUNSTREAM  FOR  KEEL  ♦  ♦ 


W0RTHAN(STMFZ#T1000,P2>  KEEL  RUN#  PROB  8*16.20  8/15/ BA 

ACCOUNT (♦*♦*♦♦)  WORTMAN  B309  X6028 

COPYSBF(INPUT#OUT°UT) 

REWIND ( INPUT ) 

COMMENT.  SPHERICAL  RESERVOIR  (50  PS  I  SHOCK  IN  STINF).  PROB  8*16.20. 
COMMENT.  START  3D  RUN  FROM  FROM  CYCLE  0  OF  PROB  8*05.09. 

COMMENT.  THIS  IS  PROBLEM  8*06.20  BEING  RERUN  WITH  NEW  HULLUP  AND  BND? 
ATTACH (TAPE**# HULLS *05P09, ID-JOW) 
DI$POSE#3UTPUT»ST-MFAIEI»UN-RJE1*01#*PR. 

ATTACH (PL# HULL# ID«JDW> 

ATTACH (HULL  IB# HULL  IB, ID-JDW) 

LIBRARY(HULLIB) 

REQUEST (TAPE*#*PF) 

REQUEST (TAPE9#*PF) 

UPDATE (P-PL#Q»L-1> 

PLANK. 

POST (COMPILE  > 

HAP(OFF) 

FTN#I-SAIL#B«KEEL#PL-500DO#nPT-2#SL«0#R-0#£L-F) 

RETURN (PL# SAIL) 

KEEL. 

CATALQG(TAPE*#HULL8*16P20#  ID-JOW) 

CATALOG (TAPE9#  STA8*16P20#  ID-JDW) 

♦EOR 

♦LIMIT  10000 
♦COMPILE  HULL 

♦  IDENT  KC  HNG 
♦NOABBREV 

♦DELETE  KEEL. *2*9, KEEL. *251 

C  THE  OLD  HULL  INPUT  FILE  AS  UNIT  **. 

C  THE  DATA  ARE  READ  INTO  ECS  AND  THE  UNIT 

♦INSERT  KEEL. *255 

CQMMON/RDWR/EQFF#  ERR 
♦INSERT  KEEL. *272 

IF(EQFF  ,NE.  0)$TOP«  HULLIN.  NO  DATA  ON  NFILOT  («**>.' 

♦INSERT  KEEL. *281 

IF(EDFF  .NE.  0)ST0P»  HULLIN.  NO  MORE  DATA  ON  NFILOT  (-**).» 
♦INSERT  KEEL. *305 

IF(FQFF  .NE.  OISTOP*  HULLIN.  NO  MORE  DATA  ON  NFILOT  (-**).' 
♦INSERT  KEEL. *327 

IF(EOFF  .NE.  0)  STOP  •  HULLIN.  NO  MORE  DATA  ON  NFILOT  {■**).» 
♦DELETE  KEEL. *3*1 

C  NORMALIZE  USING  VALUES  FROM  ATMOSP. 

♦COMPILE  KEEL 

♦  EOR 
KEEL 

PROB  8*16.20 

LBOUNO-9 

RB0UN0-9 

BBOUND-O 

TBOUND-9 

ABOUNO-O 

FBOUND-9 

ATMOS-5 

DIHEN-3 

EOS-2 

GEOM-2 

HEADER 

80  CM  GRID.  START  FROM  8*05.09.  8*16.20. 

RE  LG  AM -1.* 

RELRHO-O. 00120*12 
RELSIE-2.1037*E9 


RELP0-1.01325E6 

IHAX-9 

JHAX-7 

KHAX-8 

T-0.02 

PTST0P-0.05 

NSTN-33 

MESH 

xo-soo.o 

NX-9  DX-80.0 
NY-7  0Y-80.C 
NZ-8  DZ-80.0 
6ENERATE  AIR 
FIREIN  HULL 

SPHERE  RA0IUS-11A0.0  XCENTER-0.0  YCENTER-0.0  ZCENTER-1200.0 
STATIONS 

XS-990  YS-30  ZS-30 

XS-1250  YS-30  ZS-30 

XS-1500  YS-30  ZS-30 

XS-990  YS-30  ZS-500 

XS-1250  YS-30  ZS-500 

XS-1500  YS-30  ZS-500 

XS-990  YS-30  ZS-590 

XS-1250  YS-30  ZS-590 

XS-1500  YS-30  ZS-590 

XS-990  YS-260  ZS-30 

XS-1250  YS-260  ZS-30 

XS-1500  YS-260  ZS-30 

XS-990  YS-280  ZS-500 

XS-1250  YS-280  ZS-500 

XS-1500  YS-280  ZS-500 

XS-990  YS-280  ZS-590 

XS-1250  YS-280  ZS-590 

XS-1500  YS-280  ZS-590 

XS-990  YS-520  ZS-30 

XS-1250  YS-520  ZS-30 

XS-1500  YS-520  ZS-30 

XS-990  YS-520  ZS-500 

XS-1250  YS-520  ZS-500 

XS-15CO  YS-520  ZS-500 

XS-990  YS-520  ZS-590 

XS-1250  YS-520  ZS-590 

XS-1500  YS-520  ZS-590 

XS-118A  YS-500  ZS-30 

XS-118A  YS-500  ZS-500 

XS-118A  YS-500  ZS-590 

XS-1560  YS-30  ZS-500 

XS-1S60  YS-280  ZS-500 


♦♦  RUNSTREAM  FOR  HULL  ♦♦ 


W0RTMAN(STNFZ,T50»P1)  HULL#  PROB  8416.20  8/15/84 

ACCOUNT (♦*♦♦♦♦)  8309  X602B 

COPYSBF( INPUT# OUTPUT) 

REWIND! INPUT) 

COMMENT.  RERUN  OF  8406.20  TO  TEST  NEW  B0UND9  COOING. 

COMMENT.  SPHERICAL  RESERVOIR  (SO  PSI  STINF).  PROS  8416.20. 

COMMENT.  3D  RUN  STARTED  FROM  PROB  6405.09  THRU  F IREIN/HULLIN. 
COMMENT.  4  BOUNDARIES  SET  TO  TYPE  9.  INPUT  FROM  HULLUP. 

ATTACH (TAPE  31# P8 416 20LB# I DbJDW) 

ATTACH! TAPE32# P841620R8# I D-JDW) 

ATTACH (TAPE  34# PB41620TB# I D*JDW) 

ATTACH (TAPE  36# PB41620FB# 10- JDW) 

COMMENT.  UPDATE  CHANGES  (CORRECTIONS)  FOR  HULL  FROM  FILE  HULCORR 
BEGIN (GETMF A# FILE# LFbHULCORR#PFbHULCORR#UN>JDW) 

COMMENT.  UPDATE  CHANGES  FROM  FILE  NEWBND9  IN  MFA  TO  BND9  IN  MF2 

BEGIN (GETMF A#FILE#LF-BND9#PFbNEWBND9«UN> JDW) 
DISPQS£#0UTPUT»ST-MFAIEI#UN-RJE1401#PPR. 

ATTACH(T4»  HULLB416P20#  ID- JDW) 

ATTACH(T9,  STA8416P20#  ID-JDW) 

RE0UEST(TAPE4»*PF) 

COPYP(T4#TAPE4#iOOO) 

RETURN (T4) 

REWIND (TAPE4) 

RE0UEST(TAPE9#*PF) 

COPYP ( T9#  T  APE9# 1000 ) 

RETURN(T9) 

REWIN0(TAPE9) 

ATTACH) HULLIB# HU LLIB# ID-JDW) 

LIBRARY! HULL IB) 

ATTACH(PL#HULL#ID-JOW) 

UPDATE ( P»PL#0#  L-l) 

SWITCH (6# ON) 

PLANK. 

SWITCH (6#  OF  F ) 

POST(COMPILE) 

MAP! PART) 

FTN(  A# I*SAI L»  B«HULL  »QPT "2#  SL-0» R-0# E  L-F# LCM-I ) 

RETURN! PL#  SAIL# COMPILE ) 

HULL. 

EXIT(U) 

CATAL0G(TAPE4#HULLB416P20#  ID-JDW) 

CATAL0G(TAPE9#  STAB416P20#  ID-JDW) 

♦EOR 

♦LIMIT  10000 
♦READ  HULCORR 
♦READ  BN09 
♦COMPILE  HULL 
♦EOR 
HULL 

PROB  8416.20 

CYCLE-0 

INPUT 

ST  ABF-0.05 

MRELER-l.OE-6 

CSTOP-200 

END 


*•  HULL  CORRECTIONS.  FILE  HULCORR  ** 


•I DENT  BNDCORR 

*t  CORRECT  AN  ERRORS  CHAN6E  L BOUND  TO  ABOUND 

♦DELETE  HULL. 11059 
•KEEPTO  *1  ABOUND*  AND  LAHB 
♦/  AMEND  SHKMV  C  ALLS  * 

*/  HAKE  BOUNDARIES  MORE  GENERAL .  NO  Y-0.0  OR  X-0.0  FOR  30. 

•  /  FIRST  PHASE#  TIME  -T  FOR  EXTERNAL  CELL*  -T  ♦  0.5*DT  ON  BOUNDARY. 

•  /  SECOND  PHASE#  TIHE  ■  T  ♦  OT. 

♦DELETE  HULL. 11062 

CALL  SHKWV (T*0.5*DT»DISX» YO# ALT# UAB# VAB.HAB# EAB#RHOAB ) 

•DELETE  HULL. 11389 

CALL  SHKWV(T*0.5*DT#  DISX#  DISY* ALT#  UBB# VBB*  WBB*  EBB#RHOBB) 

♦DELETE  HULL. 1163* 

CALL  SHKWV(T*0.5*DT#X0#DISY#ALT* ULB* VLB# VLB* EL B* RHOLB) 

♦DELETE  HULL. 12*** 

CALL  SHKWV(T+DT*DISX*Y0#ALT#UA8#VAB#WAB*EAB#RH0AB) 

♦DELETE  HULL. 12600 

CALL  SHKVV{T*DT»DISX*DISY#ALT*UB8#VBB#VBB*EBB#RH0BB) 

•DELETE  HULL. 12815 

CALL  SHKWV(T«-DT»XO#DISY#ALT»ULB»VLB*WLB»ELB»  RHOLB) 

♦DELETE  HULL. 13209 

CALL  SHKWV2(T+DT*DISX*DISY*UBB#VBB*EBB*Rr40BB) 

♦DELETE  HULL. 13329 

CALL  SHKWV2( T*DT*DISX#DISY*ULB*VLB* E LB* RHOLB  ) 

♦IDENT  DTFR2D 

♦/  BETTER  2D  TINE  STEP.  FIXED  GAHHA  ONLY. 

♦INSERT  HULL. *07* 

C  COMPUTE  CS  FOR  CELL(I*J) 

C  THIS  IS  FOR  FIXED  GAMMA.  OTHERWISE*  USE  EO.  OF  STATE. 

EFORCS  •  H(N1*A)-0.5*(H(N1*2)**2  ♦  H(N1*3)**2) 

CS  -  SORT (GAMMA* (6 AMM A-l . 0 )*E FORCS ) 

C  END  CS  INSERT. 

♦IDENT  HULLCOR 
•NOABBREV 

*/  ERRORS  IN  SHORE  CELL  OPTION  (CORRECTION  THANKS  TO  PAUL  LEWIS) 

*/  POSSIBILITY  OF  ALL  FLUIO  CELLS  OMITTED. 

♦INSERT  HULL. *116 

C  CHECK  FOR  RIGHT  CELL  ALL  FLUIO 

IF(LSR  .EO.  0)  GOTO  105 
♦DELETE  HULL. *127 

105  RHOP  •  H(NR1*5) / VOLR 
♦INSERT  HULL. *1*5 

C  CHECK  FOR  CELL  ABOVE  ALL  FLUID 

IF(LSA  .EO.  0)  GOTO  115 
♦DELETE  HULL. *157 

115  RHOA  ■  H(N Al+5) / VOLA 
♦DELETE  HULL. 11256 

IF(LSI  .LE.  OJGOTO  6 
VOL  -  VOL* 0.5 
•INSERT  HULL. 11260 
6  CONTINUE 
♦INSERT  HULL. 11517 

IF (LSI  .LE.  0)G0T0  5 
♦INSERT  HULL. 11521 
5  CONTINUE 

*/  END  OF  THE  RAUL  LEWIS  CORRECTION. 

♦/  TBR1  HAS  DUPLICATE  ADDRESS  IF  NH.GT.l 

♦DELETE  HULL. 12336*HULL. 12337 
DO  17  IM-1*NN 

17  H(NA*NHN*IH>  ■  H ( NT BR  +  NHH* IM) *  DEN 
♦/  REMOVE  AN  ANOMALY  BELOW  HULL. 13661 

♦DELETE  HULL. 13661.HULL. 13662 
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**  BOUND?  COOING.  FILE  NEHBND9  •* 


♦IDENT  BOUNO? 

*/  X80UND-9  IS  INPUT  FPON  A  PREVIOUS  HULL  RUN  AT  ABOUND. 

*/  MAXIMUM  INPUT  PLANE  A  76X76  MESH  OF  POINTS.  7/24/84. 

♦INSERT  HULL. 704 

♦KEEPTO  ENOBND?  LB0UND9  OR  R BOUND?  OR  BB0UND9  OR  TBOUN09  OR  ♦ 
ABOUND?  OR  F BOUND? 

♦INCLUDE  DEEP-SIX  REZONE 
♦LABEL  ENDBND9 
*/  FILES  FOR  XBOUND*?. 

♦INSERT  HULL. 1146 
♦KEEPTO  * 1  LB0UND9 

♦  TAPE31#  TAPE61* 

♦KEEPTO  *1  P.  BOUND? 

♦  TAPE32»TAPE62» 

♦KEEPTO  +1  BB0UND9 

♦  TAPE33.  TAP£63» 

♦KEEPTO  *1  T BOUND? 

♦  T A  PE34»  TAP  E64* 

♦KEEPTO  *1  ABOUND? 

♦  TAPE35#TAPE65> 

♦KEEPTO  ♦ 1  F BOUND? 

♦  TAPE36.TAPE66* 

♦INSERT  HULL. 8161 

♦KEEPTO  *1  LB0UND9  OP  RB0UND9  OR  BB0UND9  OR  T BOUND?  OR  ♦ 

ABOUND?  OR  F BOUND? 

COMMON/SHIFT/TSHIFTj  XSHIF T» YSHI FT» Z SHIFT 
♦INSERT  HULL. 8369 

♦KEEPTO  ♦A  LBOUND?  OR  R BOUND?  OR  BBOUND?  OR  TB0UND9  OR  ♦ 

ABOUND?  OR  F BOUND? 

TSHIFT  -  0.0 
XSHIFT  -  0.0 
YSHIFT  •  0.0 
ZSHIFT  •  0.0 
♦INSERT  HULL. 9250 

♦KEEPTO  ENDBND9  0IMEN3 

♦KEEPTO  ENDBND9  LBOUND?  OR  RBOUND?  DR  BBOUND?  OR  TBOUND?  OR  ♦ 

ABOUND?  OR  F BOUND? 

C  ADD  HB  ARRAYS  TO  LCN.  USE  NEXT  LARCER  NULTIPLE  OF  8. 

C  MM  4000  AS  OF  8/21/84.  Ill  ASSUMES  MAXI  IN.  JN.KNI  .LE.  200. 
CALL  FLECS ( IB9ECS  ) 

IF( ABOUND  .EQ.  9)  IB9ECS  -  IB?ECS  ♦  4000 

IFIFBOUND  .EQ.  91  IB9ECS  -  IB9ECS  ♦  4000 

IF C BBOUND  .EQ.  9)  IB9ECS  ■  IB9ECS  ♦  4000 

IF  I TBOUND  .EQ.  9)  IB9ECS  -  IB9ECS  ♦  4000 

IF  I LBOUND  .EQ.  9}  IB9ECS  -  IB9ECS  ♦  4000 

IFIRBOUND  .EQ.  9)  IB9ECS  -  IB9ECS  ♦  4000 

CALL  SETECSI IB9ECS) 

♦LABEL  ENOBND? 

♦INSERT  HULL. 11058 
•KEEPTO  *1  ABOUND? 

CALL  ABND9D3 (T*0.5*DT»DI$X> YO/ ALT# U AB#  V AB»  WAB*  E AB.RHOAB ) 
♦SKIPTO  *4  ABOUND? 

•INSERT  HULL. 11092 
♦SKIPTO  STBND9  FB0UND9 
♦INSERT  HULL. 11102 
♦LABEL  STBND9 

•KEEPTO  ENDBN09  FB0UND9 

DISX  •  XIII  -  0.5*0X11) 

DISY  •  YIJMAX)  -  0.5*DYI JMAX) 

ALT  -  ZIK)  -  0,5*0Z(K) 

CALL  FBND9D3 (T»DISXfDISY»ALT#HINl*l)»H(Nl*2)*HCNl*3)»HINl*4)» 

♦  RHOFB) 
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H(N1  +  5 )  ■  RHOFB*DX< I»*DY< JNAX)*DZ(K> 

♦LABEL  ENDBND9 
♦INSERT  HULL.113B8 
♦KEEPTO  *1  BB0UND9 

CALL  BBND9D3  (T+0.5*DT»DISX»DISY*ALT#UBB#VBB#WBB»EBB#RrtQBB) 
♦SKIPTO  *1  BBOUND9 

♦INSERT  HULL. 11630 
♦KEEPTO  *1  L80UND9 

CALL  LBND9D3 (  T+0. 5*DT»X0>DISY>  ALT.ULB.VLB. ML  B>  ELB.RHOLB ) 

♦SKIPTO  *4  LB0UND9 

♦INSERT  HULL. 11665 
♦SKIPTO  STBND9  RB0UND9 
♦INSERT  HULL. 1167? 

♦LABEL  STBND9 

♦KEEPTO  EN0BN09  RB0UND9 

DISX  ■  X(IMAX)  -  0. 5*DX ( IHAX) 

OISY  •  Y  ( J  )  -  0. 5*PY ( J ) 

ALT  «  Z(K)  -  0.5*DZ(K) 

CALL  RBND9D3CT. DISX. OISY# ALT.H (NL*1  >.H< NL»2 > .HC NLO). HCNL  +  4)  . 

♦  RHORB) 

HlNL+5)  ■  RH0R8*DX(IHAX)*DYm*DZ<K) 

♦LABEL  EN0BND9 
♦INSERT  HULL. 12362 
♦SKIPTO  STBND9  TB0UND9 
♦INSERT  HULL. 12370 
♦LABEL  STBND9 
♦KEEPTO  EN0BN09  TB0UN09 

DISX  -  X ( I )  -  0 • 5  *DX ( I ) 

DISY  -  Y(J)  -  0. 5*0Y  < J  ) 

ALT  -  Z(KNAX)  -  0.5*DZ(KHAX) 

CALL  TBND903 t T>DISX>DISY.ALT.H(NA+l).H(NA+2).H(NA+3)fH(NA+4). 

♦  RHOTB ) 

H(NA+5I  -  RHOTB*OX( 1 1 ♦ 0 Y < J ) ♦ D Z ( KNAX ) 

♦LABEL  EN0BN09 
♦INSERT  HULL. 12443 

♦KEEPTO  *1  AB0UN09 

CALL  A8ND9D3 ( T*DT,DI $X# YO. ALT# UAB. V  AB.  W  AB. EAB. RHOAB > 

♦SKIPTO  41  AB0UND9 

♦INSERT  HULL. 12455 
♦SKIPTO  STBND9  FB0UND9 
♦INSERT  HULL. 12467 
♦LABEL  STBND9 
♦KEEPTO  EN0BN09  FB0UN09 

DISX  -  X(I)  -  0. 5*DX < I ) 

OISY  -  Y(JNAX)  -  0. 5*DY( JMAX) 

ALT  -  Z<K)  -  0. 5*DZ IK) 

CALL  FBND9D3(T*DT.0ISX.DISY.ALT»H(N1*1).H(N1*Z>»H<N1«-3)»H(N1*4> 

♦  RHOFB J 

HtNl*4>  ■  H(N1*4 )  ♦  0,5* IH(N1+1>**2  ♦  H(N1*2>**2  ♦  H(N1»3J**2) 
H(N1*  5 )  «  RHOFB* DX< I )*DY( JHAX ) ^DZ (K  ) 

♦LABEL  ENDBND9 
♦INSERT  HULL. 12599 

♦KEEPTO  *1  8B0UND9 

CALL  BBND9D3<  T+DT.DISX.DISY.ALT.UBB. VBB.WBB. EBB.RHOBB) 

♦SKIPTO  *1  BB0UN09 

♦INSERT  HULL. 12671 

♦KEEPTO  ENDBN09  FB0UND9 

IF ( I  .EO.  IH AX ) GOTO  5 
DTH  -  OTH  -  FEJ(JL) 

DNTH  -  DHTH  -  FNJ(JL) 

5  CONTINUE 
♦LABEL  EN08N09 
♦INSERT  HULL. 12672 

♦SKIPTO  STBND9  FB0UND9 


♦INSERT  HULL. 12690 
♦LABEL  ST  BN 09 

♦INSERT  HULL. 12814 
♦KEEPTO  *1  LB0UN09 

CALL  LBND9D31 T+DT»X0*DISY»ALT* ULB*VLB»WLB>ELB#RHOLB) 

♦SKIPTO  *1  LB0UN09 

♦INSERT  HULL. 12826 

♦SKIPTO  STBND9  RB0UND9 
♦INSERT  HULL. 12838 
♦LABEL  STBN09 

♦KEEPTO  EN0BND9  R BOUND 9 

DISK  «  X(INAX)  -  0. 5*DX ( INAX) 

DISY  «  Y ( J )  -  0. 5^DY  <  J ) 

ALT  -  ZCKJ  -  0.  5*DZ  ( K ) 

CALL  R8ND9D3(T+DT»DISX/DISY,»ALT/H(NL*l}/H<NL*2)»H<NL  +  3)*H(NL*-4) 

♦  RHORB) 

H(NL  +  4 1  «  H( NL*4 )  ♦  0.  5*  <  H  (NL* 1  )**2  ♦  H(NL*2>**2  ♦  H<NL+3)**2) 
HINL  +  5 }  -  RHORB*OX(INAX)*DY(J)*DZCK) 

♦LABEL  ENDBN09 
♦INSERT  HULL. 12890 
♦KEEPTO  EN0BN09  RBOUN09 
DTH  «  DTH  -  FEL 
DHTH  -  DNTH  -  FML 
♦LABEL  EN0BND9 
♦INSERT  HULL. 12898 
♦SKIPTO  STBND9  RB0UND9 
♦INSERT  HULL. 12911 
♦LABEL  STBN09 
♦INSERT  HULL. 12980 
♦KEEPTO  EN0BND9  T BOUND 9 

IF ( I  .EG.  INAX  .OR.  J  .EG.  JNAX)GOTO  5 
DTH  ■  DTH  -  FEK(KL) 

DNTH  -  DNTH  -  FNK(KL) 

5  CONTINUE 
♦LABEL  ENDBND9 

♦INSERT  HULL. 12981 
♦SKIPTO  ST8ND9  TB0UND9 

♦INSERT  HULL. 12998 
♦LABEL  STBND9 

♦INSERT  HULL. 13022 
♦SKIPTO  STBND9  TB0UND9 

♦INSERT  HULL. 13030 

♦LABEL  STBND9 

♦KEEPTO  ENDBND9  TB0UND9 

disx  -  xm  -  o.s*dx(I) 

DISY  •  Y(j)  -  0.54DYU) 

ALT  -  Z(KNAX)  -  0.5*DZ(KNAX) 

CALL  TBND9D3(T«-DT»DISX,DISY»ALT»H<NA+l)»H(NA*2),H(NAO)»H<NAt4) 

♦  RHOTB) 

H(NA*4)  -  H( N A+4 )  ♦  0. 5*( HI NA* 1)**2  ♦  HCNA*2>4*2  ♦  H(NA»3>**2> 
H(NA«-5)  «  RHOTB* DX(  Z)*DY(J)*DZ(KNAX) 

♦LABEL  ENDBN09 

♦INSERT  HULL. 11010 

♦KEEPTO  ENDBND9  LB0UND9  AND  D1NEN3 

SUBROUTINE  LBND9D3( TH#  XH»  YH»ZH#  U»V»WjE»RHO) 

C  BOUNOARY  SUBROUTINE  FOR  LBOUND-9.  INPUT  PREPARED  BY 

C  HULLUP  IS  ASSUNED  ON  FILE  TAPE31  FOR  X  •  IXIN). 

C  XH  SHOULD  BE  XO.  XH  ♦  XSHIFT  SHOULD  BE  XIN. 

C 

C  KN  ♦  2  RECORDS  PER  TINE  DUNPi 

C  HEADER--555.0*0LD  PROB* TI 1» TI 2» X IN. NNH» JN.KN 

C  GRID — (YIN(J )»J*1»JN). (ZIN(K)»K«1»KN) 

C  HYDRO  DATA  NNH/POINT  FOR  TINES  Til  AND  TI2 

C  KN  ROWS  AT  (YIN<J);ZIN<K)*J*1»JN)  FOR  K»1,KN. 


o  o  o  o 


COMM ON /SHI  FT /7SHIFT* XSHIFT, YSHIFT, ZSH1FT 
VALUES  IN  SHIFT  SET  IN  HULL IN. 

FOLLOWING  ARRAYS  SET  TO  ACCOMODATE  JN  .LE.  200. 

C0MM0N/H81/  TO  ACCOMODATE  LEVEL  2  STATEMENT. 
C0MM0N/HB1/H8K  A 000) 

DIMENSION  H( 5 ) 

DIMENSION  YIN!200),ZIN{200) 

LEVEL  2,  HB1 

NAMEL IST/HEADIN/HEAD1,  OLPRDB/T I1»TI2,XIN,NNH,JN,KN 
NAME LIST/GRID/XH»XSHIFT»X,XIN, YH,  YS  HIF  T,  Y,  YI N,  ZH,  ZSHIFT»Z,ZIN 
SET  FILE  NUMBER.  STATE  ARRAY  DIMENSIONS. 

DATA  IBND, NBND, JMX, KM  X*NMX/31> 61*200*200,4000/ 

NMX  IS  ADDED  TO  LCM  IN  HULLIN  (HULL925G). 

SET  TF1  TO  -1.0  FOR  INITIAL  VALUES. 

DATA  TF1/-1.0/ 

SHIFT  OATA  FROM  HULL  COORDINATES  TO  BOUNDARY  COORDINATES. 
T  ■  TH  +  TSHIFT 
X  ■  XH  ♦  XSHIFT 
Y  «  YH  ♦  YSHIFT 
Z  ■  ZH  ♦  ZSHIFT 

IS  THIS  AFTER  THE  INITIAL  ENTRY? 

IFITF1  .GT.  -0. 9) GOTO  30 
TF1  -  1.0 

RE  AO  FIRST  HEADING*  6RID*  AND  HYDRO  SET.  CHECK. 
REWIND  IBND 

READ ( IBND)HE ADI*  OLPROB, TI 1*TI2*  XIN*  NNH* JN*  KN 

IF(EOF(IBNO)  .NE.  0  >GDTO  100 

TTEST  •  0.00001* ( TI 2  -  Til I 

XTEST  •  0. 0000001 *X IN 

WRITE ( 6»  5) 

5  FORMAT!//'  INITIATE  LB0UND9  INPUT.  FIRST  HEADER  RECORD") 
WRITE ( b*  HE  AO  IN) 

IF! A85 (HEAD1  -  555.0)  .GT.  O.ODGDTO  110 
NNHYPR  «  NNHPJN 
NNHYPP  •  NNHYPR*2 

IF( JN  .GT.  JMX  .OR.  KN  .ST.  KMX  .OR.  NNHYPP  .6T.  NMX)60T0  120 
IF (T  +  TTEST  ,LT.  TIDGOTO  130 
0LPRB1  ■  OLPROB 

READ  GRID  FOR  INPUT  BOUNDARY  PLANE. 

READ! IBND)  < YIN ! J )* J -1* JN  ) » t ZIN IK ) » K-l, KN ) 

YTEST  -  I  YIN ! 2)  -  YINI1)  )*0.00001 
ZTEST  -  I  ZIN! 2)  -  ZINIl)  )*0. 00001 
10  IFIT*TTEST  .LE.  TI2)G0T0  20 

FIND  THE  NEXT  HEADER  RECORD. 

14  READ  I IBNDIHE AD1*QLPR0B*TI1*TI2*XIN*NNH* JN*KN 
IFIEOFIIBND)  .NE.  0)G0T0  160 

IFIOLPROB  .NE.  0LPRB1  .OR.  HE ADI  .NE.  555.0)G0T0  14 

THIS  IS  A  NEW  HEADER  RECORD.  BYPASS  THE  GRID  RECORD. 
READ! IBNOIOUMMY 
GOTO  10 

READ  HYDRO  RECORDS.  STORE  IN  FILE  NBND, 

20  REWIND  NBND 
DO  22  K-1,KN 

READ! IBND)  I HB1 t L )» L-l* NNHYPP ) 

WRITE INBND)  (HB1IL)*L»1*NNHYPP> 

22  CONTINUE 
24  REWIND  NBNO 

READ  FIRST  2  ROWS  FROM  FILE  NBND. 

READ  I NBND )  I HB1 ( L )* L-l* NNHYPP ) 

RE AO INBND)  IH81 INNHYPP+L I #L«1, NNHYPP) 


r*x»] 


JIN  ■  2 
KIN  *  2 
ISW  -  2 
IN0H1  •  0 
INDH2  •  NNH/PP 
TDENOM  •  1.0MTI2  -  TUI 
30  IF( T-TTEST  .ST.  TI2JG0T0 

IF < ASS ( X  -  XIN)  .GT.  XTESTJGOTO  140 

IF ( Y IN { 1 )  -  Y  .GT.  YTEST  .OR.  Y  -  YINTJN)  .GT.  YTEST )GOTO  140 
IF ( Z IN < 1 )  -  Z  .GT.  ZTEST  .OR.  Z  -  ZINCKN)  .GT.  ZTESTJGOTO  140 
32  IF( Z-ZTEST  .LT.  ZIN I  KIN)  J60T0  38 
IF(ISW  .EQ.  1JG0T0  34 
REAO(NBND)  < HB1 ( L >, L - 1, NNH YPB ) 

ISW  -  1 

IN0H1  •  NNHYPP 
INDH2  «  0 
GOTO  36 

34  READ ( NBND )  ( MSI ( NNHYPP+l ) , L«l> NNHYP P > 

ISW  •  2 
I NO HI  >  0 
INDH2  -  NNHYPP 
36  KIN  •  KIN*  1 

ZDENQM  ■  1.0/  (ZIN(KIN)  -  ZIN(KIN-U) 

JIN  ■  2 
GOTO  32 

38  IF ( Z+ZTEST  .LT.  ZIN(KIN-l)  IGOTD  24 
C  LOCATE  THE  CORRECT  COLUMN. 

IF C Y* YTEST  .LT.  YINUIN-1J  >  JIN  -  2 
40  IF( Y-YTEST  .LT.  YIN(JIN)  JGOTO  42 
JIN  ■  J  IN+l 
GOTO  40 

C  NOW.  SET  ADDRESSES  AND  DO  THE  INTERPOLATION. 

42  IREF  •  (JIN  -  2)*NNH 

YRATIO  «  (Y  -  YIN(J.N-l)  )/(YIN(JIN)  -  YIN(JIN-l)  ) 

TRATIO  «  (T  -  TI1J4T0EN0H 
ZRATIO  ■  (Z  -  ZIN(KIN-I)  )*ZDENON 
C  THREE  INDICES  THAT  ARE  1  OR  2  BELOW  REPRESENT  LOWER  DR 

C  HIGHER  VALUES  OF  Y,  l,  OR  T»  RESPECTIVELY. 

LF111  -  INCH  1  ♦  IREF 
LF211  -  LF 111  4  NNH 
LF112  ■  LF 111  ♦  NNHYPR 
LF212  -  LF112  4  NNH 
LF121  ■  INDH2  4  IREF 
LF221  »  LF121  4  NNH 
LF122  -  LF 12 1  4  NNHYPR 
LF222  •  LF 122  4  NNH 
DO  50  L*l,  5 
Fill  •  HBKLF1U  4  L) 

FY11  «  Fill  4  (HBKLF2114LJ  -  F111)*YRATI0 
F121  •  HB1 ( L  F121+  L ) 

FY21  ■  F121  4  (H8KLF2214L)  -  F121)*YRATI0 
FYZ1  •  FY11  4  (FY21  -  F Yll )*ZR ATIO 
F112  •  HBKLF112  4  L) 

FY12  «  F112  ♦  (HB1UF2124L)  -  F112)*YRATI0 
F122  •  HBKLF1224L) 

FY22  -  F122  4  ( HB 1( LF 22 2* L )  -  F122)*YRATI0 
FYZ2  •  FY12  4  (FY22  -  FY12)*ZRATI0 
H(LJ  -  FYZ1  4  ( F YZ2  -  FYZ1 J4TRATI0 
50  CONTINUE 
U  •  H  ( 1 ) 

V  ■  H  ( 2  J 
W  »  H  (  3 ) 

£  •  H  ( 4 1 
RHO  ■  H(  5) 
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I NPOSED-SOLUT I ON  BOUNDARIES  FOR  THREE-DINENSIONAL  HULL 
<U>  ARHV  BALLISTIC  RESEARCH  LAB  ABERDEEN  PROVING  GROUND 
HD  J  D  HORTHAN  AAR  86  BRL-HR-3499 


RETURN 

C 

100  WRITE (6*109) 

STOP*  ABORT/  LBND9D3.  NO  DATA  ON  FILE.* 

109  FORMAT!//'  ABORT/  LBND9D3.  NO  DATA  ON  FILE.'t 

110  WRITE!6/115) 

STOP*  ABORT  LBND9D3.  FIRST  FILE  WORD  NOT  959. 0* 

119  FORMAT!//*  ABORT/  LBN09D3.  HE  ADI  NOT  995.0*) 

120  WRITE (6/ 125 )JN/ JNX/KN/KMX/NNHYPP/NMX 

STOP*  ABORT*  LBND9D3.  DIMENSION  TOO  SMALL.* 

125  FORMAT!//*  ABORT/  LBND9D3 .  DIMENSION  TOO  SMALL.*/ 

♦  •  JN/JMX/KN/KMX,NNHYPP/NMX»*»6I8) 

130  WRITE(6/139)TH/TSHIFT/T/TI1 

STOP*  ABORT/ LBND9D3.  REOUESTED  T  LT  INITIAL  TIME  IN  FILE.* 

139  FORMAT!*  ABORT*  LBN09D3.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.*/ 

♦  *  TH/TSHIFT/T/TIi-  '/1P4E15,?) 

140  WRITE (6*149) 

WRITE (6/GRID) 

STOP*  ABORT/  LBND9D3 .  SPATIAL  VARIABLE  OUTSIDE  RANGE. • 

145  FORMAT!//*  ABORT*  LBN09D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.*) 

160  WRITE! 6/1 65 )TH/TSHIFT*T*TI1/TI2 

STOP*  ABORT/  LBN09D3.  END  OF  FILE  BEFORE  T.* 

165  FORMAT!*  ABORT/  LBND9D3.  END  OF  FILE  BEFORE  T.*/ 

♦  '  TH/TSHIFT.T/TI1/TI2-  *#1P5E15.7) 

END 

* LAB EL  ENDBN09 

•KEEPTO  ENDBN09  RB0UND9  AND  DIMEN3 

SUBROUTINE  RBND9D3! TH/ XH, YH/ ZH/U/ V/ W /£/ RHO ) 

C  BOUNDARY  SUBROUTINE  FOR  RBQUND-9.  INPUT  PREPAREO  BY 

C  HULLUP  IS  ASSUMED  ON  FILE  TAPE32  FOR  X  >  (XIN), 

C  XH  SHOULD  BE  X (INAX )-0.3*DX! IMAX ) .  XH  ♦  XSHIFT  SHOULD  BE  XIN. 

C 

C  KN  ♦  2  RECORDS  PER  TIME  DUMP* 

C  HEADER— 555.0/OLD  PR0B/TU/TI2/  X IN* NNHz  JN/KN 

C  GRID—  IYIN!J)/J«l#JN)#!2INIlO/K-l/KN) 

C  HYDRO  DATA  NNH/POINT  FOR  TIMES  Til  AND  TI2 

C  KN  ROWS  AT  ( YIN! J )/ ZIN(K)/ J«i/ JN)  FOR  K-1»KN. 

C 

COMMON/SHIFT/TSHI FT* XSHIFT. YSHIFT/ZSHIFT 
C  VALUES  IN  SHIFT  SET  IN  HULLIN. 

C 

C  FOLLOWING  ARRAYS  SET  TO  ACCOMODATE  JN  .LE.  200. 

C  C0MN0N/HB2/  TO  ACCOMODATE  LEVEL  2  STATEMENT. 

CONHON/HB2/HB2! 4000) 

DIMENSION  H( 5 ) 

DIMENSION  YIN!200)/ZIN(200) 

C 

LEVEL  2/  HB2 

NAMELIST/HE  ADIN/HEAD1* OLPROB*  T II, TI 2/XIN/NNH/ JN,KN 
NAME L IST/GRIO /XH, XSHIFT* X* XIN/ YH/YSHIFT/Y/ YIN* ZH/ZSHI FT* Z/ZIN 
C  SET  FILE  NUMBER.  STATE  ARRAY  DIMENSIONS. 

OATA  IBN0,N8N0,JMX,KHX/NMX/32* 62,200/200/4000/ 

C  NMX  IS  ADDED  TO  LCM  IN  HULLIN  (HULL9250). 

C  SET  TF1  TO  -1.0  FOR  INITIAL  VALUES. 

DATA  TF1/-1.0/ 

C 

C  SHIFT  DATA  FROM  HULL  COORDINATES  TO  BOUNDARY  COORDINATES. 

T  •  TH  ♦  TSHIFT 
X  •  XH  ♦  XSHIFT 
Y  •  YH  ♦  YSHIFT 
Z  •  ZH  ♦  Z SHI  FT 

C  IS  THIS  AFTER  THE  INITIAL  ENTRY T 

IF (TF 1  .CT.  -0.9»GOTO  30 
TF1  •  1.0 


C  READ  FIRST  HEADING#  GRID#  AND  HYDRO  SET.  CHECK. 

REWIND  IBND 

RE  ADC IBND) HE  ADI# QLPROB# TI 1# TI2# X I N» NNH# JN# KN 
IF  I E  3  F I  IBND)  .NE.  0  JGOTO  100 
TTEST  •  0. 00001* ( TI2  -  Til) 

XTEST  -  0. 0000001 *X IN 
WRITE ( 6#  5 ) 

5  FORMAT!//*  INITIATE  R80UND9  INPUT.  FIRST  HEADER  RECORD* • ) 
WRITE ( 6#HE ADIN) 

IF( ABSIHEAD1  -  555.0)  .GT.  0.01IG0T0  110 
NNHYPR  •  NNH* JN 
NNHYPP  -  NNHYPR*2 

IFIJN  .GT.  J NX  .OR.  KN  .GT.  KNX  .OR.  NNHYPP  .GT.  NHX)60T0  120 
IF(T*TT£5T  .LT.  Til ) GOTO  130 
0LPRB1  ■  OLPROB 

C  READ  GRID  FOR  INPUT  BOUNOARY  PLANE. 

READ  (  IBND )  < YIN< J )» J-l# JN)#  (ZIN(K)#K*1#KN) 

YTEST  •  ( YIN(  2)  -  YINU)  )*0.00001 
2TEST  -  ( Z IN  1 23  -  ZIN(l)  )*0.00001 
10  IF( T+TTE ST  .LE.  TI2)G0T0  20 
C  FIND  THE  NEXT  HEADER  RECORD. 

14  READ (IBND) HE  ADI, 0LPR08# TI 1, TI2 # XIN. NNH# JN#KN 
IF( EOF ( IBND)  .NE.  0)GQT0  160 

IFIOLPROB  .NE.  0LPR61  .OR.  HEA01  .NE.  555.0)G0T0  14 
C  THIS  IS  A  NEW  HEADER  RECORD.  BYPASS  THE  GRIO  RECORD. 

REAO( IBND) DUHNY 
GOTO  10 

C  READ  HYDRO  RECORDS.  STORE  IN  FILE  NBND. 

20  REWIND  NBND 
DO  22  K-1#KN 

READ) IBND)  ( HB2 ( L )• L-l# NNHYPP ) 

WRITE ( NBND )  (HB2( L) » L-l# NNHYPP) 

22  CONTINUE 
24  REWIND  NBND 

C  READ  FIRST  2  ROWS  FROM  FILE  NBND. 

REAO(NBND)  ( H82 ( L )# L-l# NNHYPP ) 

READ  I NBND)  !  HB2  C  NNHYPP *L ) # L-l#  NNHY*P ) 

JIN  -  2 
KIN  -  2 
ISW  -  2 
INDH1  •  0 
INDH2  •  NNHYP» 

TDENOM  -  1.0/ITI2  -  Til) 

30  IF(T-TTE5T  .GT.  TI2)G0T0  14 

IF ( ABS ( X  -  XIN)  .GT.  XTEST)GOTO  140 

IF<YIN<1)  -  Y  .GT.  YTEST  .OR.  Y  -  YINUN)  .GT.  YTE$T)GOT 0  140 
IF(ZIN(1)  -  Z  .GT.  ZTEST  .OR.  Z  -  ZIN(KN)  .GT.  ZTEST)GOTO  140 
32  IF ( Z-ZTEST  .LT,  ZIN(KIN)  )GOTO  38 
IF ( I S W  .£0.  1IG0T0  34 
READ! NBND)  (HB2(l )#l-l, NNHYPP) 

ISW  -  1 

INDH1  -  NNHYPP 
INDH2  -  0 
GOTO  36 

34  READ ( NBND)  < HB2 ( NNHYPP*L )# L-l# NNHYPP ) 

ISW  •  2 
INDH1  -  0 
IN0H2  -  NNHYPP 
36  KIN  -  KIN* 1 

Z DEMON  -  1.0/ (ZIN(KIN)  -  ZIN(KIN-l)) 

JIN  -  2 
GOTO  32 

38  IF(Z*ZTEST  .LT.  ZIN(KIN-l)  )GOTO  24 
C  LOCATE  THE  CORRECT  COLUMN. 
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IF(Y+YT£ST  .LT.  YZN(JIN-l)  )  J  IN  •  2 
40  IF( Y-YTEST  .LT.  YIN(JIN)  )G0T0  42 
JIN  -  JIN+l 
GOTO  40 

C  NOW*  SET  ADDRESSES  AND  DO  THE  INTERPOLATION. 

42  IREF  -  (JIN  -  2 ) *NNH 

YRATIO  •  (Y  -  YIN(JIN-I)  ) / (YIN( J IN  >  -  YIN(JIN-l)  I 
TRATIO  «  (T  -  TI1)*TDE'(0H 
ZRATIO  -  (Z  -  Z IN (KIN— 1 }  )*ZDENOM 
C  THREE  INDICES  THAT  ARE  1  OR  2  BELOW  REPRESENT  LOWER  OP 

C  HIGHER  VALUES  OF  Y,  Z*  OR  T*  RESPECTIVELY. 

LF111  «  INDH1  ♦  IREF 
LF21I  •  IF111  ♦  NNH 
LF112  «  LF 111  ♦  NNHYPR 
LF212  •  LF 112  ♦  NNH 
LF121  ■  INDH2  ♦  IREF 
LF221  •  LF121  ♦  NNH 
LF122  «  LF 121  ♦  NNHYPR 
LF222  •  LF 122  ♦  NNH 
DO  50  L-1,5 
Fill  •  HB2 ( LF 111  ♦  L) 

FY11  -  Fill  ♦  (HB2(LF211+L)  -  F111)*YRATI0 
F121  -  HB2 ( L  F121+L) 

FY21  •  F121  ♦  (HB 2( L F221  +  L  )  -  F12l)*YRATI0 
FYZ1  •  FY11  ♦  (FY21  -  FY11)*ZRATI0 
F 112  -  HB2UF112  ♦  L> 

FY12  •  F112  ♦  (HB2(LF212*L)  -  F112)*YRATI0 
F122  -  HB2 (L  F122  +  L) 

FY22  «  FI 22  ♦  ( HB2( LF2224L)  -  F122)*YRATI0 
FYZ2  -  FY12  *  ( F Y 22  -  FY12)*ZRATI0 
H(U  •  FYZ1  ♦  ( FYZ2  -  FYZ1)*TRATI0 
50  CONTINUE 
U  >  H  ( 1 ) 

V  •  H  (  2 ) 

W  «  H(3) 

E  •  H  ( 4) 

RHO  -  H(5I 
RETURN 
C 

100  WRITE (6*105) 

STOP'  ABORT*  RBND903.  NO  DATA  ON  FILE.' 

105  FORHATI//'  ABORT*  RBND9D3.  NO  DATA  ON  FILE.') 

110  WRITE ( 6*115) 

STOP'  ABORT  RBN09D3.  FIRST  FILE  WORD  NOT  555.0* 

115  FORHAT(//'  ABORT*  RBND9D3.  HE  ADI  NOT  555.0') 

120  WRITE (6* 125) JN*  JNX»KN»KHX»NNHYPP*NNX 

STOP'  AB0RT*RBND903.  DIMENSION  TOO  SHALL.' 

125  FORMAT!//'  ABORT*  RBND9D3.  DIMENSION  TOO  SMALL.'/ 

♦  •  JN* JM  X*  KN»  KMX* NNH YPP#NHX» '*  618 ) 

130  WRITE (6*135) TH»  TSHIFT,T,TI1 

STOP*  ABORT, RBN09D3.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.' 

135  FORMAT ( •  ABORT*  RBN09D3.  REQUESTED  T  LT  INITIAL  TINE  IN  FILE.'/ 

♦  •  TH,TSHIFT,T,TI1-  **1P4E15.7) 

140  WRITE (6* 145) 

WRITE (6,GRI0) 

STOP'  ABORT*  RBN09D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.' 

145  FORMAT! //•  ABORT*  R6ND903.  SPATIAL  VARIABLE  OUTSIDE  RANGE.') 

160  WRITE(6,165)TH,TSHIFT,T»TIl,TI2 

STOP'  ABORT*  RBN09D3.  END  OF  FILE  BEFORE  T.' 

165  FORMAT! •  ABORT*  RBND9D3.  END  OF  FILE  BEFORE  T .'/ 

♦  •  TH»TSHIFT,T,TI1,TI2-  '*1P5E15.7J 
END 

•LABEL  ENDBN09 

•KEEPTO  END8N09  B BOUND 9  AND  DIMEN3 
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SUBROUTINE  B8ND9D 3( TH# XH* YH* ZH, U# V# W#E#  RHO > 

BOUNDARY  SUBROUTINE  FOR  BBOUND-9.  INRUT  RRE PARED  BY 
HULLUP  IS  ASSUMED  ON  FILE  T APE33  FOR  Z  -  EZIN). 

ZH  SHOULD  BE  ZO.  ZH  ♦  ZSHIFT  SHOULD  BE  ZIN. 

JN  ♦  2  RECORDS  PER  TIME  DUHPi 
HEADER— 335. 0#0LD  PR0B*TI1*TI2*  ZIN# NNH#IN#JN 
GRID— EXINEI)»I«1*IN)* EYINEJ J# J»l# JN) 

HYDRO  DATA  NNH/POINT  FOR  TINES  Til  AND  TI2 
JN  RONS  AT  EXIN(I)*Y1NEJ)*1«1#IN)  FOR  J»1#JN. 

CONN ON /SHI FT /TSHI FT#  XSHIFT#  YSHIFT#  ZSHIFT 
C  VALUES  IN  SHIFT  SET  IN  HULUN. 

C 

C  FOLLOWING  ARRAYS  SET  TO  ACCOMODATE  IN  .LE*  200. 

C  COMMON/ HB3/  TO  ACCOMODATE  LEVEL  2  STATEMENT. 

COMNON/HB3/HB3EAOOO) 

DIMENSION  H( 5 ) 

DIMFNSION  XINE200)#YINE200> 

C 

LEVEL  2#  HB3 

NAME  LIST/HE AO  IN/ HEADl# OLPROB# T II. TI2*ZIN»NNH, IN# JN 
NAMF L IST/GRI D/XH* XS HI F T# X#XIN*YH,YSHIFT*V# YIN# ZH# ZSHIFT* Z# ZIN 
C  SET  FILE  NUMBER.  STATE  ARRAY  DIMENSIONS. 

DATA  I END#  NBND# IMX# JMX#NNX/33#  63#  200*200* 4000/ 

C  NMX  IS  ADDED  TO  LCM  IN  HULL  IN  EHULL9250). 

C  SET  TF1  TO  -1.0  FOR  INITIAL  VALUES. 

DATA  TF1/-1.0/ 

C 

C  SHIFT  DATA  FROM  HULL  COORDINATES  TO  BOUNDARY  COORDINATES. 

T  -  TH  ♦  TSHI FT 
X  ■  XH  ♦  XSHIFT 
Y  •  YH  ♦  YSHIFT 
Z  -  ZH  ♦  ZSHIFT 

C  IS  THIS  AFTER  THE  INITIAL  ENTRYT 

IF (TF 1  .GT.  -0. 9) GOTO  30 
TF1  -  1.0 

C  READ  FIRST  HEADING*  GRID*  AND  HYDRO  SET.  CHECJ. 

REWIND  I BNO 

READ(IBND)HEAD1*0LPR0B*TI1*TI2#ZIN#NNH*IN*JN 
IFIEOF(IBND)  .NE.  0  )G0T0  100 
TTEST  -  0.00001* ( TI2  -  Til) 

ZTEST  -  0. 00000001*ZIN 
WRITE  E  6#  5) 

9  FORMATE//*  INITIATE  BB0UND9  INPUT.  FIRST  HEADER  RECORD*') 
WRITE (6* HE  AD IN) 

1FE A9SEHEAD1  -  599.0)  .GT.  0.0D60T0  110 
NNHYPR  -  NNH* IN 
NNHYPP  •  NNH YPR*2 

IFE  IN  .GT.  IMX  .OR.  JN  .GT.  JMX  .OR.  NNHYPP  .GT.  NMX)GOTO  120 
IF ( T ATTEST  .LT.  TIUGOTO  130 
0LPRB1  -  OLPROB 

C  READ  GRIO  FOR  INPUT  BOUNDARY  PLANE. 

REAO( IBND)  (  XIN(I)» I -1* IN  )# E YINE J ) # J -l# JN > 

XTEST  •  ( X  IN  <  2 )  -  XIN(I)  )*0.00001 
YTEST  ■  EYINE 2)  -  YIN(l)  1*0.00001 
10  IFET*TTE$T  *  LE.  TI2)G0T0  20 
C  FIND  THE  NEXT  HEADER  RECORD. 

1 A  RE AD ( IBND ) HE  A 01#  OLPROB# Til# TI2* ZIN*  NNH* IN* JN 
IFEE3FE IBND)  .NE.  OJGOTO  160 

IFEOLPROB  .NE.  0LPRB1  .OR.  HEAD1  .NE.  555.0EG0T0  1* 

C  THIS  IS  A  NEW  HEADER  RECORD.  BYPASS  THE  GRID  RECORD. 

RE  AD E IBND) DUMMY 
GOTO  10 
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READ  HYDRO  RECORDS.  STORE  IN  FILE  NBND. 

20  REWIND  NBNO 
DO  22  Jal*  JN 

READ ( IBND)  (H83(L)»L»1*NNHYPP) 

WRITE (NBNO)  (H83(L)»L"1* NNHYPP) 

22  CONTINUE 
24  REWINO  NBNO 

READ  FIRST  2  ROWS  FRON  FILE  NBND. 

READ ( NBND)  (HB3(L)#tal»NNHYPP) 

READ(NBND)  ( HB3  (NNHYPPK)  *  L«l»  NNHYPP ) 

I  IN  ■  2 

JIN  •  2 

ISW  ■  2 

INDH1  -  0 

INDH2  -  NNHYPP 

TOENOH  •  1.0/ITI2  -  Til) 

30  IF(T-TTEST  .ST.  TI2)fi0T0  14 

IF(ABS(Z  -  ZIN)  .ST.  ZTEST (GOTO  140 

IF(XINU)  -  X  .ST.  XTEST  .OR.  X  -  XIN(IN)  ,ST.  XTEST)S0T0  140 
IF(YIN(1)  -  Y  .ST.  YTEST  .OR.  Y  -  YIN(JN)  .ST.  YTEST)60T0  1A0 
32  IF(Y-YTEST  ,LT.  YIN(JIN)  )S0T0  36 
1FI2SW  .EO.  1)S0T0  34 
READ (NBNO)  (HB3(L)#L-l*NNHYPP) 

ISW  -  1 

INDH1  ■  NNHYPP 
INDH2  •  0 
SOTO  36 

34  READ ( NBNO)  (HB3(NNHYPP«L)»L*1>NNNYPP) 

ISW  -  2 
INDH1  >  0 
INDH2  -  NNHYPP 
36  JIN  ■  JIN+1 

YDENON  -  1.0/(YIN(JIN)  -  YIN(JIN-D) 

I  IN  •  2 
SOTO  32 

38  IF(Y*YTEST  , LT«  YIN(JIN-l)  )60T0  24 
LOCATE  THE  CORRECT  COtUHN. 

IF( X+XTEST  ,LT.  XIN(IIN-l)  )  IIN  •  2 
40  1FU-XTEST  .LT.  XIN(IIN)  )60T0  42 
IIN  -  IXN+1 
SOTO  40 

NOW*  SET  ADORESSES  AND  DO  THE  INTERPOLATION. 

42  IREF  -  (IIN  -  2 ) *NNH 

XRATIO  •  (X  -  XIN(IIN-I)  ) / ( XI N( I IN  )  -  XIN(IIN-l)  ) 

TRATIO  «  (T  -  TIDPTOENON 
YRATIO  -  (Y  -  YIN (JIN-1)  )• YDENON 

THREE  INDICES  THAT  ARE  1  OR  2  BELOW  REPRESENT  LOWER  Oi. 
HISHER  VALUES  OF  X*  Y#  OP.  T#  RESPECTIVELY. 

LF111  -  IN0H1  ♦  IREF 
LF211  -  LF111  ♦  NNH 
LF112  •  LF HI  ♦  NNHYPR 
LF212  •  LF 112  ♦  NNH 
LF121  *  INDH2  ♦  IREF 
LF221  ■  LF121  ♦  NNH 
LF122  •  LF 121  ♦  NNHYPR 
LF222  -  LF122  ♦  NNH 
00  50  L*l»  5 
Fill  •  HB3 ( LF 111  ♦  LI 

FX11  -  Fill  ♦  (HB3(LF211*L)  -  F111)*XRATI0 
F121  •  HB3 ( LF121+L) 

FX21  •  F121  ♦  (HB3( LF221+L)  -  F121)*XRATI0 
FXY1  •  FX11  ♦  (FX21  -  FX11)*YR ATIO 
F112  •  HB3(LF112  ♦  LI 

FX12  •  F112  ♦  (HB3(LF212+L)  -  F112)*XRATI0 
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F122  -  H63 I LF122U) 

FX22  «  F122  ♦  <HB3(LF222*L)  -  F122)PXRATI0 
FXY2  •  FX12  ♦  ( FX22  -  FX12WRATI0 
H(l)  •  FXY1  ♦  (FXY2  -  FXYDPTRATIO 
50  CONTINUE 
U  «  H<1) 

V  •  H ( 2) 

W  «  H( 3) 

E  »  H  ( 4) 

RHO  »  H( 5 1 
RETURN 


100  WRITE (6*105) 

STOP'  ABORT#  BBN09D3.  NO  DATA  ON  FILE.' 

103  FORMAT!//'  ABORT#  BBND9D3.  NO  DATA  ON  FILE.') 

110  WRITE (6#115) 

STOP'  ABORT  B8ND9D3.  FIRST  FILE  WORD  NOT  555.0* 

115  FORMAT!//'  ABORT#  BBND9D3.  HE  4D1  NOT  555.0') 

120  WRITE (6# 125) IN# IMX# JN# JMX#NNHYPP»NMX 

STOP'  AB0RT#BBND9D3.  DIMENSION  TOO  SMALL.' 

125  FORMAT!//'  ABORT#  BBND9D3.  DIMENSION  TOO  SMALL.'/ 

♦  •  IN#IHX#JN»JNX#NNHYPP»NMX-'#6I8) 

130  WRITE(6»135ITH,  TSHIFT»T,TU 

STOP'  ABORT# BBND903.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.' 

135  FORMAT!'  ABORT#  BBND9D3.  REQUESTED  T  LT  INITIAL  TINE  IN  FILE.'/ 

♦  •  TH#TSHIFT#T»TU-  '#1P4E15.?> 

140  WRITE  16,145) 

WRITE I6#GRI0 ) 

STOP'  ABORT,  BBND903.  SPATIAL  VARIABLE  OUTSIDE  RANGE.' 

145  FORMAT!//'  ABORT,  BBND9D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.') 

160  WRITE(6#165)TH,T$HIFT,T#TI1#TI2 

STOP'  ABORT#  BBND9D3.  END  OF  FILE  BEFORE  T. • 

165  FORMAT!'  ABORT#  8BND9D3.  END  OF  FILE  BEFORE  T.'/ 

♦  •  TH,TSHIFT,T#TI1#TI2-  *#1P5E15.7) 

END 

♦LABEL  ENDBND9 

♦KEEPTO  EN0BND9  TB0UND9  AND  DIMEN3 

SUBROUTINE  TBND903! TH# XH# YH# ZH#U# V# W# E# RHO ) 

C  BOUNDARY  SUBROUTINE  FOR  TBOUND-9.  INPUT  PREPARED  BY 

C  HULLUP  IS  ASSUMED  ON  FILE  T APE 34  FOR  Z  -  IZIN). 

C  ZH  SHOULD  BE  Z!KHAX)-0.5PDZ!KNAX).  ZH  ♦  ZSHIFT  SHOULD  BE  ZIN. 

C 

C  JN  ♦  2  RECORDS  PER  TIME  DUMP: 

C  HEADER— -555. 0»0L0  PR0B»TI1 # T 12# Z IN# NNH# IN# JN 

C  GRID— !X  IN!  I)#I«1#IN)#  (YIN(J)»J*1»JN) 

C  HYDRO  DATA  NNH/POINT  FOR  TIMES  Til  AND  TI2 

C  JN  ROWS  AT  !XINII)#YIN(J)#I-1#IN>  FOR  J-1#JN. 

C 

COMHON/SHI FT /TSHI FT# XSHIFT# YSHIFT# ZSHIFT 
C  VALUES  IN  SHIFT  SET  IN  HULL  IN. 

C 

C  FOLLOWING  ARRAYS  SET  TO  ACCOMODATE  IN  .LE.  200. 

C  C0MM0N/HB4/  TO  ACCOMODATE  LEVEL  2  STATEMENT. 

COMMON /HB4/HB 4! 4000 ) 

DIMENSION  HIS) 

DIMENSION  XINI200)#YIN!200) 

C 

LEVEL  2#  HB4 

NAME  LI ST /HE  AO IN /HE ADI#  OLPROB# T I1#T 1 2#Z IN#NNH# IN#  JN 
NAME L IST/GRIO /XH, XSHIFT# X#XIN, YH# YSHIFT# Y# YIN# ZH# ZSHIFT# Z# ZIN 
C  SET  FILE  NUMBER.  STATE  ARRAY  DIMENSIONS. 

DATA  IBND#  NBND# IMX#  JMX#NMX/34#  64,200,200#  4000/ 

C  NMX  IS  ADDED  TO  LCM  IN  HULLIN  (HULL9250). 

C  SET  TF1  TO  -1.0  FOR  INITIAL  VALUES. 
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DATA  TF1/-1.0/ 

SHIFT  DATA  FROM  HULL  COORDINATES  TO  BOUNDARY  COORDINATES. 
T  •  TH  ♦  TSHIFT 
X  ■  XH  ♦  X  SHI  FT 
Y  •  YH  ♦  YSHIFT 
2  •  ZH  ♦  ZSHIFT 

IS  THIS  AFTER  THE  INITIAL  ENTRY? 

IF(TF1  .GT.  -0.9IG0T0  30 
TF1  -  1.0 

READ  FIRST  HEADING/  GRID/  AND  HYDRO  SET.  CHECJ. 
REWIND  IBND 

READ ( IBND)HEA01/ OLPROB/ TI 1/TI2.ZIN/ NNHz IN/ JN 

IF(EOF (IBND)  .HE.  0  )GOTO  100 

TTEST  •  0.00041* ( TI 2  -  Till 

ZTEST  -  0. 00000001* Z IN 

WRITE (6/51 

5  FORHAT ( / /•  INITIATE  TB0UND4  INPUT.  FIRST  HEADER  RECORD* •) 
WRITE (6/ HE AD IN) 

IF(ABS(HEAD1  -  S5S.0)  .GT.  O.ODGOTO  110 
NNHY°R  •  NNH* IN 
NNHYPP  -  NNHYPR*2 

IF( IN  .GT.  IHX  .OR.  JN  .GT.  JNX  .OR.  NNHYPP  .GT.  NHX)GOTO  120 
IF(T +TTEST  .LT.  TIDGOTO  130 
0LPRB1  -  OLPROB 

READ  GRID  FOR  INPUT  BOUNDARY  PLANE. 

READ ( IBND)  (XIN(I)/I>1/IN)/(YIN(J)/J-1/JN) 

XTEST  -  ( X IN ( 2)  -  XIN(1)  )*0. 00001 
YTEST  •  ( YIN( 2)  -  YIN(l)  )*0.00001 
10  IFCT^TTEST  .LE.  TI2)G0T0  20 

FIND  THE  NEXT  HEAOER  RECORD. 

14  READ ( IBND ) HE ADI* OLPROB/ T 1 1/ T 12/ Z IN/ NNH/ IN/ JN 
IF(EOFdBND)  .NE.  0)G0T0  160 

IFIOLPROB  .NE.  0LPRB1  .OR.  HE ADI  .NE.  555.0)G0T0  14 

THIS  IS  A  NEW  HEADER  RECORD.  BYPASS  THE  6RID  RECORD. 
READ( IBNDIDUNHY 
GOTO  10 

READ  HYDRO  RECOROS.  STORE  IN  FILE  NBND, 

20  REWIND  NBND 
00  22  J-l/JN 

READ) IBND)  ( HB4 ( L ) / L -1/ NNHYP* ) 

WRITE (NBND )  ( HB4( L) / L" 1/ NNHYPP ) 

22  CONTINUE 
24  REWIND  NBND 

RE  AO  FIRST  2  ROWS  FROH  FILE  NBND. 

READ(NBND)  ( H84 ( L )/ L-lz NNHYPP ) 

RE  AO  (NBND)  (HB4(NNHYPP-»L)/L"1/  NNHYPP) 

I  IN  •  2 

JIN  -  2 

ISW  -  2 

1NDH1  -  0 

INDH2  «  NNHYPP 

TDENON  •  1.0/ (TI2  -  Til) 

30  XF( T-TTEST  .GT.  TI2)GOTO  14 

IF(ABS(Z  -  ZIN)  .GT.  ZTEST)GOTO  140 

IF (X  IN ( 1 )  -  X  .GT.  XTEST  .OR.  X  -  XIN(IN)  .GT.  XTEST)60T0  140 
IFIYIN(l)  -  Y  .GT.  YTEST  .OR.  Y  -  YINCJN)  .6T.  YTEST)60T0  140 
32  IF( Y-YTEST  .LT.  YIN(JIN)  1G0T0  38 
IF(ISW  .EO.  1 (GOTO  34 
READ ( NBND)  (  HB4( L )/ L  «1/ NNHYPP) 

ISW  -  1 

INDH1  •  NNHYPP 
INDH2  •  0 
GOTO  36 
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34  READ ( NBNO )  t HB4 (NNHYPP*L ) ,L-1, NNHYP P ) 

!SW  -  2 
INDH1  •  0 
INPH2  •  NNHYPP 
36  JIN  -  JIN*1 

YDENON  •  1,0/  (YIN(JIN)  -  YIN(JIN-l)) 

I  IN  •  2 
GOTO  32 

38  IF(Y-YTEST  .LT.  YIN(JIN-l)  IGOTO  24 
C  LOCATE  THE  CORRECT  COLUMN. 

IF(X*XTEST  .LT.  XIN(IIN-I)  )  IIN  ■  2 
40  IF(X-XTEST  .LT.  XIN(IIN)  IGOTO  42 
IIN  •  IIN-1 
GOTO  40 

C  NOW.  SET  AOORESSES  AND  DO  THE  INTERPOLATION. 

42  IREF  ■  (IIN  -  2I«NNH 

XRATIO  ■  (X  -  X  IN ( I  IN-1 )  )/(XIN( IIN)  -  XIN(IIN-l)  ) 

TRATIO  •  (T  -  Til )*TDENOM 
YRATIO  •  (Y  -  YIN (JIN-1 )  )*YDENOM 
C  THREE  INOICES  THAT  ARE  1  OR  2  BELOW  REPRESENT  LOWER  OR 

C  HIGHER  VALUES  OF  X,  Y,  OR  T,  RESPECTIVELY. 

LF111  •  IN0H1  ♦  IREF 

LF211  •  LF111  ♦  NNH 

LF112  •  LF111  ♦  NNHYPR 

LF212  ■  LF 112  ♦  NNH 

LF121  •  IN0H2  ♦  IREF 

LF221  -  LF121  ♦  NNH 

LF122  •  LF 121  ♦  NNHYPR 

LF222  •  LF 122  ♦  NNH 

00  50  L»l/ 5 
Fill  -  HB4(LF111  ♦  LI 

FX11  •  Fill  ♦  (HB4( LF211+L)  -  FU1)*XRATI0 
F121  •  HB4( LF121+L) 

FX21  •  F121  ♦  (HB4(LF221U>  -  F121)*XRATI0 
FXY1  -  FX11  ♦  ( FX21  -  FX11IPYRATI0 
F112  •  HB4 ( LF112  ♦  LI 

FX12  -  F112  ♦  <HB4(LF212*LI  -  FU2)*XRATI0 
F122  -  HB4 (LF122+L) 

FX22  •  F122  ♦  (HB4(LF222U>  -  F122)*XRATI0 
FXY2  -  FX12  ♦  ( FX22  -  FX12I*YRATI0 
H(L)  -  FXY1  ♦  (FXY2  -  FXY1)*TRATI0 
50  CONTINUE 
U  -  H  ( 1 1 
V  -  H  ( 2 1 
W  -  HOI 
E  •  H  ( 4) 

RHO  •  H(5) 

RETURN 

C 

100  WRITE (6.103) 

STOP*  ABORT.  TBN0903.  NO  DATA  ON  FILE.' 

105  FORMA T(  // '  ABORT/  TBND9D3.  NO  DATA  ON  FILE.') 

110  MRITE (6.115) 

STOP'  ABORT  TBN09D3.  FIRST  FILE  WORD  NOT  333.0' 

115  FORMAT ( // '  ABORT.  TBN0903.  HE  ADI  NOT  335.0') 

120  WRITE(6.125)IN. IMX.JN.JMX.NNHYPP.NMX 

STOP*  ABORT. TBN0903.  DIMENSION  TOO  SMALL.' 

125  FORMATS/'  ABORT.  TBND9D3.  DIMENSION  TOO  SMALL.'/ 

♦  •  IN, IMX.JN.JMX.NNHYPP.NMX-', 618) 

130  WRITE (6, 135ITH, TSHIFT.T.TIl 

STOP'  ABORT, TBND903.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.' 

135  FORMAT ( •  ABORT,  TBND9D3.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.'/ 

♦  •  TH.TSHI FT.T.TIl-  '.1P4E15.7) 

140  WRITE ( 6. 145) 


oooo  onooonnonn 


WRITE (6* GRID ) 

STOP'  ABORT*  T8ND9D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE* ' 

IAS  FORMAT!//'  ABORT*  TBN09D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.') 
160  WRITE (6*165)TH*TSH1FT»T»TI1»TI2 

STOP*  ABORT*  TBND9D3.  END  OF  FILE  BEFORE  T.' 

163  FORMAT ( •  ABORT*  TBN0903.  END  OF  FILE  BEFORE  T.'/ 

♦  «  TH*TSHIFT»T,TI1*TI2«  '»1P5E15.7) 

END 

•LABEL  ENDBND9 

•KEEPTO  EN0BND9  AB0UN09  AND  DIHEN3 

SUBROUTINE  A8N09D3( TH*  XH*  YH*  ZH*  U*  V*  W*  E*  RHO ) 

BOUNDARY  SUBROUTINE  FOR  ABOUND-9.  INPUT  PREPARED  BY 
HULLUP  IS  ASSUMED  ON  FILE  TAPE35  FOR  Y  ■  (YIN). 

YH  SHOULD  BE  YO.  YH  ♦  YSHIFT  SHOULD  BE  YIN. 

KN  ♦  2  RECORDS  PER  TINE  DUMP' 

HEADER— 555.0* OLD  »R0B*TI1*TI2* YIN*NNH* IN*KN 
GRID— <XIN< I),I-1, IN), (ZIN(K)*K«1,KN) 

HYDRO  DATA  NNH/POI NT  FOR  TINES  Til  AND  TI2 
KN  ROWS  AT  <X1N(I)*ZIN(K)>I*1#1N)  FOR  K»1*KN. 

CONHON/SHIFT/TSHIFT, XSHIFT* YSHIFT* ZSHI FT 
VALUES  IN  SHIFT  SET  IN  HULL  IN. 

FOLLOWING  ARRAYS  SET  TO  ACCOMODATE  IN  »LE.  200. 

C0HM0N/HB5/  TO  ACCOMODATE  LEVEL  2  STATEMENT. 

COMMON /HB5/HB 5 ( 4000 ) 

DIMFNSION  H( 5  ) 

DIMENSION  XIN(200)*ZIN(200) 

C 

LEVEL  2*  HB5 

NAMEL IST/HE ADIN/HEAD1»QLPR0B»TI1»TI2,YIN,NNH» IN*KN 

nahelist/grio/xh,xshift*x*xin,yh,yshift,y,yin,zh*zshift,z,zin 

C  SET  FILE  NUMBER.  STATE  ARRAY  DIMENSIONS. 

OATA  IBND,NBND,IMX,KHX*NMX/35, 65, 200*200, 9000/ 

C  NMX  IS  ADDED  TO  LCH  IN  HULLIN  (HULL9250) • 

C  SET  TF1  TO  -1.0  FOR  INITIAL  VALUES. 

DATA  TF1/-1.0/ 

C 

C  SHIFT  DATA  FROM  HULL  COORDINATES  TO  BOUNDARY  COORDINATES. 

T  •  TH  ♦  TSHIFT 
X  ■  XH  ♦  XSHIFT 
Y  •  YH  ♦  YSHIFT 
Z  ■  ZH  ♦  ZSHI FT 

C  IS  THIS  AFTER  THE  INITIAL  ENTRY! 

IF(TF1  .GT.  -0. 9 ) GOTO  30 
TF1  -  1.0 

C  RE  AD  FIRST  HEADING*  GRID*  AND  HYDRO  SET.  CHECK. 

REWIND  IBND 

READ<IBNO)HEA01,OLPROB*TIi*TI2*YIN,NNH,IN,KN 
IF ( EOF ( IBND)  .NE.  0  ) GOTO  100 
TTEST  •  0.00001* (TI2  -Til) 

YTEST  -  0. 00000001* YIN 
WRITE  ( 6,  5) 

3  FORMAT (//*  INITIATE  AB0UND9  INPUT.  FIRST  HEADER  RECORDi') 
WRITE ( 6*HEADIN) 

IF(ABS(HEAD1  -  355.0)  .GT.  O.ODGOTO  110 
NNHYPR  -  NNH* IN 
NNHYPP  -  NNHYPR* 2 

IF(IN  .GT.  IMX  .OR.  KN  .GT.  KMX  .OR.  NNHYPP  .GT.  NMXIGOTO  120 
IFIT+TTEST  .LT.  TIDGOTO  130 
0LPRB1  -  OLPROB 

READ  GRID  FOR  INPUT  BOUNDARY  PLANE. 

READ ( IBND)  ( XIN ( I )* I»l, IN )* (ZIN(K ), K»1,KN ) 


C 


XTEST  •  ( X  IN  (  2)  -  X IN ( 1 )  )*0. 00001 
ZTEST  -  (ZIN(2)  -  ZIN(l)  )*0. 00001 
10  IF(T+TTEST  .LE.  TI2JG0T0  20 

FIND  THE  NEXT  HEADER  RECORD. 

14  READ ( IBND ) HE ADI » OLPROB. T I 1» TI 2 , Y IN. NNH# IN.  KN 
IF ( EOF (IBND)  .NE.  OIGOTO  160 

IF  (  OL  PROB  .NE.  0LPRB1  .OR.  HEAD1  .NE.  555.0)G0T0  14 

THIS  IS  A  NEW  HEADER  RECORD.  BYPASS  THE  GRID  RECORD. 
READ( IBND) DUNNY 
GOTO  10 

READ  HYDRO  RECORDS.  STORE  IN  FILE  NBND. 

20  REWIND  NBND 
DO  22  K-l.KN 

READ! IBND)  ( HB5 < L >. L -1/ NNHYPP ) 

WRITE ( NBND )  ( HB 5 ( U » L- 1. NNHYPP > 

22  CONTINUE 
24  REWIND  NBND 

RE  AD  FIRST  2  ROWS  F  RON  FILE  NBND. 

READ ( NBND )  ( H85 ( L )» L»l. NNHYPP } 

RE  AD (NBND >  ( H85 ( NNHYPP+L ). L -1. NNHYPP ) 

IIN  •  2 
KIN  •  2 
ISW  -  2 
INDH1  -  0 
INDH2  »  NNHYR® 

TDEN3M  ■  1.0/ (TI2  -  Til) 

30  IF ( T-TTEST  .GT.  TI2)G0T0  14 

IF( ABS ( Y  -  YIN)  .GT.  YTEST)GOTO  140 

IF(XIN(1)  -  X  .GT.  XTEST  .OR.  X  -  XIN(IN)  .GT.  XTEST)GOTO  140 
IF(ZIN(1)  -  Z  .GT.  ZTEST  .OR.  Z  -  ZIN(KN)  .GT.  ZTESTIGOTO  140 
32  IF( Z-ZTEST  .LT.  ZIN(KIN)  )60T0  3B 
IFdSW  .EQ.  1IG0T0  34 
READ(NBND)  ( HB5 ( L ). L* 1. NNHYP«M 
ISW  »  1 

INDH1  •  NNHYPp 
INDH2  «  0 
GOTO  36 

34  RE AD (NBND)  ( HB5( NNHYPP 4L ) >  L-l.  NNHYPP  ) 

ISW  -  2 
INDH1  -  0 
INDH2  •  NNHYPP 
36  KIN  «  KIN+1 

ZDENOM  -  1.0/(ZIN(KIN)  -  ZIN(KIN-l)) 

IIN  -  2 
GOTO  32 

38  IF ( Z*  ZTEST  .LT.  ZIN(KIN-l)  )GOTO  24 
LOCATE  THE  CORRECT  COLUMN. 

IF ( X+XTEST  .LT.  XIN(IIN-I)  )  IIN  •  2 
40  IF ( X-XTEST  .LT.  XIN(IIN)  )GOTO  42 
IIN  -  IIN*1 
GOTO  40 

NOW.  SET  ADDRESSES  AND  DO  THE  INTERPOLATION. 

42  IREF  •  (IIN  -  2 ) *NNH 

XRATIO  •  (X  -  XIN(IIN-I)  >/(XIN(IIN>  -  XIN(IIN-l)  ) 

TRATIO  •  (T  -  Til )*TDENOM 
ZRATIO  •  (Z  -  Z IN (K IN-1 )  )*ZDENON 

THREE  INDICES  THAT  ARE  1  OR  2  BELOW  REPRESENT  LOWER  OR 
HIGHER  VALUES  OF  X.  Z»  OR  T,  RESPECTIVELY. 

LF111  •  INDH1  ♦  IREF 

LF211  -  LF111  ♦  NNH 

LF112  •  LF111  ♦  NNHYPR 

LF212  •  LF 112  ♦  NNH 

LF121  •  INDH2  ♦  IREF 

LF221  •  LF 12 1  ♦  NNH 


oooo  onnoonnnno 


LF122  ■  LF121  ♦  NNHYPR 
LF222  •  LF 122  ♦  NNH 
00  50  l-li5 
Fill  •  HB5 ( IF  11 1  ♦  L) 

FXll  -  Fill  ♦  (HB5( LF211+L )  -  F111)PXRATI0 
F121  «  HB5 ( LF 121  +  L) 

FX21  «  F121  ♦  ( H8 5( l F221  +  L  )  -  F121)*XPATI0 
FXZ1  ■  FXll  ♦  (FX21  -  F  Xll )*ZR AT  10 
F112  •  HB5 ( L F 11 2  ♦  L) 

FX12  ■  F 112  ♦  (HB5(LF212*L>  -  F112)*XRATI0 
F122  ■  HB5UF122+L) 

FX22  *  F122  ♦  (HB5(LF222»L>  -  F122)*XRATI0 
FXZ2  -  FX12  ♦  ( F X 22  -  FX12)*ZRAT10 
H(L )  -  FXZ1  ♦  ( FXZ2  -  FXZ1  )*TRAT 10 
50  CONTINUE 
U  »  H  ( 1 ) 

V  ■  H  (  2 ) 

W  •  H  <  3) 

E  •  H  (  4 ) 

RHO  «  H(  5 } 

RETURN 

C 

100  WRIT  E ( 6, 105 ) 

STOP'  ABORT,  ABND9D3.  NO  DATA  ON  FILE.' 

105  FORMAT!//'  ABORT,  ABN09D3.  NO  DATA  ON  FILE.') 

110  WRITE (6,115) 

STOP'  ABORT  A8ND9D3.  FIRST  FILE  WORD  NOT  555.0' 

115  FORMAT!//'  ABORT,  ABND403.  HE  ADI  NOT  555.0') 

120  WRITE (6,125) IN, I  NX,  KN»KMX#NNHYPP, NM X 

STOP'  ABORT, ABND9D3.  DIMENSION  TOO  SHALL . 1 
125  FORMAT ( // •  ABORT,  ABND903.  DIMENSION  TOO  SMALL.'/ 

♦  •  IN»IHX»KN»KNX»NNHYPP»NMX«*,6I8) 

130  WRITE(6,135)TH,TSHIFT,T,TI1 

STOP'  ABORT, ABND903.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.' 

135  FORMAT ( '  ABORT,  ABN09D3.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.'/ 

♦  •  TH»T$HIFT»T»TI1*  '#1P4E15.7) 

140  WRITE (6,145) 

WRITE ( 6, GRID  ) 

STOP*  ABORT,  ABN09D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.' 

145  FORMAT ( //•  ABORT,  ABN0903.  SPATIAL  VARIABLE  OUTSIDE  RANGE.') 

160  WRITE (6,165)TH,TSHIFT,T,TI1,TI2 

STOP*  ABORT,  ABN0903.  END  OF  FILE  BEFORE  T.' 

165  FORMAT ( '  ABORT,  ABND9D3.  END  OF  FILE  BEFORE  T.'/ 

♦  •  TH,TSHIFT,T,TI1,TI2«  *,1P5E15.7) 

END 

♦LABEL  ENDBND9 

♦KEEPTO  ENDBND9  FB0UND9  AND  DIME N3 

SUBROUTINE  FBN0903(TH,XH,YH,ZH,U, V, V»E,RHO) 

BOUNDARY  SUBROUTINE  FOR  FBOUND-9.  INPUT  PREPARED  BY 
HULLUP  IS  ASSUMED  ON  FILE  T APE 36  FOR  Y  •  (YIN). 

YH  SHOULD  BE  Y ( JMAX ) -0 . 5PDY UMAX ) .  YH  ♦  YSHIFT  SHOULD  BE  YIN 

KN  ♦  2  RECORDS  PER  TIME  DUMP: 

HEADER  — 5 55.0, OLD  »R0B,TI1,TI2, Y IN, NNH, I N»KN 
GRID— (XIN(I)»I«1, IN), (ZIN(K),K-1,KN) 

HYDRO  DATA  NNH/POINT  FOR  TIMES  Til  AND  TI2 
KN  ROWS  AT  (XIN(I),ZIN(K),I-1,IN1  FOR  K»1,KN. 

COMMON/SHIFT/TSHIFT. XSHIF T, YSH IFT, Z SHI  FT 
VALUES  IN  SHIFT  SET  IN  HULL  IN. 

FOLLOWING  ARRAYS  SET  TO  ACCOMODATE  IN  .LE.  200. 

CQNM0N/HB6/  TO  ACCOMODATE  LEVEL  2  STATEMENT. 

COMMON /HB6/HB6( 4000 ) 


DIMENSION  H( 5 ) 

DIMENSION  XIN(200)#ZIN{200) 

LEVEL  2*  HB6 

NAMELIST/HEADIN/HEAD1#  OLPROB# TI1»TI2#YIN#  NNH# IN»KN 
NAMELIST/GRID/XH, XS  H I F  T»  X*  XIN#  YH, YSH IFT# Y» YIN# ZH» ZSHI F T» Z# Z I N 
SET  FILE  NUMBER.  STATE  ARRAY  DIMENSIONS. 

DATA  I  BSD#  NBND> IM X,  KMX# NMX/36# 66# 200#200» 4000/ 

NMX  IS  ADDED  TO  L CM  IN  HULLIN  (HULL9250). 

SET  TF1  TO  -1.0  FOR  INITIAL  VALUES. 

DATA  TF1/-1 . 0 / 

SHIFT  DATA  FROM  HULL  COORDINATES  TO  BOUNDARY  COORDINATES. 
T  «  TH  ♦  TSHIFT 

X  •  XH  ♦  XSHIFT 

Y  •  YH  ♦  Y SH I CT 

Z  •  ZH  ♦  ZSHIFT 

IS  THIS  AFTER  THE  INITIAL  ENTRYT 
IFITFl  .GT.  -0, 9 ) GOTO  30 
TF1  «  1.0 

READ  FIRST  HEADING#  GRID#  AND  HYDRO  SET.  CHECK. 
REWIND  I BND 

R€AD(  IBND) HE ADI# OLPR OB # TI 1»TI2 # YIN# NNH# IN# KN 

IF(EOFdBND)  .HE.  0  (GOTO  100 

TTEST  ■  0. 00001* ( TI 2  -  Till 

YTEST  *  0. 00000001*Y IN 

WRIT  E  (  6#  5 ) 

5  FORMAT ( / / •  INITIATE  FB0UND9  INPUT.  FIRST  HEADER  RECORD**) 
WRITE ( 6#  HE  AD  IN) 

IF(  ABS  (HE ADI  -  555.0)  .GT.  O.ODGOTO  110 
NNHYPR  «  NNH* IN 
NNHYPP  ■  NNH Y PR* 2 

I F( IN  .GT.  I M X  .HR.  KN  .GT.  KMX  .OR.  NNHYPP  .GT.  NMX)GOTO  123 
IF  (  T  +  TTEST  .LT.  TIDGOTO  130 
0LPRB1  -  OLPROB 

READ  GRID  FOR  INPUT  BOUNDARY  PLANE. 

READ  (  IBND  )  (  XIN  (  I  )>  1  *1#  IN  J#  (  ZI N  ( K  ) »  K*l»  KN  ) 

XTEST  «  ( X IN ( 2)  -  XIN(l)  )*0. 00001 
ZTEST  «  ( Z IN ( 2)  -  ZINt 1)  1*0.00001 
10  IF ( T*TTE5T  .LE.  TI2JG0T0  20 

FIND  THE  NEXT  HEADER  RECORD. 

1 A  RE  AD  (  IBND)HEAD1,0LPR08»TU,TI2»YIN#NNH#IN#KN 
IF(EOFdBND)  ,NE.  0)GOTO  160 

IF ( OL  PROB  .NE.  0LPRB1  .OR.  HEAD!  .NE.  555.0)G0T0  14 

THIS  IS  A  NEW  HEAOER  RECORO.  BYPASS  THE  GRID  RECORD. 
READ! IBND) DUMMY 
GOTO  10 

READ  HYDRO  RECORDS.  STORE  IN  FILE  NBND. 

20  REWIND  NBND 
DO  22  K«1#KN 

RE AD ( I BND )  < HB6 < L )> L -1# NNHYPP > 

WRITE (NBND)  ( H86 ( L ), L-l# NNHYPP ) 

22  CONTINUE 
24  REWIND  NBND 

READ  FIRST  2  ROWS  FROM  FILE  NBND. 

READ(NBND)  ( HB6 < L )# L -1# NNHYPP ) 

RE  AO ( NBND )  ( HB6 ( NNH YPP*L )# L-l# NNHYPP ) 

I  In  •  2 

KIN  -  2 

ISW  -  2 

INDH1  «  0 

INDH2  •  NNHYPP 

TDENOM  «  1.0MTI2  -  Til) 

30  I F ( T-TTEST  • GT.  TI2)G0T0  14 


112 


t 

V 


If ( A8S  (Y  -  YIN)  .GT.  YTESTIGOTO  140 

IF(XINd)  -  X  .ST.  XTEST  .OR.  X  -  XIN(IN)  .ST.  XTEST)GOTO  140 
IF(ZIN(1)  -  Z  .ST.  ZTEST  .OR.  Z  -  ZIN(KN)  .ST.  ZTESTIGOTO  140 
32  IF(Z-ZT£ST  . LT.  ZIN(KIN)  >SOTO  38 
IFJISW  .EO.  1 ) SOTO  34 
READ(NBNO)  ( H86( L I. L-l.NNHYPP ) 

ISW  «  1 

IN0H1  ■  MNHYPP 
IN0H2  »  0 
GOTO  36 

34  READ(NBND)  ( HB6( NNHYPP+L) *  L“l.  NNHYPP  J 
ISW  >  2 
INDH1  •  0 
INDH2  ■  NNHYPP 


I 


i 


36  KIN  ■  KIN+1 

ZDENON  ■  1.0/ (ZIN(KIN)  -  ZIN(KIN-l)) 

IIN  -  2 
GOTO  32 
38  IF(Z+ZTEST 
LOCATE 
IF  C X+XTEST 
40  IF ( X-XTEST 
IIN  -  IINd 
SOTO  40 

NOW*  SET  ADDRESSES  AND  DO  THE  INTERPOLATION. 
IREF  •  (IIN  -  2 ) *NNH 

XRATIO  -  (X  -  XIN(IIN-l)  »  /  €  X  Z  N  ( IIN )  -  XIN(IIN-l) 
TRATIO  -  (T  -  TI1JPT0EN0H 
ZRATIO  ■  (Z  -  ZIN(KIN-I)  )*ZDENOH 


i  LT.  ZINC KIN-1 )  ISOTO  24 
THE  CORRECT  COLUMN. 

•  LT.  X IN ( IIN-1 )  )  IIN  •  2 
i LT.  XIN(IIN)  (GOTO  42 


42 


THREE  INOICES  THAT  ARE 
HIGHER  VALUES  Of  t,  Z. 


LF111 
LF211 
LF112 
LF212 
LF121 
LF221 
LF122 
LF222 
DO  ?0 
Fill 
FX11 
F121 
FX21 
FXZ1 
F112 
FX12 
F122 
FX22 
FXZ2 
H(L) 

30  CONTINUE 
U  •  H (II 
V  -  Hf 2) 

W  -  H  (  3) 

E  -  H ( 4) 
RHO  -  H( 5 ) 
RETURN 


IREF 

NNH 

NNHYPR 

NNH 

IREF 

NNH 

NNHYPR 

NNH 


-  IN0H1 

•  LF111 

-  LF111 

-  LF 112 
«  INDH2 

•  LF121 

•  LF 121 

•  LF122 
L-1.5 

HB6(LF 111  ♦  LI 

Fill  ♦  ( HB6( LF211»L  I  -  F111)*XRATI0 
HBb( LF121+L) 

F121  ♦  (HB6(LF221+L)  -  F121>*XRAT10 
FX11  ♦  (FX21  -  FX11I*ZRATI0 
HB6ILF112  ♦  L) 

F 112  ♦  ( HB6( L  F212  +  L  I  -  F112J*XRATI0 
HB6(LF122*L> 

F122  ♦  (HB6( LF222+L)  -  F122)»XRATI0 
(FX22  -  FX12IPZRATI0 
( FXZ2  -  FXZ1)*TRATI0 


1  OR  2  BELOW  REPRESENT  LOWER  OR 
OR  T#  RESPECTIVELY. 


FX12 

FXZ1 


100  WRITE (6.105) 

STOP*  ABORT.  FBND9D3.  NO  DATA  ON  FILE.1 
105  FORHATU/'  ABORT.  FBND9D3.  NO  DATA  ON  FILE.1  I 
110  WRITE (6.115) 

STOP*  ABORT  FBND9D3.  FIRST  FILE  WORD  NOT  555.0* 


113 


115  FORMAT (//•  ABORT*  FBND9D3.  HEAOl  NOT  555.0*  ) 

120  WRITE (6*125) IN* IHX*KN*KNX»NNHYPP»NMX 

STOP*  ABORT*  FBND9D3.  DIMENSION  TOO  SHALL.* 

125  FORHATt//*  ABORT*  FBND903.  DIMENSION  TOO  SHALL.*/ 

♦  •  IN*IHX*KN*KHX*NNHYPP*NHX-**6IB) 

130  WRITE ( 6* 135) TH*TSHIFT*T»TI1 

STOP*  ABORT, FBND9D3.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE.* 
135  FORMAT ( •  ABORT,  FBN09D3.  REQUESTED  T  LT  INITIAL  TIME  IN  FILE 

♦  '  TH,TSHIFT,T»TI1-  •» 1P4E15.7) 

140  WRITE (6*145) 

WRITE (6*GRI0 ) 

STOP*  ABORT*  FBN09D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.' 

145  FORMAT ( //•  ABORT*  FBND9D3.  SPATIAL  VARIABLE  OUTSIDE  RANGE.*) 
160  WRITE(6,165)TH,TSHIFT,T»TI1,TI2 

STOP*  ABORT*  FBN09D3 •  END  OF  FILE  BEFORE  T.* 

165  FORMAT ( •  ABORT*  FBND903.  END  OF  FILE  BEFORE  T.*/ 

♦  •  TH,TSHIFT*  T*TI1*TI2»  *,1P5E15.7) 

END 

*1  ABEL  ENDBND9 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 
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